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A comprehensive overview of kaon decays is presented. The Standard Model predictions are 
discussed in detail, covering both the underlying short-distance electroweak dynamics and the 
important interplay of QCD at long distances. Chiral perturbation theory provides a universal 
framework for treating leptonic, semileptonic and nonleptonic decays including rare and radiative 
modes. All allowed decay modes with branching ratios of at least lO"'^^ are analyzed. Some decays 
with even smaller rates are also included. Decays that are strictly forbidden in the Standard Model 
are not considered in this review. The present experimental status and the prospects for future 
improvements are reviewed. 
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I. INTRODUCTION 



Kaon decays have played a key role in the shaping of 
the Standard M odel (SM) ( |Glashow[ |1961[ [Salam] |1968 



Weinberg 1967|), from the discovery of kaons (Rochester 
and Butler 1947) until today. Prominent examples are 



the introduction of strangeness ( Gell-Mann 
1952), parity violation {K — >■ 27r, Stt puzz. 



1953^ Pais 



1954 



Lee and Yang 



1956[), quark mixing (Cabibbo 1963 



e) ( Dalitz 



Kobayashi and Maskawa[ 1973|), the discovery of CP vi 



olation (Christenson et aZTl964), suppression of flavor- 



changing neutral currents (FCNC) and the GIM mecha- 
nism (Glashow et al. 1970). Moreover, kaon decays con- 



tinue to have an important impact on flavor dynamics in 
constraining physics beyond the SM. 

The aim of this review is a comprehensive survey of 
kaon decays allowed in the SM with branching ratios of at 
least 10~^^. Some decays are included with even smaller 
decay rates. We do not cover decays that are strictly 
forbidden in the SM such as lepton-number or lepton- 
flavor violating decays. 

Kaon decays involve an intricate interplay between 
weak, electromagnetic and strong interactions. A major 
theoretical challenge has to do with the intrinsically non- 
perturbative nature of the strong interactions in kaon 
physics. The last 25 years have seen the development 
of a systematic approach to low-energy hadron physics 
in the framework of chiral perturbation theory (we re- 
fer to Sec. III.B for references). The approach provides 



a systematic expansion of decay amplitudes in terms of 
momenta and meson masses. The momenta of particles 
in the final states are sufficiently small so that the kaon 
mass sets the scale for the quality of the expansion. The 
relevant dimensionless ratio is AI'^/{AtiFt^Y — 0.18. 

Looking in more detail into the predictions of chiral 
perturbation theory (CHPT) for kaon decays, one ob- 
serves a rather wide range, from cases where great preci- 
sion can be achieved [especially in (semi)leptonic decays] 
to processes where hadronic uncertainties remain the ma- 
jor obstacle for comparison with experiment (mainly in 
nonleptonic decays). However, an important feature of 
CHPT in general is that it parametrizes the intrinsic 
hadronic uncertainties by a number of parameters, the 
so-called low-energy constants (LECs). The quality of 
theoretical predictions hinges to a large extent on the 
available information on those LECs. Over the years, a 
lot of progress has been achieved in the theoretical under- 
standing and the phenomenological knowledge of LECs. 
In recent years, lattice simulations have made important 
contributions to this field. For the time being, la ttic e 
QCD gives access mainly to strong LECs (see Sec. HI). 



This is one of the reasons why much better precision can 
be achieved in semileptonic than in nonleptonic decays. 
Kaon decays form a fascinating chapter of particle 
physics in themselves but they also give access to fun- 
damental parameters of the SM such as the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements V^s and 
Vtd- Although we stay strictly within the SM for this 



review, kaon decays have of course the potential to catch 
glimpses of New Physics, especially in those cases where 
the SM makes precise predictions. The constraints on 
New Physics are most effective in combination with other 
information from the low-energy /high-intensity frontier 
that are out of scope for this review. However, the im- 
pact of such constraints depends crucially on the quality 
of SM predictions. A state-of-the-art survey of kaon de- 
cays in the SM is therefore expected to be very useful 
also for physics beyond the SM. 

We review the status of K decays at a time when both 
theory and experiment progress at a somewhat slower 
pace than ten years ago. Nevertheless, there is still con- 
siderable activity in the field. The experimental program 
for kaon decays concentrates on specific channels of par- 
ticular interest for physics beyond the SM but results for 
other channels can be expected as by-products. The main 
players in the nea r future on the experiment al side will be 
NA62 at C ERN (|Collazuoll[2009}, K°TO (|Nanjo[ |2009| ) 



and TREK (Kohl 2010 



KLOD (Bolotov et al. 



Paton et al. 2006) at J-PARC, 



KLOE-2 at DX^NE~(|Amelino-Camelia et al] |2010D 



2009| ) and OKA ( [Kurshetsov 



2009 ) at IHEP Protvino, and the proposed Project-X at 



Fermilab (Bryman and Tschirhart 2010) 



Kaon decays have been treated in several reviews and 
lecture notes during the past 20 years ( Artuso et al. 2008 



Buchholz and Renk 
D'Ambrosio et al 



2008 



Barker and KetteU', '2000^ 'Ba ttiston et aLi|1992[ [Bryman 
1989; Buchalla, 2001 
1996; iBuras et al 



Winstein and Wolfenstein 1993 ) 



tenberg and Valencia 1993 Ritchie and Wojcicki[ 1993 



1997 Buras 



19941 LiH 



The review is organized along the following lines. We 
start with a brief summary of the theoretical framework, 
including both the short-distance aspects and the low- 
energy realization in terms of CHPT. For the latter, we 
also recapitulate the present knowledge of the coupling 
constants in the chiral Lagrangians (LECs). As already 
emphasized and in contrast to most previous reviews of 
the field, we then discuss essentially all kaon decays al- 
lowed in the SM that have either already been measured 
or that may become accessible experimentally in the not 
too distant future. We divide the actual review of differ- 
ent channels into three parts: leptonic and semileptonic 
decays, dominant nonleptonic decays {K — >■ 27r, 37r) and, 
finally, rare and radiative decays. 

In Sec.|nj we recall the effective Lagrangians in the SM, 
for both semileptonic and nonleptonic decays, after inte- 
grating out the W and Z bosons and the heavy quarks 
{t,b,c). For nonleptonic transitions, the leading QCD 
corrections are summed up with the help of the operator 
product expansion (OPE) and the renormalization group. 
The second part of this section contains a brief introduc- 
tion to CHPT, concentrating on the various chiral La- 
grangians relevant for K decays. Estimates of LECs are 



presented in Sec. HI In addition to phenomenological de- 
terminations, the large- iVc limit of QCD provides a use- 
ful framework for theoretical estimates. Strong LECs are 
dominated by the exchange of meson resonances. The re- 



suiting numerical estimates are in general in good agree- 
ment with upcoming lattice determinations. Through 
the hadronization of short-distance operators appearing 
in the effective Lagrangian for nonleptonic transitions at 
the quark level, electroweak LECs can be expressed in 
terms of strong and electromagnetic couplings. 



In Sec. IV we review the current status of leptonic and 
semileptonic modes (including radiative channels). K(i2 
and Ki^ decay rates can be predicted with great accuracy, 
providing non-trivial tests of the SM and allowing for the 
extraction of Vus ■ Ku decays are used as sensitive probes 
of chiral dynamics in tttt scattering whereas Ki,^ decays 
are included for completeness only. 

The dominant nonleptonic decays K -^ 27r, Stt are 
treated in Sec. |V] The two-pion modes are used to ex- 
tract the leading-order (LO) nonleptonic LECs Gg, 6*27 
on the basis of a next-to-leading-order (NLO) calcula- 
tion. To determine the s-wave tttt phase shift difference 
Sf){Mx) — 52{Mx), inclusion of isospin- violating correc- 
tions is mandatory. This applies also to the SM predic- 
tion of the CP-violating ratio e'/e although the theoreti- 
cal precision still does not match the experimental accu- 
racy. For the three-pion modes, NLO corrections signifi- 
cantly improve the agreement between theory and exper- 
iment. In addition to the CP-violating decay Ks — >■ Stt", 
special attention has been given to CP-violating asym- 
metries in the linear Dalitz plot parameter for the three- 
pion decays of charged kaons. Somewhat unexpectedly, 
K —i' Stt decays with at least two 7r° in the final state 
allow for a precise extraction of s-wave tttt scattering 
lengths by analyzing the cusp near threshold. 



Rare and radiative decays are considered in Sec. VI We 
first summarize the status of the rare decays K]^ — > Tr^vi? 
and K^ — >■ ir^i/D that can be predicted with a precision 
surpassing any other FCNC process involving quarks. 
The modes K -^ 'K'kvv can also be predicted with good 
accuracy but experimental limits are still far above the 
rates expected in the SM. Most of the remaining radia- 
tive modes K -^ ^(*)-f(*), K ^ i+£-, K -^ 7r^+^-, 



K -^ 



7r77 



(*) and K 



TTTT^^*' are dominated by long- 
distance dynamics to be analyzed in the CHPT frame- 
work. All channels have been calculated to NLO but 
in many cases estimates of the dominant next-to-next- 
to-leading-order (NNLO) effects are also available. Al- 
though of minor phenomenological interest at present, 
the decays K^ -> 87, Kl — > jj£'^i~, Kl -^ l^^v, 
Ks -^ itf-1^2^2^ K° -^ n^TT^j, K° -^ TT°TT°jj and 
K -> 3ttj are included for completeness. Conclusions 
and an outlook are presented in Sec. |VII| Some one-loop 
functions are collected in the appendix. 



II. THEORETICAL FRAMEWORK 

A. Short-distance description 

The SM predicts strangeness-changing transitions with 
AS* = 1 via W exchange between two weak charged cur- 



rents. At the kaon mass scale, the heavy W boson can be 
integrated out and the interaction is described in terms 
of effective four-fermion operators. 

Semileptonic transitions are mediated by the effective 
Lagrangian 

+ h.c. (2.1) 

where Vj^ denot e s the element ij of the CKM m atrix 
( Cabibbo] 1963[ Kobayashi and Maskawa[ 1973) and 
Gf = 1.1663788 (7) x IQ-" GeV-^ ( [WebbeTe^ at, ,2011] ) 
is the Fermi constant as extracted from muon decay. The 
universal short-distance factor 



2a 



^EW — 1 

= 1.0223 ±0.0005 



vr V An/ 



M„ 



7r2 J 



(2.2) 



encodes electroweak corrections not included in Gp ( Sir- 



lin 1978 1982|) and small QCD effects (jMarciano and 
SirlinI |1993P 



W exchange between two quark currents generates the 
AS* = 1 four-quark operator 



Q2 = [s7''(l-75)w] [u7p(l-75)d] 



(2.3) 



mediating nonleptonic K decays. Gluonic corrections 
bring further AS' = 1 operators, which mix under renor- 
malization ( Altarelli and Maiani 1974| Gaillard and Lee 



1974a[ [Shifman et a/.[|1977||Vainshtein et aL[|i975| T 

Qi = [r7''(l-75)^^] [zi^4l-75)d"], 
Q3 = [57^(1 -75)d] Y. [97/.(l-75)g], 

q—u,d^s 

Q4 = [s-"7^(l-75)rf^] E [9Am(1- 75)91, 

q—u.d.s 

Qs = [sr{l-l5)d] Y. [97^(1+75)9], (2.4) 

q—u.d^s 

Qe = [s-"7^(l- 75)^1 E [9Am(1 + 75)91, 

q—u.d.s 

where a, j3 denote color indices and color-singlet currents 
are understood whenever color labels are not explicit 
[qVq = g"rg"). 

Owing to the presence of very different mass scales 
(Afjr < Mk •€. Mw)j the QCD corrections are ampli- 
fied by large logarithms. The short-distance logarithmic 
corrections can be summed up using the OPE and the 
renormalization group, all the way down from My/ to 
scales fj. < rric (Oilman and Wise 1979 1980). One gets 



in this way an effective Lagrangian, defined in the three- 
flavor theory (Buras 1998), 



£, 



AS=1 
off 



~V^,V:,YG,{f^)Qdf^), (2.5) 



which is a sum of local four-fermion operators Qi, con- 
structed with the light degrees of freedom (m < fi), 
modulated by Wilson coefficients Ci{fi), which are func- 
tions of the heavy masses {Mz, Mw, rrit, mb, rric > /i) and 
CKM parameters. The unitarity of the CKM matrix, 



A„ + A, + A, -0, X, = V^^V;,, (2.6) 

allows to write the Wilson coefhcients in the form 



Ci(M) = Ziip) + Tyiip), 



(2.7) 



where t = —Xt/X^- The CP-violating decay amplitudes 
are proportional to the components yi{n)- 

The Wilson coefficients are known at the NLO (Bu ras| 
eTaTl |1993a[ [19961 [T993bl |Ciuchini ^ToTI [19931 ilOMT 
1995). This includes all corrections of 0{a^t'^) and 
©(af+^t"), where i = In (M1/M2) refers to the logarithm 
of any ratio of heavy mass scales Mi, M2 > fJ,- Moreover, 
the full rrit/Mw dependence (at lowest order in as) has 
been taken into account. 

The combination Q^ = Q2 — Qi and the 'penguin' 
operators Qi, i — 3,4, 5, 6, induce pure A/ — 1/2 transi- 
tions and transform like (8^, 1^) under chiral SU(3)l x 
SU(3)r transformations in flavor space, while Q'--'^'^^ — 
2Q2 + 3Q1 — Q3 transforms like a {27l, Ifl) operator that 
induces both A/ = 1/2 and A/ = 3/2 transitions. 

The inclusion of virtual electromagnetic interactions 
brings in the additional electromagnetic penguin opera- 



tors (Bijnens and Wise 



Lusignoli 1989 Sharpe 



1984 Buras and Gerard 1987 



1987 



q—u.d,s 

Qs = I [rrii -To)df'] J2 ^1 [9%.(l + 75)g1, 

q^u,d.s 



Q9 



[S7^(l-75)d] J2 e, [97a. (1- 75)9], (2-' 



q—u.d,s 



q—u,d.s 



where Cq denotes the corresponding quark charges in 
units of e = ^/ina. Their Wilson coefficients get also 
higher-order electroweak contributions from Z-penguin 
and W^-box diagrams. Under the action of the chiral 
group SU(3)l X SU(3)r the operators Qr and Qg trans- 
form like combinations of (8l, In) and (8l, 8^) opera- 
tors, while Qq and Qiq transform like combinations of 
(8l, 1r) and (27^, 1,^). 



Three more operators need to be considered in pro- 



cesses with leptons in the final state (Buchalla et al. 




W 



,c,t 



" 1j 



< — ^aaaaaaJ — > 

w 



qk 



FIG. 1 Penguin and box topologies. 



1996[|Gilman and Wise[|1980[ ): 



ni 



Qjv = [s7^(l 

Ql2 = Q7A - [S7''(l 



75)d] E [^"7m^] ' 

75)d] E [^"7m75^] , (2.9) 

Qi3 = Q{Di^) - [s^t^{l-j^)d][Djf,{l-j5)'^]. 

The mixing between the K^ and its antiparticle is in- 
duced through box diagrams with two W exchanges 
(Gaillard and Lee 1974b I, which in the three- flavor the- 



ory generate the effective Lagrangian ( Gilman and Wise 
T983| 



-AS=2 
-off 



(-f2 t\/t2 
'^|^CAS=2(A^)gAS=2(/i), (2.10) 



with 



Qa5=2 = [S7''(l-75)d] [57^(1-75)^] 



(2.11) 



There are in addition long-distance contributions from 
two f-^g^^ insertions. 

The box and penguin topologies, shown in Fig. [T] in- 
troduce contributions proportional to Xq {q = u, c, t) 
from virtual up-type quark exchanges. The hierarchy of 
CKM factors is better visualized through the | Wolfenstein| 
(1983) parametrization: A^ ~ A = \Vus\, Ac ~ — A and 

At A^X^, where A = \Vcb\IX^. The sum of the three 

contributions cancels for equal up- type quark masses, ow- 
ing to the unitarity relation (2.6). The quark-mass de- 



pendence of the associated loop functions breaks the GIM 
cancellation and strongly enhances the relevance of the 
top contribution in short-distance dominated processes. 
This is especially true for CP-violating effects, which are 
proportional to ( [Buras et al.\ 1994b) 

Ini(At) = -lm(Ac) ~ r^X^ A^ . (2.12) 

The corresponding CP-conserving factors are given to a 
very good approximation by 

Rc(A,) = -A (^1 - ^ j , 

Rc(A,) = -A'X'' (l - ^ Vl - p). (2.13) 



The value of A can be extracted from Kg^ decays and 
is given in Eq. (4.58). The remaining CKM parameters 



are obtained from global SM fits to flavor-changing data 
(jBona ei al.\ |2006| |Charles et al.\ 120051). We will adopt 



the values (Charles et al. 2010) 



A = 0.812t°:°^?, 



p = 0.144 ±0.025, 



1- — ) r/ = 0.342 ±0.016. 



(2.14) 



In order to fully calculate the kaon decay amplitudes, 
we also need to know the matrix elements of the operators 
Qi between the initial and final states, which involve non- 
perturbative dynamics at low energies. These hadronic 
matrix elements should cancel the renormalization-scalc 
dependence of the Wilson coefficients Ci (/i) . Methods to 
tackle this cha llenging task include lattice gauge theory 
( Boucaud et a/.^ 2005; Boyle and Garron[ 2010 Colangelo 



et al. 2011[ Goode and Lightman, 2010; Kim and Sachra- 



]da||2010[pLaiho and Van de Water 2010[pu||2010D, the 



1/JVc expansion (Bardeen et al 



bye et al. 2000 1998 1999), QCD sum rules (Guberina 



et aH|1985||Jamin and Pichl [1994 



and de Rafael 1985 1987 



1996 



1987a|bt |1988[ |Ham 



|Pich et a/.l|1986[|Pich 



Prades et al. 1991), 



functional bosonization techniques (jBruno anc 
|2004a| 



1993 Friot et al. 



Hambye etd] 2003 



etal] |1999b[ |2001| [Peris and de Rafaelj |2000| 



de Rafael 1991 ) and dynamical models ( Antonelli et al. 



Prades 



Knecht 



'ich and 



1996a|b||i997|[Bertolini et aL[ [19951 [1998[[Bijnens etal. 



2001[|2d66[|Bijnens and Prades[|1995a|b[|1999[|2000a|b| . 
CHPT provides a solid effective field theory framework 
to analyze systematically the long-distance dynamics, 
which we discuss next. Nevertheless, non-perturbative 
techniques remain necessary for matching the effective 
low-energy theory with the underlying QCD Lagrangian. 



B. Chiral perturbation theory 

Most kaon decays are governed by physics at long dis- 
tances. Below the resonance region one can use symme- 
try considerations to define another effective field the- 
ory in terms of the QCD Goldstone bosons. CHPT de- 



scribes (Gasser and Leutwyler 1985a Weinberg 1979) 



the pseudoscalar-octet dynamics through a perturbative 
expansion in powers of momenta and quark masses over 
the chiral symmetry-breaking scale A^ ^ 1 GeV. Chiral 
symmetry fixes the allowed operators, providing a com- 
prehensive framework for both semileptonic and nonlep- 
tonic kaon decays including radiative corrections. For 



general reviews of CHPT, see Ecker (1995), Pich (19951, 



de Rafael (^1995), Pich (1998), Scherer (2003), Gasser 



( |2004[ ) and |Bijnens| ( |2007D. The standard textb ook on 



the SM at low energies is Donoghue et al. ( 1992 ) 



The effective mesonic chiral Lagrangian in use today 
is listed schematically in Table |l] The strong chiral La- 
grangian relevant for (semi)leptonic decays to NNLO is 



TABLE I The AS = 1 chiral Lagrangian describing 
(semi)leptonic and nonleptonic kaon decays. Leptons 
are incorporated for radiative corrections in (semi)leptonic 
decays. The numbers in brackets denote the number of LECs. 



^-chir 



(# of LECs) 



loop order 



£,2(2) + £°f (0) + /:gf=.^(i) + c^'^=l{i) 



rr^xis) + <T;r(5) 

£^6(90) 



L = 



L< 1 



L<2 



iven by 




-^strong - 


= ^{D^UD^U^+XU^+X^U) 




^p2(2) 




+ Y^L.QP + £wzw 




^ * . r°dd|-ni 




V^io) 




i i 



(2.15) 



-Cp6(90) 



£°dd(23) 



Fq is the pion decay constant in the limit of chiral SU(3), 
the SU(3) matrix field U contains the pseudoscalar fields, 
D^ is the covariant derivative in the presence of exter- 
nal vector and axial-vector fields, the scalar field x ac- 
counts for explicit chiral symmetry breaking through the 
quark masses rriu, md, rrig, and (. . .) stands for the three- 
dimensional flavor trace. Explicit forms of the higher- 



nens et al. (1999 2002) and Ebertshauser et al. (2002). 



order L a grangians can be found in Gasser and Leutwyler 
([1985a|), [Wess and Zurnino] (|1971[), [Witten[ ([1983|), |Bij 



The nonleptonic weak Lagrangian to NLO is given by 



/2AS = lf-|-l 
^G8p2 U) 



G27FQ ( ^^23^11 + 0^^21-^13 



(2.16) 



GsF^ Y^ N,Of + G27F^ Y ^^Of + h.c. 



£-^--(22) 



£^'■' = 14(28) 

G27P 



The matrix L^ = iW D^U represents the octet oiV — A 
currents to lowest order in derivatives; A = (Ag — iA7)/2 
projects onto the s ~> d transition. The terms propor- 
tional to Gg and G27 transform under chiral transforma- 
tions like (Sijlfl) and (27i,,lfl), respectively, providing 
the most general effective realization of the corresponding 
short-distance operators in E g. (|2.5|. T he Lagrangians 
can be found in Cronin (19671, Kambor et al. (1990) and 
Ecker et aL| ( |1993[ ) 

To include electromagnetic corrections for both 
(semi)leptonic and nonleptonic decays, we also need the 
chiral Lagrangians 



^Icpton — •^kin (^) 



^FiY.^,&r 



)ton 



(2.17) 



the generating functional. The matrix elements and 
form factors with at most six external particles (an 
external photon counts as two particles) , at most one W 
(semileptonic decays) and with at most three propaga- 
tors in the one-loop amplitudes were presented in closed 
form by Unterdorfer and Ecker ( 2005 1 . This includes 
most of the processes covered in this review. Except 
for electromagnetic and anomalous contributions, all 
amplitudes to 0{p^) can be obtained from a Mathemat- 
ica program written by Rene Unterdorfer. The code 
along with several examples can be downloaded from 
http : / /homepage . univie. ac.at/Gerhard.Ecker/CPT- | 
amp. html, 



-,2„2 



(5) 



-'clm 



e'ZF^iQU^QU) + e^Gggcwki^t? (AC/^gC/) 



^Xo(l) 



/•cmwcak ( -y \ 



-t- e 



^nY^^fif^ 



•^"2' 2 (13) 



'G^F^Y. ^^Of^ + h-c, (2.18) 



Gge2p2V I 

where Q is the diagonal matrix of quark charges. The 
corresp o nding Lagrangians can be found in Ecker et~a[\ 
(|1989b i),iBijne ns and Wise] (| 1984i),iG rinsteine^aZ.|(|1986D, 
Urech, (1995), lEcker et al. (2000a) and IKnecht ei al. 



III. ESTIMATES OF LOW-ENERGY CONSTANTS 

A first-principle calculation of LECs requires to per- 
form the matching between CHPT and the underlying 
SM. This is a very difficult task. In many cases one re- 
sorts to phenomenology to determine the values of the 
LECs. For instance, most of the O(p^) couplings of the 
strong chiral Lagrangian, L^, are rather well known (Bij 



nens and Jemos 2011[ Ecker 2007p from low-energy data 
(tttt scattering, tt and K decay constants and masses, tt 
electromagnetic radius, tt — >■ 61^7, r decay). The electro- 
magnetic LEG Z arising at ©(e^p^) can be expressed in 
terms of the squared-mass difference of the pious: 



Ml± 



M^o = SnaZF^. 



(3.1) 



The matching procedure is of course simpler whenever 
the hadronic dynamics gets reduced to quark currents. 
For instance, in the GHPT language, the short-distance 
enhancement of semileptonic decays is encoded in the 



( |2000| )- The terms proportional to gcwk and Zi provide factor (|Descotes-Genon and Moussallam 



the low-energy realization of the electromagnetic penguin 
operators ( 2.8 ), while those with couplings Xj account for 



et al. . 2000 1 



2005 Knecht 



the operators in Eq. (2.9 1 with explicit lepton fields. 

The chiral realization of the AS — 2 effective La- 
grangian contains a single 0{p'^) operator: 



-AS=2 
-off 



G^M^ 



(4^) 



g^s^.F^ {XU^D^U){XU^D^U). 



l--(X6^-4Xy^l-y^6"'"^ 
which is related to S'ew in the Lagrangian (|2.1|) by 



2 yphys 



The (p'^) operators were discus sed by Kambor et al. 
(pool) andlEsposito-Faresel {|1991|). 



(2T91 



(3.2) 



(Mp) = ^EW - 1 + e'xf^'iMp), (3.3) 



where X^ ^^{Mp) denotes a residual long-distanc e con- 
tribution ( Descotes-Genon and Moussallam 2005). 



The GHPT framework determines the K decay am- 
plitudes in terms of the LEGs multiplying the relevant 
operators in the chiral Lagrangian. These LEGs encode 
all information about the short-distance dynamics, while 
the chiral operators yield the most general form of the 
low-energy amplitudes compatible with chiral symmetry. 
Ghiral loops generate non-polynomial contributions, with 
logarithms and threshold factors as required by unitarity. 

Practically all kaon decays discussed in this review 
have been calculated at least up to 0{p'^). The strong 
and nonleptonic parts of the corresponding NLO am- 
plitudes can be given in compact form in terms of 



A. The large-Nc limit of QCD 

The limit of an infinite number of quark colors is a 
very useful starting point to understand many features 
of QCD ('t Hooftl |1974a|b| |Witten[ [19791 ). Assuming 
confinement, the strong dynamics at Nc — >■ 00 is given 
by tree diagrams with infinite sums of hadron exchanges, 
which correspond to the tree approximation for some lo- 
cal effective Lagrangian. Hadronic loops generate correc- 
tions suppressed by factors of 1/Nc- Resonance chiral 



theory ( Cirigliano et al. 2006 Ecker et al. 1989a|b I pro- 



vides an appropriate framework to incorporate the mas- 
sive mesonic states ( Pich 2002 ) . Integrating out the res- 



onance fields, one recovers the usual CHPT Lagrangian 
with explicit values for the LECs, parametrized in terms 
of resonance masses and couplings. The resonance chi- 
ral theory generates Green functions that interpolate be- 
tween QCD and CHPT. Analyzing these Green functions, 
both for large and small momenta, one gets QCD con- 
straints on the resonance couplings and, therefore, infor- 
mation on the LECs. 

Truncating the infinite tower of meson resonances to 
the lowest states with '", O"*""*", 1 and 1+"'" quan- 
tum numbers, one gets a very successful prediction of the 
0{p'^Nc) strong CHPT couplings in terms of only three 
parameters: My, Ms and Fq. This provides a theoreti- 
cal understanding of the role of resonance saturation in 
low-energy pheno menology, which was recently extended 
to 0{p^) ( jCirighano et al. 2006[ ). Of particular inter- 
est for our present purposes are the resulting predictions 
for the 0{p^) couplings C12 and C34, which govern the 
amount of SU(3) breaking in the Ki^ form factor at zero 



momentum transfer ( Cirigliano et al. 2005 2006 1 



C12 
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8M|' 
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16M| 
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Ml 



(3.4) 

The large- A'^c limit turns out to be very useful to an- 
alyze the nonleptonic weak Lagrangian, because the T- 
product of two color-singlet quark currents factorizes: 



{J-J) ^ (J) (J) \i + o 



NcJ 



(3.5) 



Since quark currents have well-known CHPT realizations, 
the hadronization of the short-distance operators Qi can 
then be done in a straightforward way. As a result, in 
this large- A'c framework the electroweak chiral couplings 
can be related to strong and electromagnetic LECs of 
order p^, p'^, p^ and e^p^, respectively. The lowest-order 
electroweak LECs take the following values at large Nc 



(Cirighano et al. 2004a Pallante et al. 2001): 



-^ Ci(m) + I C2(m) + C^ip) - l6L,B{fi) Ceifi), 



9T7 = f[^i(m)+C2(a^)], 
o 

{e^gsge^kr = -3B(M)C8(/i) 
16 

and 



(3.6) 



gZ.. = CAs=2i^i), (3.7) 

where the dimensionless couplings gs , 527 sue defined as 

Gp ^ 



G8,27 - -—p^ViidYus ff8,27- 



(3.8) 



The operators Qi {i ^ 6,8) factorize into products 
of left- and right-handed vector currents, which are 
renormalization-invariant quantities. Thus, the large- A^c 
factorization of these operators does not generate any 
scale dependence. The only anomalous dimensions that 
survive when A'c — ^ 00 are the ones corresponding to 
Qg and Qg (IBardeen et al\ |1987bl IBuras and Gerard 



1987). These operators factorize into color-singlet scalar 
and pseudoscalar currents, which are [i dependent. The 
CHPT evaluation of the scalar and pseudoscalar currents 
provides, of course, the right /i dependence, since only 
physical observables can be realized in the low-energy 
theory. What one actually finds is the chiral realization 
of the renormalization-invariant products rriq 9(1,75)9. 
This generates in Eq. (3.6 1 the factors 



B{p.) = 



Ml 



p2 



{ttIs +md){p)F.„ 
^(2L8-i5) + 8 



1 



8M2 



(3.9) 



2Ml 



p2 



Af2 

-^ (3L4 - 4^6; 



which exactly cancel the /i dependence of C6,8(m) ^-t large 
Nq. There remains a dependence at NLO. 

The large- A"c expressions imply the numerical values 

g^ = (1.13 ± 0.05^ ±0.08l,±0.05„J 

+ T(0.64±0.15^±0.20L,t°;?^^J, 

0.46 ±0.01,,, 



ff27 
(ffsSewk) 



-1.60 ± 0.86^ ± 0.25^.1 ^JL 



+ 9.1 



r 25.0 ± 4.5^ ±1.0/^^ 



(3.10) 



where r is the ratio of CKM matrix elements in Eq. (2.7) 
and the main sources of uncertainty are indicated. 

While the CP-even part of g^ gets contributions from 
all (8l,1_r) short-distance operators, its CP-odd com- 
ponent is completely dominated by the strong penguin 
contribution, proportional to ryg(/i). The CP-odd com- 
ponent of 5cwk is dominated by the electroweak pen- 
guin contribution, proportional to Tyg,{ii), while the CP- 
even part receives contributions of similar size from both 
strong (Qe) and electroweak (Qs) penguin operators. Ex- 
plicit predictions for the higher-order electroweak LECs 
can be found in Cirigliano et al. (2004a). 

It is important to stress that the large- A'c' limit is only 
applied to the matching between the three-flavor quark 
theory and CHPT. The evolution from the electroweak 
scale down to /i < rric has to be done without any un- 
necessary expansion in powers of I/Nq] otherwise, one 
would nuss large corrections of the form \n{M/m)/ Nq, 
with M 3> m two widely separated scales. Similarly, 
the long-distance rescattering of the final pious in the 
K decay generates large logarithmic corrections through 
chiral loops, which are of higher order in both the mo- 
mentum and 1/A"c expansions (Pallante et al. 2001 1. 



These next-to-leading contributions, which give rise to 
the large s-wave tttt strong phases, are rigorously incor- 
porated through the CHPT framework. 



B. Lattice determinations 

Lattice results for strong LO and NLO LECs, both for 
SU(2) and for SU(3), have recently been summarized by 
the Flavianet Lattice Averaging Group -FLAG- ( Golan- 



gelo et al. 2011 ). We refer to this comprehensive review 



for a detailed discussion. 



IV. LEPTONIC AND SEMILEPTONIC DECAYS 

Purely leptonic and semileptonic modes are among the 
theoretically cleanest K decays. Using GHPT and lattice 
QGD, Ki2 and Ki^^ decay rates can be predicted with 
high accuracy and provide non-trivial tests of the SM 
and its extensions. On the other hand, K^^ decays can 
be used as probes of chiral dynamics in tttt scattering. 
In this section we review the current status oi K ^ iv 
{K12), K -> -Kiv {K13), and K — > TTTrii' {K14) decays, as 
well as their radiative counterparts. 



A. K12 (and 1112) decays 

1. Electromagnetic corrections 

The discussion of electromagnetic contributions to K12 
decays serves as the simplest example of the treatment of 
electromagnetism in (semi)leptonic processes. We start 
with the parametrization of the inclusive P — >■ iv^(^) 
decay rate proposed by |GirigHano and Rosell ( 2007a|b ) 



(here P = K 



± ^± 



), 






F{mt/Mp) 



(4.1) 



1- 



3 Mp 
2^"Mp 



SP) 



M^ 



p ~(P) 1 ^'^l 
c\ In — ^ 



M2 



mi 



SP + cPim./Mp) 



F{x) can be found in Marciano and Sirlin (19931. The 



JfT 



structure-dependent coefficients c\ are independent of 
the lepton mass mg and start at C(e^p^) in GHPT. The 
other coefficients appear only at higher orders in the chi- 
ral expansion. The one- loop result (order e^p^) for c\ 
was given by Knecht et al. ( 2000 ) , 



„(^: 



= -E^{Mp 



■>-i{ 



2 In 



Ml 
Ml 

M2 






(4.3) 



wher e Z is the {e^p ^) electromagnetic coupling given in 
Eqs. (2.18) and (3.1). E^{Mp) is a certain linear combi- 
nation of LEGs appearing in the Lagrangians C^Tr,2 and 

/^lcpton_ 



E'^ 



An^ 



-Kl 



-/a 



20 



k: 



-XI - AX^ 
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AXl 



K'. 



20 



Xf""'] . (4.4) 



The results in Eq. (4.3) are a nice example demon- 



strating the power of effective field theory methods. We 
see how the electromagnetic corrections to 1112 and K12 
of 0{e^p^) are related. In particular, c]^ and ci con- 
tain the same combination of LEGs E^ . Taking the ratio 
^Ke2M f^-^e^i-y) ' ^^^ coupling constant E^ cancels and the 
remaining expression 



(K) (tt) Z Mj^ 

=1 -'' ^J^^'m 



(4.5) 



is uniquely determined in terms of measurable quantities 



(Knecht et al. 2000). 



Finally, we note that using the matching calculation 
of Descotes-Genon and Moussallam ( 2005 1 for the LEG 



combination E^{Mp), one obtains 



„M 



-2.4 ±0.5, 



„(^) 



-1.9±0.5, 



(4.6) 



which is a modified version of the expression given in 
Marciano and Sirlin (1993). The decay rate in the ab- 



sence of radiative corrections is given by 



.(0) 
Pe2 



A-K 



Mp mi 1 - 



m,1 
Ml 



(4.2) 



where 14 = Vud, Vk — Vus and F^±,Fk± denote the 
pseudoscalar decay constants in pure QGD including 
strong isospin breaking.^ The first term in curly brack- 
ets is the universal long-distance correction for a point- 
like meson. The explicit form of the one-loop function 



^ For a discussion of subtletie s involved in this separation, see 
Gasser and Zarnauskas I 2010 1. 



where the errors given here are based on naive power 
counting of unknown contributions arising at 0{e'^p'^). 



2. Extraction of Vus/Vnd 



As suggested by Marciano (2004), a determination of 



|Kis/Kid| can be obtained by combining the experimen- 
tal values for the decay rates K — )■ nv{j) and tt — >■ ^iv('^) 
with the lattice determination of Fk/ Fj^, which is cur- 
rently performed in the isospin limit (?7i„ = m^) of QGD. 
The relation to be used to extract |I4s/Kid| reads 

rK..(,, \Vus? Fl Mk± (1 - mjlMl^f 



■^t2{-,} 



Kd?F^ MM^-mjlMl^Y 

1 + '5em+'5su(2)), (4.7) 



where F-j^ and Fk denote the decay constants in the 
isospin hmit, i5em is the long-distance electromagnetic 
correction (the short-distance part cancels in the ratio) 
and the strong isospin-breaking correction i5su(2) is de- 
fined by 



Cirigliano and Neufeld (2011) obtained 



WusWl 



1/2 



K 



\Vud\F^ 



0.23922 (25) x 



r 



Tf2(T) 



0.2763 ±0.0005. 



Fl^ 



Fl 



F\ Fi 



(l + '5sU(2)) 



(4.8) 



The electromagnetic correction is given by (Cirigliano 
and Neufeldl |20lTl iKnecht et al\ l2000| 



(4.16) 

Finally, taking as reference value for Fk/F-^ the FLAG 
average in Eq. ( 4.121, the ratio of the CKM ma trix ele- 
ments is given by (Cirigliano and Neufeld[ 2011) 



(^EM — 



F{me/MK) - F{me/M^) + 



^^In^ 
4 M2 



IK, 



= 0.2316 ±0.0012. 



(4.17) 



-0.0069 ±0.0017. 



(4.9) 



3. The ratio R' 



/m 



The 25% uncertainty of the numerical value is an esti- 
mate of corrections arising to higher order in the chi- 
ral expansion. The correction parameter 5su(2) reads 



In a first systematic calculation to 0{e'^p^), the coeffi- 



cients c. 



(P) JP) JP) s(^) 



(Cirigliano and Neufeld 2011) 



andRosell (2007; 



2- 



Tafb ). 



C2 were determined by Cirigliano 



0SU(2) 

where 

e(2) _ ^ 
4i?' 



V3e(2) 



This allowed the determination of 
rPe2(^)/rp^2(^) (P = TT, K) with an un- 



{Fk/F^ 



(4.10) 



3{AnrF^ 
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Ml 
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(4.11) 



the ratios R , — ^ j- ,, , , ^ j- ,, , 

precedented theoretical accuracy. In the SM, the ratios 

( p) 
R^, are helicity suppressed as a consequence of the V—A 

structure of the charged currents, constituting sensitive 

probes of New Physics. The two-loop effective theory 

results were complemented by a large- A^p calculation of 

an associated counterterm and a summation of leading 

logarit hms a" ln"(?7ip/me) ( Marciano and Sirlin 1993) 

giving ( Cirigliano and Rosell^ |2007a|b| 



With the FLAG ( jColangelo et al.\ |2011[ ) averages of lat 



i?i';|, = (1.2352 ± 0.0001) X 10" 



tice calculations with A^^^ = 2 + 1 dynamical fermions, 

i? = 36.6 ±3.8, Jif/F^ = 1.193 ±0.006, (4.12) 
Cirigliano and Neufeld ( |2011 ) obtained 



R 



{K) 



(2.477 ±0.001) X 10" 



(4.18) 



?(^) 



In the case of i?S the uncertainty arising from matching 
was increased by a factor of four to account for higher- 
order chiral corrections of ©(e^p^). The central value 



fSU(2) 



0.0044 ± 0.0005 ± O.OOllhighcr orders, (4-13) Qf ^(y jg jj^ agreement with the results of previous cal 



'e//i 



where the uncertainty due to higher-order corrections in 
the chiral expansion was estimated to be at a level of 
25%. We note that the strong isospin-breaking correc- 
tion (5su(2) is of the same order of m agni tude as the elec- 
tromagnetic correction ^em in Eq. (4.91 and should not 



culations (Finkemeier 1996 Marciano and Sirlin |1993 l, 
pushing the theoretical uncertainty below lO""*. The dis- 
crepancy with a previous determination of i? , can be 



e/M 



traced back to inconsistencies in the analysis of Finke- 
(p96l). 



be neglected in the extraction of the ratio |Ks/Kid|- 

Combined with the measured values for the leptonic 
widths of the pion ( Nakamura et al. . 2010[) 



meier 



The above theoretical results are compatible with 



V2(7) 



38.408 ±0.007 (/zs)-^ 



and of the kaon ( Antonelli et al. 2010b|) 



K„ 



51.25 ±0.16 ifis)- 



(4.14) 



(4.15) 



current experimental measurements (Ambrosino et al. 
2009a| [Bntton et alj \T99^ [19921 [Czapek et al.[ |1993 
Lazzeroni et al. 2011) averaging to 



e/fi 



R 



7a' 



cxp 



(1.230 ±0.004) X lO"-*, 
(2.488 ±0.012) X 10"^ 



(4.19) 



• Based on the result s from|Alltoii et aL|l|2008|;|Aoki et a;.l||20 09|; 
lAubin et al.\ (|2008| |2004|); [Bazavov et aLTl2010| [ 2009); Bcane 



et al. 



et al. 



louch 



They provide a clean basis to detect or constrain non- 
standard physics in these channels by comparison with 
upcoming measurements, which will pus h the fractional 



2007t ; [Bernard ei g/.l (120071); IBl ossier et a;.| 1,2009 J; .Diirr 
2010 i;|Follana et al. (2008); Go ckeler eral] \2006\ ; |Lel- 



2009 ;|Mawhinncy|(|2009|);|Noaki et al. (20091 



uncertainty from 0.3% to 0.05% in R^ ,' (Pocanic et al. 



(tt) 



2009 



Sher et al. 
1[ |2009[ 



azuo 



e/^i 



2009 ) and from 0.5% to 0.4% in i?^f ^ 



e/fi 
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B. K12-1 

The radiative leptonic decays -^'^27, 



£^e,n (4.20) 



allow one to probe the low-energy structure of QCD and 
its anomalous couplings. The matrix element reads 



T = -teGpV:, e* {FkL^ - H^'^l,) 



(4.21) 



where e^ denotes the photon polarization vector and the 
other quantities are given by {W — p — q) 



L^' 



miuip^){l + 75) 



2p^ 



2p^' -t- > 



2p- q 2pi-q 

l^' = u{p,)r{l-l5)v{pi), 

H'"' = J^ \iV{W'^)e^"'''^qo,pp 

- A{W^) {q ■ Wg"" 



vipe), 
(4.22) 



V2Mk 



Wq") 



V{W'^) and A(Ty^) denote the vector and axial-vector'^ 
Lorentz-invariant form factors characterizing the general 
decomposition of the correlator of weak and electromag- 
netic currents, 



^/7 = s-ff,u, 



^M = «7m75U, 



V^ 



(2uu — dd — ss)/3, 



(4.23) 



between the vacuum and the one-kaon state at q^ = 0, 



i.e., for real photons ( [Bijnens et al. 1993): 



(4.24) 



d'xe^'^^Q\T{V^^{x){V:{Q)-AZm) \K+{p)). 



The exphcit relation of A{W'^) and ^(W^^) to the above 
correlator can be found in the appendix of |Bijnens et al.\ 
(1993). The amplitude (4.21) is the sum of two terms: 



the so-called "inner-bremsstrahlung" (IB) term (propor- 
tional to i^), uniquely determined in terms of the non- 
radiative amplitude, and the "structure-dependent" (SD) 
term proportional to H^^'^ . Correspondingly, the spin- 
averaged differential decay distribution can be decom- 
posed into three terms (IB, SD, and INT, the latter de- 
noting the interference of IB and SD amplitudes) : 



d^T _ 
dxdy 

AsD 

2 

AiNT 



AmhB{x,y) 



{V + Af /sD+ (x, y) + {V- Af /sD- (x, v) 



[{V + A) /iNT+ (x, y) + {V- A) /iNT- (x, y)] 



(4.25) 



^ Note that we are using here the normahzation of V and A 
adopted in Nakamura et al. 1 20101. 



Here 



AiNT 



Af 2 



p\*' us 

327r2 
nii 



Fk 

Mk 



A 



SD, 



'-Mh 
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the independent kinematical variables are 
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y 



2p-pe 
Ml ■ 



(4.26) 



(4.27) 



and the functions / ib, /sd±, and /int± can be found in 
Bijnens et al. ( 1993 1. Note that the inner-bremsstrahlung 
and interference terms in the rate (4.25) are propor- 
tional to the helicity suppression factor {me/Mx)'^, so 
that Ke2'Y is dominated by the structure-dependent con- 
tribution. On the other hand, -^^^27, while sensitive to in- 
terference and structure-dependent terms, is dominated 
by internal brcmsstrahlung. 

Internal brcmsstrahlung is fixed in terms of a single 
hadronic input, namely the kaon decay constant Fk. 
The form factors V{W'^) and A(W'^) can be calculated in 
CHPT, with y(W2) arising from the anomalou s sector. 
The first non-trivial contributions arise at O(p^) (Bijnens 
et al 



1993): 



A,-. 



4V2M, 
^0 



K 



(LS+iio)= 0.042, 



V2Mk „ ... 
Vr,^ — „ o„ = 0.096, 



(4.28) 



where the numerical values were obtained by using the 
central values Fq = _FV, F[^{p- = Mp) — 6.9 x lO^'^, and 

The contributions of 0{p^) to the vector form factor 
involve one-loop graphs in the anomalous sector as well 
as 0(p^) counterterms. They were worked out first by 



Ametller et al. ( 1993 ) and later confirmed by Geng et al. 
]2004). To this order, V{W'^) acquires a non-trivial de- 



pendence on W = Mi(l — x), which can be approxi- 
mated by the linear parametrization 



ViW"^) ^Vq l + A(l-a;) 



(4.29) 



Geng et al. (2004) use estimates of the LECs in the 



anomalous sector from vector meson dominance and the 
chiral quark model. Their results can be summarized as: 



(Vo) 



p4_|_p6 



^ 0.078 ±0.005, 



A = 0.3 ±0.1. 



(4.30) 
(4.31) 



Our estimate of the uncertainty was obtained by looking 

at the spread in results obtained within different models. 

The contributions of 0{p^) to A{W^) were worked out 

first by Geng et al. (2004). They involve two- loop graphs 



11 



with insertion of the lowest-order CHPT Lagrangian £2, 
one-loop graphs with one insertion of £4, and tree-level 
graphs from L^. |Geng tt al\ ( |2004[ ) find that A^S/V^) is 
essentially flat in W'^ and receives a correction of about 
20% in the normalization, with an uncertainty induced 
by the LECs at the percent level [the subscript indicates 
that we quote the ly-independent part of A(W^^)]: 



(Ao)p4+p6 = 0.034. 



(4.32) 



The OijP") corrections are of the expected size for both 
V and A form factors. 

Experimentally, the ^^"£27 mode provides the best con- 
straint on the combination Vq -|- ^0 (Nakamura et al. 



2010 ). The data ar e dominated by the 1484 ± 63 events 
collected by KLOE fAmbrosino et al. 2009a I in the range 



10 MeV < -B™^" < 250 MeV. A fit to the measured spec- 
trum leads to the results 



Vo + Ao = 0.125 ±0.007stat±0.001,y,t, 

A = 0.38±0.20stat±0.02syst, (4.33) 

which are in excellent agreement with the 0{p^) CHPT 
predictions [from the work of Geng et al. ( 2004 ) we ex- 
tract (Vb + Ao)p4 



^p6 — 0.112 ± 0.005, where the uncer- 
tainty captures the model dependence in the LECs en- 
tering the p^ calculation]. However, the Kf^2'y extraction 
of the combination Vo + Aq differs from the Ke2'y value 



(Adler et al. 2000): 



^0 + ^0 = 0.165 ± 0.007stat ± 0.011 



systi 



(4.34) 



assuming V + A = Vo + Aq, i.e., no dependence on W^. 
Until very recently, the combination Vq — Aq was only 
loosely constrained by K£2j measurements, the best de- 
termination, Vq — Aq = 0.077 ± 0.028 (consistent with 
CHJ'T), coming from K 
[2002[) . The situation has chan ged with the rec e nt re- 



C. Ki3 

The photon-inclusive decay rate for all four K — > TTiv 
modes {K — K^, K'^, K'^; £ — e, fi) can be written as 



Kr. 



128^r^ '^^wl^+ ^^>^ 



X /' 



(A.) ( 



^ + "em + ''SU(2) 



(4.37) 



The Clebsch-Gordan coefficient Ck differs for neutral 
and charged kaons {Ck = 1 for K^^ and Ck = l/-\/2 
for K^^), while I^^ is a phase space integral depending 
on slope and curvature of the form factors f^^{t), de- 
fined by the QCD matrix elements 



{TT{p.^)\sjf,u\K{pK)) = 
\2 



(4.38) 



f"(t) + (pK-p.)M/^"W, 



where t = {pk — Ptt)^ — [pi+PvY- As usual, the vector 
form factor of the K''^ decay at zero momentum transfer 



has been pulled out in Eq. (4.37). The strong isospin 



breaking correction is defined as 

'5s''u(2)-(/+^(0)//f'^"(0))'-l. 
The long-distance electromagnetic corrections 






Kl , 



5^1,{V,) + 5^{^{V,,,) 



(4.39) 



(4.40) 



receive contributions from three-particle and four- 
particle final states. In the following sections we review 
the theoretical quantities appearing in the expression of 
the decay rate and then combine this information with 
experimental input to extract the CKM element Vus- 



1. Electromagnetic effects in A'f3(^) decays 

The calculation of the electromagnetic contributions 
to ©(e^p^) in Ki^ decays using the methods of CHPT 



suit s from ISTRA-H reported in Tchikilev et al. (2010) 

and|DukeiaL|(l20lTl). They studied the if ^27 mode in a ^1*^ virtual pho tons and leptons (jKnecht et al.\ |2000D 



kinematic region in which the interference terms can be 
extracted, thus providing direct sensitivity to the com- 
bination Vq — Aq. Different analyses assuming constant 
form factors lead to (]Tchikilev et al. 2010) 



was presented in Cirigliano et al. (2002). Based on this 



analysis, full numerical results on the Kf.3 decay modes 
were given in Cirigliano et al. ( 2004b I , adopting a specific 



Vo-Aa = 0.126 ± 0.027,tat ± 0.047, 



systi 



and (Duk et al. 2011) 



Fo - Ao = 0.21 ± 0.04stat ± 0.04, 



:syst: 



(4.35) 



(4.36) 



prescription for treating real photon emission and a spe- 
cific factorization scheme for soft photons. This approach 
resulted in the partial inclusion of higher-order terms in 
the chiral expansion. In a recent publication (Cirigliano 



respectively, about 1.5 and 3 sigmas above the CHPT 
prediction Vq — Aq = 0.044 ± 0.010, where we take as 
uncertainty the difference between the central values at 
0{p^) and 0{p'^). At this stage it is clearly premature to 
claim a serious tension between data and theory. 



et al. 2008b I , the numerical analysis of electromagnetic 
corrections was extended to a complete study of K^^ de- 
cays. At the same time, previous results for the K,^^ 
modes were update d using the new estimates of elec- 
tromagnetic LECs ( Ananthanarayan and Moussallam 



2004 Descotes-Genon and Moussallam 2005 1 that af 



feet the structure-dependent electromagnetic contribu- 
tions. Rathe r than using the soft-p hoton factorization 
procedure of Cirigliano et al. (|2002|), the analysis was 



performed at fixed chiral order (rjr. 
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TABLE II Electromagnetic corrections to Kes, decay rates 
(Cirigliano et al. 2008b I. D3 refers to three- particle phase 



space, I?4 to four-body kinematics. 





6^Ml>3)i%) 


5i^(04-3)(%) 


5SU%) 


K% 


0.50 


0.49 


0.99 ± 0.22 


K% 


-0.35 


0.45 


0.10 ± 0.25 




1.38 


0.02 


1.40 ± 0.22 


^i 


0.007 


0.009 


0.016 ± 0.25 



Table |lT] summarizes the numerical results for the 
long-distance radiative corrections. Two characteristic 
features can be understood by qualitative arguments. 
Firstly, the electromagnetic corrections for the neutral K 
decays are expected to be positive and sizable on account 
of the final-state Coulomb interaction between £+ and 
w~ producing a correction factor of 7ra/w^+^_ ~ 2% 
over most of the Dalitz plot. While the exact correction 
and the relative size of K^^ and K^;^ depend on other 
effects such as the emission of real photons, the qual- 
itative expectation based on the Coulomb interaction 
is confirmed by the detailed calculation. Secondly, 
the large hierarchy (5gj^(I?4_3) <^ S^^CD^-^) admits 
a simple interpretation in terms of bremsstrahlung off 
the charged lepton in the final state. The probability of 
emitting soft photons is a function of the lepton velocity 
vi which becomes logarithmically singular as vi — > 1, 
thus enhancing the electron emission. For typical values 
of V£ in I?4_3, the semiclassical emission probability 
implies <5f^(p4-3)/<5EM(2'4-3) ^ 20 ^ 40. 

The theoretical uncertainties assigned to the 6^^^ in 
Table III] arise from two sources: the input parameters 
(LECs and form factor parameters) used in the calcula- 
tion and unknown higher-order terms in the chiral ex- 
pansion (the latter would require a complete analysis at 
order e^p'^). For a detailed discussion of the error esti- 
mate, we refer to 



The differentia 
the form 



Cirigliano et al. (2008b). 



decay distribution can be written in 



dydz 



Gp \Vus 



2 7i,,r5 /^2 

128^3 ^Ewl/+ (0)1 



p(o)(y,z) + <5pE^(y,z) , (4.41) 



where the Lorentz invariants y — 2pK-Pi/M^ = 2Ei/Mk 
and z = 2p7r • Pk/M]^ = 2E-^/Mk are related to the en- 
ergy of the charged lepton and of the pion, respectively, 
measured in the rest frame of the kaon. Here p^^'{y,z) 
represents the Dalitz plot density in absence of electro- 
magnetic corrections, while 6p^^{y,z) accounts for elec- 
tromagnetic effects. Illustrative figures displaying the 
relative size of the electromagnetic corrections over the 
Dalitz plot can be found in Cirigliano et al. (2008b). It 



turns out that the corrections to the Dalitz distributions 
can be locally quite large (~ 10%) and do not have a def- 
inite sign, implying cancellations in the integrated total 
electromagnetic corrections. 



2. Quark mass ratios and Ktj decays 

The isospin-breaking correction to K^r^ decays takes 
the form (Cirigliano et al. 2002) 



"SU(2) 



2V3(e 



- e 



(4) _L _(4) 
=-EM 



(4.42) 



It is dominated by the lowest-order tt^-t] mixing angle 
£*^^^ defined in Eq. (4.11). The NLO corrections e ^ (or- 
der p'^) and £pA (orde r e'^p ^) were computed in Gasser 
and Leutwyler ( |l985b ) and Neufeld and Rupertsberger 
'( 1995 ) , respectively. The explicit expressions for these 



quantities can be found in [Cirigliano et al. (2002). The 
dots refer to NNLO contributions (arising at order p^), 
for which the late st results can be found in |Bijnens and| 
Ghorbanil ( |2007b[ ). Working strictly to ©(p^JlnCHPT 

and neglecting the tiny contribution e^f^, one can relate 



-'SU(2) 



to quark mass ratios by 



''SU(2) 



3 1 

2Q2 



Ml 

Mi 



2 



1 



(4.43) 



where Q^ = {''^1 — ™^)/("T-d ~ ™u) = Rims/fh + l)/2, 
MK,Tr are the isospin-limit meson masses and Xp^ = 0.219 
is a calculable loop correction. Thus, Sg^,2) i^ essen- 
tially determined by Q^ (the uncertainty in nis/m affects 



negligibly due to the smallness of Xp-^ ) ■ 



"SU(2) 

The standard strategy up to now has been to use 
all known informa tion o n light quark masses to predict 
'^SXJ(2) using Eq. (4.43). The double ratio Q^ can be 
expressed in terms of meson masses and a purely electro- 
magnetic contribution (Gasser and Leutwyler 1985a), 



Aif^M|.(l + 0(m2)) 



M2[A 



where A; 



(A 



KOK+ 



(4.44) 
Mq. The electromagnetic term 
iro)EM vanishes to lowest order 



KOK + 



Ml 



(|Dashen[|1969[). A calculation at 0{e^ p^) gives (Neufeld 
and Rupertsberger 1995 Urech 1995[ ) 



eV 



{^K"K+ + A7r+7ro)EM — 6 M j^ 



lezLiif,) 



-H -(^5 + KeYi^i) - 8iK,o + K^.n^i) (4.45) 



^31n^-4 + 2Zlni^ 



47r2 



f^ 



M 



+ Oie'M^) 



The current estimates ( Ananthanarayan and Moussal- 



1am 2004) of the electromagnetic LECs appearing in 
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this expression imply a large deviation from Dashen's 



limit, (A;^o^+ + A^+^o)em 



-1.5 A. 



, which im- 
plies Q = 20.7 ± 1.2 ( |Kastner and Neufeld[ [2008| . Such 



a small value of Q [compared to Q = 22.7 ± 0.8 given by 
Leutwyler (1996) ] was also supported by o ther studies 
\Q = 22.0 ±0.6 inJBijnens and Pradesl ( |1997[ ) and Q ~ 20 

however, 



Amoros et al. 



( |200ip j. It should be noted 
that the rather large value Q = 23.2 was obtained from 



an analysis of ?/ -^ 3n at two loops (Bijnens and Ghor- 
bani 2007a I. On the other hand, the non-lattice deter- 
minations of the second input parameter ms/fh ^ 24 
have remained rather stable over the last years. Com- 
2 0.7 ± 1 .2 with iris/m 
((2008| found R = 33.5 ± 4.3 and finally 



bining Q 



and Neufeld 

S^^f^^ =0.058 ±0.008. 

It is worth stressing that the present precision of the 
decay rates and of the radiative corrections permits also 
an "experimental" determination of <^§uf2') ? which can be 
used as a constraint on the quark mass ratio Q via the 
formula (4.43). Combining recent Kg^ data (Antonelli 



et al. 2010b) with the expression 



''SU(2) 



2r 



^i 



'Kfe 



r, 



u 



M, 



KO 



M. 



1 



{si 



EM 



<5^ 



EM 



KU3 J-K+e \^^^K+ . 

(4.46) 
one obtains Sg^,2)exp ~ 0.054 ± 0.008, in perfect agree- 
ment with the value obtained from quark mass ratios. 

Alternatively, one may use the Nf = 2 + 1 lattice av- 
erage (Colangelo et al. 2011) nis/fh = 27.4 ± 0.4 being 



considerably larger than the values obtained with non- 
lattice methods. Combined with Q = 22.8 ±1.2 from 



>;K^ 



4.43|) yields .3^(2) 



the same data compilation, Eq. 
0.048 ±0.006, still consistent with the experimentally de- 
termined result. 



3. Form factors and phase space integrals 

Calculation of the phase space integrals Iki requires 
knowing the momentum dependence of the for m fact ors. 



The vector form factor f^'^{t) defined in Eq. (4.38) re- 



presents the p-wave projection of the crossed-channel ma- 
trix element (0|s7'^u|-?ir7r) whereas the s-wave projection 
is described by the scalar form factor 



/o(t)-/+W + 772 



t 



Ml - M2 



f-it)- 



(4.47) 



It is convenient to normalize all the form factors to 
/+" ^ (0) (denoted /+(0) in the following). In terms 
of the normalized form factors /j(i) = fi{t)/fj^{0), the 
phase space integrals read 



24.7±1.1, Kastner UTd. 



with A = (t - {Mk + M^f){t - {Mk - M^f). 

Traditionally, a polynomial parametrization has been 
used for the form factors, 



/+,o(t) = l + A' 






t 



AP 



(4.49) 



where A^ p ^-^d A" q are the slope and curvature, re- 
spectively. Fits to the experimental distributions of K^^ 
decays allow to extract the parameters A'^, A", and Aq. 
The resulting uncertainty on the phase space integrals is 
at the level of 0.12% for Ikc and 0.30% for Ik^ (jAntonelli 
2010bl) 



This affects the extraction of V^s at the 
level of 0.06% (Kes) and 0.15% (if^g). 

Other form factor parametrizations have been pro- 
posed, in which, by using physical inputs, specific rela- 
tions between the slope, the curvature and all the higher- 
order terms of the Taylor expansion (4.49) are imposed. 



This allows to reduce the correlations between the fit- 
ted slope parameters: only one parameter is fitted for 
each form factor. Explicit examples used to analyze data 
include the pole parametrization, dispersive parametriza- 



tions (Abouzaid ei aL 2010 Bernard ei aL 2006 20091, 



and the so-called z-parametrization (Hill 2006) 



4. The Ki3 scalar form factor 

SM predictions for the slope parameter Aq of the scalar 
form factor of Ki^ decays were obtained by using differ- 
ent approaches. In the isospin limit, the combination 



of a two- loop result in chiral perturbation theory (Bij 



nens and Talavera 2003j ) with an updated e stimate of 
the relevant p" l ow-energy couplings base d on Cirigliano 
et al. (2005) and Cirigliano et al. j 2003b ) gave the result 



AT 



(T3:9To;j^±0.4)xlO-^ (jKastner and Neufeldl 20081. 



Dispersive methods were employed by several authors.^ 
Typical numbers for the resulti ng scalar slope para meter 



are: A'q = (14.7 ±0 4) x IQ-^, (Jamin et al. 2006), and 



a;, = 13.71 X 10-3 ( [Bernard etal. \\20Ul 

The low-energy theorem of Callan and Treiman 



( Callan and Treiman 1966 Dashen and Weinstein 1969 ) 



predicts the size of the scalar Kg^ form factor at the (un- 
physical) momentum transfer t — Aktv, 



fo{AK^)=FK/F^ + AcT, 



(4.50) 



with a correction term of 0{mu,md,e'^). In the isospin 
limit (m„ = md, e = 0), and at first non-leading order, 
the tiny value Aqt = —3.5 x lO"'^ was computed by 
Gasser and Leutwyler| ( |1985b ) . A discussion of higher- 
order effects on this quantity can be found in Bijnens and 



Ixe 



*« dt 
.? Ml 

flit) 



A3/2 1 + 



2t 



^"'^^^^/^(O), 



{2t + m'j)X'' 



~2t 



(4.48) 



* See for instance Jamin et al. (2002), Jamin et aZ.|||2004l,|janiin 

I et al. (|2006|l, |Bernard et al., (2006^, |Bernard and Passemar 

(2008 . [Bernard et a;.| ( |2009| l, |Bernard et a;.| ( |2011[ |, [Abbas etHi 

mi6\. 
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TABLE III Experimental results for the slope parameter of 
the scalar Ka form factor, in units of 10^'^. 



Experiment 


a;, 


Ref. 


1STRA+ (k;^) 
KTeV (Kl^s) 

KTeV {Kl^s + KLe-i) 
NA48 (A-L^a) 
KLOE [Kl^^) 

KLOE {Kl^3 + i^t.a) 


17.1 ±2.2 
12.8±1.8 
13.7±1.3 
9.5 ±1.4 
9.1 ±6.5 
15.4±2.2 




Yushchenko et al. 


20041 


Alexopoulos et al. 


20041 


Alexopoulos et al. 


(20041 




Lai et al. 


20071 


Ambrosino et al. 


20071 


Ambrosino et al. 


20071 



Talavera (2003) and Kastner and Neufeld (2008). Note 



that the constraint at the Callan-Treiman point plays an 
essential role in the dispersive analysis of the scalar form 



Jamin et al. 2004 2006). 



factor ( Bernard et al. 2006 Bernard and Passemar 2008 



The effect of isospin violation and electromagnetic cor- 



rections has also been considered (Kastner and Neufeld 
2008). This introduces an additional uncertainty for 



the values of the slope parameters of at most ±10"'^, 
mainly due to not yet fully determined low-energy cou- 
plings. Combining the loop results given by |Bijnens and| 
Ghorbani ( 2007b ) with an estimate of the relevant low- 



energy couplings, the difference of the slope parameters 



of the neutral and charged kaons was found ( Kastner and 



Neufeldl 



2008) to be confined to the rather small range 



KiK^s) ^ 



10- 



The present experimental situation is displayed in Ta- 
ble lllll Note that the numbers^ shown here are those 



where the quadratic parametrization (4.49 ) has been used 



for the simultaneous determination of the vector form 
factor . The K TeV data were also reanalyzed by |Abouzaid| 
et al. ( |2010 ) using the dispersive parametrization. The 
' ^ (13.22 ± 1.39) X 10-3 (KLt,3 only) and 



results Ag 
A^ 







(12.95 ± 1.17) X 10-3 (X 



Lfi3 



K 



Le3> 



are rather 



close to the corresponding numbers given i n Table |III| 
The a nalogous procedure for the KLOE data (Ambrosino 



etal^2W7] gave AJ, = (14.0±2.1) x 10-^ (i^L^a+i^LeSj- 

The experimental results for the scalar slope parameter 
found by 1STRA+, KTeV and KLOE are in agreement 
with the predictions of the SM. On the other hand, the 
value found by NA48 can hardly be reconciled with the 
theoretical numbers. Furthermore, an isospin violation 
of the size as it would be suggested by the simultaneous 
validity of the ISTRA+ and NA48 results cannot be ex- 
plained within the SM (iKastner and Neufeldl |2008D . A 



clarification of the origin of these puzzling resu 



decays would be highly welcome ( Leutwyler 2009 ) 



ts in Kfj,3 



5 The ISTRA^ 



result has been rescaled by M^^ l^'^a • 



5. SU(3) breaking effects in jf "" (0) 

The value of the Kf^ form factor at zero momentum 
transfer, ]^ '^ (0) = /+(0), is the missing theoretical 
ingredient for the extraction of Vus- Within CHPT we 
can break up the form factor according to its expansion 
in quark masses: 

/+(0)-l±/p4+/p«±... . (4.51) 

Deviations from unity (the octet symmetry limit) are of 



second order in SU(3) breaking (Ademollo and Gatto 



1964 Behrends and Sirlin 19601. The first correction 



arises at Oij)^) in CHPT: a finite one-loop contribu- 



tion ( Gasser and Leutwyler 1985b Leutwyler and Roos 



1984[ ) determines /p4 = -0.0227 in terms of F.^, Mk 
and M^, with essentially no uncertainty. The jf' term 
was first estimated by Leutwyler and Roos ( 1984 ) in the 



quark model framework, leading to 

/+(0)lr = 0.961 ±0.008. 



(4.52) 



Within CHPT, /j,6 receives contributions from two- 
loop diagrams, one-loop diagrams with insertion of one 
vertex from the p^ effective Lagrangian, and tree-level 
diagrams with two insertions from £^4 or one from £p 



( [Bijnens and Talavera[ |2003[ [Post and Schilcherl |2002[ ) : 



phi xloop 



/p^ = !;.^{i^) + /,r'°°p(M) + iT^y^) ■ (4-53) 



Individual components depend on the chiral renormaliza- 
tion scale /i, their sum being scale independent. Using 



M„ 



= 0.77 GeV, one has (Bijnens and Talavera 
2003|) jl^^°°'^{Mp) + jh''^°°'^{Mp) = +0.0093 ± 0.0005. 



J pa 



■' p" 



The p^ constants appearing in /'e^^ could be determined 



phenomenologically (Bijnens and Talavera 2003), pro- 



vided the experimental errors on the slope and curvature 
of the scalar form factor reach the level AAq ^ 0.001 and 
AAq ~ 0.0001, which is unfortunately not achievable in 
the near future. Therefore, further theoretical input on 
/+(0) is needed. 



In Cirigliano et al. (2005) a (truncated) large- A^c es- 



timate of /S'^'' was performed. It was based on match- 
ing a meromorphic approximation to the (SPP) Green 
function (with poles corresponding to the lowest-lying 
scalar and pseudoscalar resonances) onto QCD by im- 
posing the correct large-momentum fall-off, both off-shell 
and on one- and two-pion mass shells. The uncertainty 
was estimated by varying the matching scale in the range 
^ e [M.fj, IGcV], leading to 



/+(0)cHPT+i/Nc = 0.984 ± 0.012. 



(4.54) 



Finally, starting with the work of Becirevic et al. ( 2005 ) 



it has been reahzed that lattice QCD is a powerful tool 
to estimate /+(0) at a level of accuracy interesting for 
phenomenological purposes. Unquenched re sults are now 
available with both Nf = 2 + 1 and Nj = 2 ( Boyle et al. 



Lt 



20081 [20101 [Brommel eTa!] [2007{ [Dawson et al.\ |2006 
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FIG. 2 Compilation of values for [^^^[/^-(O) extracted from 
all Ki3 channels. The vertical band denotes the average. 



From AntoneUi et al. (2010a I 



Lubicz et al. 2009 Tsutsui et al. 2006). The lattice re- 



sults agree quite well with the Leutwyler-Roos estimate, 
while the analytic approaches tend to be higher as a con- 
sequence of incl uding the large and p ositive (~ 0.01) 
two-loop effects (Bijnens and Talavera 2003). In com- 



ing years lattice calculations will be performed closer 
and closer to the physical light quark masses and im- 
proved chiral extrapolations can be expected. For our 
phenomenological extraction of Vus we take as reference 
value the most recent lattice result with Nf = 2 -f 1 dy- 



namical flavors (Boyle et al. 2010 1, which is also the most 



:47)(i4), 



(4.55) 



precise result for /+(0) 

/+(0)rbc/ukqcd = 0.9599(34) (: 

where the first error is statistical, the second is due to the 
uncertainties in the chiral extrapolation and the third is 
an estimate of discretization effects. Currently the domi- 
nant systematic uncertainty arises from the extrapolation 
of lattice results, obtained with unphysical quark masses. 



to physical light quark masses (Bernard and Passemar 
2010). More progress on /+(0) is expected soon from 



other lattice collaborations. 



6. Determination of Vus and CKIVl unitarity tests 

The combination |T4s|/-i-(0) can be extracted from 
both charged and neutral K decays and its value is dom- 
inated by K^ modes (see Fig. [2 ) . Using the experimental 
averages from AntoneUi et al. (2010a I, one obtains 



\Vus\f+{0) ^ 0.2163 ± 0.0005. 



(4.56) 



Taking as a reference value for /+(0) the lattice result of 
Boyle et al. (2010), \Vus\ is given by 



\Vus\{Ki3) = 0.2255 ± 0.0005exp ± 0.0012th, (4.57) 



|Vus 



0.227 


unitarity ».\ 


-^'|Vi,7Vud|(Kp2/lTii^; 


0.226 


M^ \ 


i fit 


0.225 


r-^ 


(^^^ |Vus|(K,3) 


0.224 

0.223 


|Vud|(0"-» 0")^> 


\ : 



|Vud 



FIG. 3 Graphical representation of the current status of 
IKid], I Vus I and the corresponding CKM unitarity test. The 
horizontal band represents the constraint from Kg-^ decays, 
the thin vertical band the constraint from 0''' — >■ 0"^ nuclear 
decays, the oblique band the constraint from A'^2/7r;j2, and 
the ellipse is the \a fit region. From [Cirigliano and Neufeld| 
(|20TI|. 



where we have explicitly displayed the current experi- 
mental and theoretical uncertainties [dominated by 
/_|.(0)]. A smaller value would be obtained with the 



analytical /+(0) value in Eq. (4.54). As discussed by 



AntoneUi et al. (2010a), one can perform a fit to \Vud\ 
and \Vus\ us ing a s input the values of \Vus\ fro m K j^ 
decays [Eq. (|l57f], |K./Kd| from K^/Trn [Eq. ^Jj^], 
and |Kd| = 0. 97425 ± 0.00022 from super allowed nuclear 
/3 transitions (Hardy and Towner 2009). The outcome 
is (Cirigliano and Neufeld 2011) 



\Vud\ = 0.97425 ± 0.00022, 
\Vus\ = 0.2256 ±0.0009, 



(4.58) 



with x^/ndf — 0.012 and negligible correlations between 
\Vud\ and \Vus\. Fig. ^provides a graphical representation 
of the various constraints in the \Vus\ - \Vud\ plane and 
the Icr fit region. These values [together with the neg- 
ligible contribution from \Vub\ = 0.00393(36) (AntoneUi 
et al. 2010b l] can be used to perform a very stringent 
test of CKM unitarity or, equivalently, of the universality 
of quark and lepton weak charged-current couplings. For 
the first-row unitarity sum we find 



iCKM 



= iKdl 



IK, 



iKfcl'- 1 = 0.0001(6). (4.59) 



This constraint allows one to set bounds on the effective 
scale of operators that parametrize new-physics contri- 
butions to AcKM (Cirigliano et al. 20101. The effective 



scale is constrained to be A > 11 TeV (90% C.L.), which 
puts this low-energy constraint at the same level as the 
bounds from Z-pole measurements. 



16 



7. T violation in A'^a decays 

The transverse muon polarization in i^T^s decays 

Pt = '^' ^^"^ "" f'^ (4.60) 

violates T in the absence of final-state interactions (FSI) 



(D'Ambrosio and Isidori 1998). In the case of K^ 



TT~ ii'^v, with two charged particles in the final state, the 
electro magnetic interac tion g enerates (Pt)fsi ^ 10"'^ 
(lOkun and Khriplovich 1968 ) . In K^^ decays, this effect 



does not exceed 10 '^ { Efrosinin et al. 2000 Zhitnitsky 



1980) and T-violating effects could be important. 

The S M CP-violating contribution to Pt is very small 



10 ' (Bigi and Sanda 



2000 



Cheng, 1983). Therefore, 



the measurement of the transverse polarization of muons 
in iC„, is regarded as a sensitive probe for physics beyond 



/j3 



the SM (Kohl 2010 Paton et al. 20061 



The present experimental value ( Abe et al. 2006 ) 



Pt 



-0.0017 ±0.0023stat± 0.0011 



syst 



(4.61) 



is consistent with no T violation and corresponds to the 
limit \Pt\ < 0.0050 (90%C.L.). The sensitivity of TREK 



(Kohl 2010 



Paton et al. 2006) will be able to improve 



this current upper limit by at least a factor 20. 



D. Kti-, 

The radiative Ki-^-f decays {(. = e, /i) 

K+{p) -^ ^°(p')^+(w)^£(p.)7('Z), 

K\p) -> T:-{p')t+{pt)v,{p,)-i{q) (4.62) 

allow us to perform quantitative tests of CHPT, thanks to 
theoretical developments over the past couple of decades 
as well as recent and ongoing high-statistics experimental 
studies. The decay amplitude can be written as [we focus 
for definiteness on K^^^, t he generalization is st raightfor- 
ward and can be found in Bijnens et al. (1993)]: 



Gr f 






(4.63) 



+ 



2pe • q 



'—u{P,y)l'' (^ - 75) ("^e -|>p - ^j It^viPe) 



where the first and second terms correspond to diagrams 
a) and b), respectively, in Fig. H The hadronic matrix 
elements are defined hy {J ^ V, A) 



^ liu — ^ 



dxe^'-(7r-(p')|r(\^-"(x)JJ(0))|i^°(p)), 



F„, 



{^-{p')\v:{x)\K\p)), 



(4.64) 



with the weak and electromagnetic currents defined in 
Eq. ( |4.23| ). The Ward identities g^ V"^^ = F^ and 
9^ A^v — guarantee the gauge invariance of the total 
amplitude. 




a) b) 

FIG. 4 Diagrams describing the if^^ amplitude. 



The total amplitude can be decomposed into an "inner- 
bremsstrahlung" (IB) and a "structure-dependent" (SD) 
part, both gauge invariant. IB captures the infrared sin- 
gularities according to the |Low| ( |1958J ) theorem and the 
SD part con tains terms oi 0{q) and higher. In Gasser 



et al. ( 2005 ) this decomposition was performed in such a 



way as to guarantee that the SD amplitude is regular in 
the Mandelstam plane, except for the branch points re- 
quired by unitarity. In this treatment the tensor V^i, has 
an IB component calculable in terms of Kg^^ form factors 
f± and a purely SD component, while the tensor A^i, is 
purely SD. The SD amplitudes can be parametrized in 
terms of eight structure functions Vi,Ai {i = 1,...,4). 
Defining W — p ~ p' ~ q, one has ( Gasser et al. 2005 ): 



v: 



SD 



K^ = ie^,upa{Alp"'q'' +A2q''W) 

+ ie^xpaP^q'W^ {A^W, + Aip'^) , (4.65) 



Vi (Pf^qi^ - di^i^p -q) +V2 {w^qi, - g^j^vW ■ q) 



+ V3{q-Wp'^W,-p'-qW^W,) 
+ Vi{q-Wp'^p',^p' -qW^p',). 



(4.66) 



Early theoretical calculations of Kn^^ were ba sed on 



Low's theorem and current algebra (Fearing et al. 1970). 



Modern calculations (Bijnens et al. 



2005 Holstein 1990 



Kubis et al. 2007) 



1993 



Gasser et al. 



lave been per- 
formed within CHPT, which provides a natural frame- 
work to systematically expand the hadronic amplitudes 
-P)ii A^^, and V^^,. The chiral expansion for V^^, contains 
both IB and SD terms and starts at ©(p^). A^y starts 
at ©(p^), with the leading contribution generated by the 
WZW functional ( jWess and Zumino[|i97l||Witten[ [19831 ) 
accounting for the chiral anomaly. 

The first complete analysis to 0{p'^) was performed by 
Bijnens et al. (1993) who calculated the branching ratios 



for all if^3-y modes for given cuts in the photon energy 
and in the photon-electron opening angle in the kaon rest 
frame: E*^ > ^™*, 6**^ > 6*™*. Recently, the CHPT anal- 
ysis was revisited and extended to 0{p^ 



for K2i^. ( Gasser 



et al. 



2005 



represent t 



and Kf^^ decays (|Kubis et al. 2007|), which 



e3T 



le theoretical state of the art. To 0{p ) the 



axial form factors Ai are constant, while the vector form 
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factors Vi receive contributions from both LECs (Lg j^q) 
and loops. Since all cuts in loop functions lie far outside 
the physical region, the Vi are constant to good accu- 
racy. In order to gain control on the size of higher-order 
corrections, Gasser et al. ( 20051 ) and Kubis et al. (2007) 
performed a complete p" analysis of the axial terms Ai 
(one loop) and determined the 0{p^) Li x Lj contribu- 
tions to the vector terms Vi. Despite the appearance of 
cuts in the physical region, to 0{p^) the real parts of Vi 
and Ai are well approximated by smooth functions, with 
dominant uncertainties coming from the 0{p^) LECs. 



The updated theoretical analysis (Gasser et al. 2005 



Kubis et al. 2007 ) leads to very stable predictions for the 



relative branching ratios defined by 



n / pcut 



ncut 

67 



Tk^,^{e;>e'^^\ 



> 



gcut\ 



K, 



(4.67) 
It turns out that this ratio is very insensitive to the details 
of the non-radiative K^^ form factor f+{t). Moreover, the 
SD terms contribute only a ~ 1% correction to the IB 
result for R. The final predictions are (for representative 
cuts^ for which experimental data are available): 



Rki^^ (-B™*= 30MeV, 6i™'= 20°) = 0.0096(1), (4.68) 
i?^± (^™*= lOMeV, 26°< 6l*^< 53°) = 0.00559(6). 

(4.69) 



'-K- 



37 



In order to measure the SD terms one needs to resort 
to differential decay distributions. A natural first observ- 
able to consider is th e photon spectrum dV/dE *. It was 
found (Gasser et al. 



( Alavi-Harati et al. 2001b Alexopoulos et al. 2005 Am- 



brosino et al. " p008 Lai et ar| "2005) that are not fully 
consistent (scale factor 5* = 1.9): 

i?^_L^ (£;;"*= 30 MeV, 61™*= 20°) = 0.935(15) x 10"^ 

(4.71) 



nil ^ T 



This PDG average (Nakamura et al. 2010) is about 



one signia below the theoretical prediction given in 
Eq. (4.68 1. The most recent individual measurement 



from KLOE (Ambrosino et al. 2008), 



i?^!, ^ - (0.924 ± 0.023stat ± 0.016,yst) x 10"^ (4.72) 

does not have sufficient precision to resolve the ten sion 
between the NA48 resuh R = (0.964 ±0.013) x 10"^ (|LaI 
20051) and the KTeV resu lt R = (0.916 ± O.OlTj 



et al. 



10 '' (Alexopoulos et al. 2005). By m easuring the pho- 
ton spectrum, the KLOE collaboration ([Ambrosino et al. 



2008 ) performed a fit to the SD contribution {X), finding 



(^> 



-2.3±1.3stat±1.4, 



syst 1 



(4.73) 



in agreement (within the large uncertainty) with the 
CHPT prediction {X) = -1.2 ± 0.4. 



For the charged mode, the PDG average (Nakamura 
et al. |2010 l i s dominated by three measurernents (jAkF 
menko[ |2007{ [Barmin et aL[ |1991[ [Bolotov et aH|1986| ), 
again not fully consistent (scale factor S* = 1.3): 



^Kt^S""'. 



lOMeV, 26° < 61* < 53°) = 0.00505(32). 



(4.74) 



This value is about one sigma below the theoretical pre- 



2005 Kubis et al. 2007) that in diction given in Eq. (4.69) 



both charged and neutral modes the SD terms correct 
the IB spectrum by one single function f{E*): 



dV 

dE* 



dV 



IB 



dE* 



(x) I{e; 



(4.70) 



Despite the tremendous progress in Ki^^ decays, from 
the comparison of theory and experiment we conclude 
that more accurate experimental data are desirable in 
order to perform definite tests of the CHPT predictions 
and to detect SD contributions. 



where {X) is a mode-dependent linear combination of 
phase space averages of the SD terms {Vi) and {Ai). For 
the K%^ mode, {X) = -1.2 ± 0.4 is dominated by {Vx). 
In contrast, in the Kf^ mode {X) ~ —2.2 ± 0.7 receives 
contributions of similar size from (Vi) and {Ai), thus 
making possible a detection of the effect of the chiral 
anomaly. More complicated angular distributions might 
be used to disentangle the dominant SD terms (Gasser 



et al.\ 2005 Kubis et al. 2007) 



In parallel to new theoretical developments, there has 
been considerable experimental progress in these modes 
i n the past decade. While we refer to [Nakamura et al.\ 
(2010) for a complete experimental summary we focus 
here on the Ke3-y decays. For the neutral mode, the 
world average is dominated by four recent measurements 



E. Km 

Ku is the shorthand notation for the decays 

K+{p) ^ 7r+{p,)7r-{p2)e+{pi)Mp,), (4.75) 
K+{p) -> n"{p,)TT'>{p2)e+{pi)i^i{p.), (4.76) 
K°{p) ^ n%pi)Tr-{p2)e+{pi)iyeip,) (4.77) 

and their charge-conjugate modes. In the isospin limit 
(ttIu = TTid, e — 0), the amplitude for (4.75) is given by 



T = ^k; u{p^)r(.i - l5)v{Pi) {V, A,), (4.78) 

where the last factor contains the hadronic matrix el- 
ements of the strangeness-changing vector and axial- 
vector currents, 



For more choices see Gasser et 



[ 2005 I and Kubis et al. 1 2007 1 



Vp - ^p = {'^^{Pl)'^ iP2) \sjf,U - S7^75u|if+(p)), 



(4.79) 
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with the form factor decompositions 



V^ 



H 



e^,.pa{Pi+PuYiPl +P2Y[pi -P2V, (4.80) 



^t^=-Jf- [^(Pl +P2)^ + G{pi - P2)p + Ripi + Pz.)m] ■ 



The matrix elements for the channels (4.76) and (4.77) 



can be obtained by isospin symmetry. The form factors 
F, G, R and H depend on the variables s = {pi +P2Y, 
t — {pi ~ pY and u — (p2 — pY ■ Quite often, 



s = {pi+p2f, si = {pt+Pyf 



cost 



(4.81) 



are used instead, where Ot^ is the angle of the tt"'" in the 
center-of-mass system of the two charged pions relative 
to the dipion line of flight in the rest system of the kaon 



(Cabibbo and Maksymowicz 1965 1968 



The chiral expansion of the Ku form factors was stud- 



ied in CHPT at the one- loop level (Bijnens 1990 



bach et al. 1991) and beyond (Amoros et al 



Riggen- 
2000a|b 



Bijnens et al. 1994 ). Comparison with experimental data 



yields information on several LECs. 

The analysis of K^ — )■ n^ir^e^v is a very efficient and 
clean approach to study pion-pion scattering at low ener- 
gies. In the limit of isospin symmetry, one identifies the 
TTTT phase shifts in the matrix element in a standard man- 
ner by performing a partial-wave expansion using uni- 
tarity and analyticity. Details of this procedure can be 



found in Pais and Treiman (1968a) and Berends et al. 



(1967 1968). At the end, it boils down to the following 



parametrization of the form factors. 



F 
G 
H 



Gpe'^" 



A6. 



F„e"'"cose^ 



(4.82) 



where the dots refer to neglected d-wave contributions. 
Note that the third axial form factor R gets multiplied by 
a factor ml/ Se and cannot be measured in KeA decays. 
One is therefore left with one phase difference 5 — 5s~5p 
and the four real form factors -Fs, Fp, Gp, Hp. 

Under the assumption of isospin symmetry, the form 
factors can be expanded in a series of the dimensionless 



invariants q^ = s/AM^ — 1 and Se/4M^ (Amoros and 
Bijnens 1999). Within the currently available statistics 
(Batley et al. 2010c[ ), a constant term, two slope pa- 



rameters and one curvature parameter are sufficient to 
describe the variation of the form factor F, , 



■/ 2 I i'"„4 I r/ "e 



Fs = fs + ly + fjq' + f, 



(4.83) 



4M2' 

while two terms are needed to describe Gp, 

Gp=gp + g'pq\ (4.84) 

Fp and Hp can be described by two constants. Based 
on an analysis of 1.13 x 10^ decays, the NA48/2 collab- 
oration reports the following results for the form factor 



measurements (Batley et al. 2010c I: 



fjfs = 0.152 ±0.007 ±0.005, 

Is I Is = -0.073 ±0.007 ±0.006, 

f'e/fs = 0.068 ±0.006 ±0.007, 

fplJs ^ -0.048 ± 0.003 ± 0.004, 

g^/f^ = 0.868 ±0.010 ±0.010, 

g'^/f^ ^ 0.089 ±0.017 ±0.013, 

hp/fs ^ -0.398 ±0.015 ±0.008, 



(4.85) 



with statistical (first) and systematic (second) errors. 

A preliminary branching ratio based on the NA48/ 2 
data has recently been presented ( Bloch-Devaux 2011 1, 



BR(if=' 



-e^v) 



(4.279 ±0.035) x 10-^ (4.86) 



which improv es by a factor of 3 the P PG average (4.09 ± 
0.10) X 10"^ (|Nakamura et aL[|2010|. The final analysis 



will allow to give precise absolute values for all decay 
form factors (Bloch-Devaux 2011). 



The extraction of the tttt scattering lengths from the 
phase shift measurements oi 5 = 5s — 5p requires ad- 
ditional theoretical ingredients. Forty years ago, |Roy| 
( 1971 ) established an integral equation based on analy- 



ticity, unitarity and crossing that allows to predict the 
TTTT phase values close to threshold by using experimen- 
tal data above the matching point {^/s = 0.8 GeV) and 
two subtraction constants gq and 02, the I = 0, 2 s- 
wave scattering lengths. On the other hand, using mea- 
surements and the Roy equations, one can determine 
the corresponding values of the scattering lengths. Nu- 
merical solutions of the Roy equations obtained by two 



et al. 



20100. 



groups ( Ananthanarayan et al. 2001 Descotes-Genon 



2002 ) were employed for this purpose ( Batley et al. 



Isospin-br eaking turns out to be quite substantial in 
Ku d ecays ( [Cuplov and Nehme[ |2003[ |Nehme[ |2004a|b 



2005[) . Trigge red by the precise results ofNA48/2 ( [Batley 



et al. 2008b ) , a new theoretical procedure for the treat- 



ment of isospin-breaking on K^i phase measurements 
was suggested by Colangelo et al. (2009). The measured 
phase of the 1 = s-wave is no longer (5jJ but 



ipo = 



1 



327rF2 
+ 0(p4 



(4A,+,o + s) a± + (s - Af2„)(i + _)ao 

(4.87) 



In this expression, F is the pion decay constant in the 
limit of chiral SU(2), A^+^o = Af^± - M^o, R is the 
quark-mass ratio defined in Eq. (4.11) and 



l-4A/2,/s 



(a; = ±,0). 



(4.88) 



Although the difference between the mass-symmetric 
angle (A^+^o = 0, 1/R = 0, a± = ctq) and ^0 is rather 
small (10 to 15 mrad) over the whole range accessible 
in Kei decays, there are non-negligible effects on the ex- 
traction of the scattering lengths (Batley et al. 2010c I, 
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35 NA48/2 Ke4 (2003-2004) 
Q 3 : o no isospin corrections 

0.25 - 

? 0.2 

1 

'0.15 
0.1 - 
0.05 - 




0.28 



• with isospin corrections 



0.3 



0.32 



0.34 



0.36 0.38 0.4 

M„ (GeV/c') 



-0.01 

02 

-0.02 



-0.03 - 



-0.04 - 



-0.05 - 



-0.06 - 



NA48/2 Ke4 (2003-2004) 




no isospin corrections 



0.2 



witli isospin corrections 



0.21 0.22 0.23 0.24 0.25 0.26 0.27 



FIG. 5 Measurements of the phase shift 5 without (open 
circles) and with (full circles) isospin mass effect corrections 
from NA48/2 K^i data. From |Batley et fflZ.|pOlOc| , Copyright 
CERN for the benefit of the NA48/2 Collaboration 2010. 



as can be seen in Figs. [5] and [6] Note that these ef- 
fects had not been taken into account in the analyses of 
the older Geneva-Saclay fRossclct et al. 1977) and BNL 



E865 (Pislak et al. 2001. 2003) experiments. 



FIG. 6 Fits of the NA48/2 KeA data in the (00,02) plane 
without and with isospin mass effects. The wide band (UB) 
refers to the "universal band" [see Batley et al. (2010c I and 
[Descotes-Genon et al.\ ( |2002| | for details]. The large ellipses 
are 68% C.L. contours of the two-parameter fit leading to 



( |2010c| ). 



Eq. (4.891 and the circles the one-para meter fit of Eq. (4.911 
imposing the CHPT constraint of Eq. ( |4.90[ |. The small el- 
lipse corresponds to Eq. (4.921. Adapted from Batley et al 



A two-parameter fit of the NA48/2 data leads to the 
result (iBatley et al.\ |2010c|) 



ao = 0.2220 ±0.0128 ±0.0050 ±0.0037, 

02 = -0.0432 ±0.0086 ±0.0034 ±0.0028, (4.89) 

where the errors refer to statistics, systematics and the- 
ory. Alternatively, using the additional theoretical con- 
straint IColangelo et al.\ |2000[ |2001a|b[ ) 



02 



-0.0444 + 0.236(ao - 0.220) - 0.61(ao - 0.220)^ 
-9.9(ao- 0.220)3 (4.90) 



(with a width of ±0.0008), a one-parameter fit was per- 
formed giving (Batley et al.\ 2010c I 



ao = 0.2206 ± 0.0049 ± 0.0018 ± 0.0064, 



(4.91) 



corresponding to the value 02 = —0.0442 from Eq. ( |4.90 |. 
These results should be compared with the most precise 
prediction to NNLO in CHPT ( [Colangelo eFam2001a|b| : 



ao = 0.220±0.005, aa = -0.0444±0.0010. (4.92) 

A measurement of K+ -^ -K^ii^e+v by the KEK-E740 
Collabo ration with a data sam ple of 216 events was pub- 
lished in Shimizu et al. (2004 1. However, due to large sys- 



tematic errors, this result was not included in the PDG 



fit 'QR{K+ 



7r"7rOe+z.) 



(2.2 ±0.4) X 10-5 (Nakamura 



et al. 2010|), which uses data from an old low-statistics 



experiment (Barmin et al. [1988) . Based on the analysis 
of 4.4 X 10^ V^^ events, the NA48/2 Collaboration has 
recently reported the preliminary branching ratio 



BR(is:^ 



^o^o„± 

TT TT e 



v) = (2.595 ± 0.042) x 10-^ (4.93) 



which corresponds to a factor 10 improvement compared 
to the PDG value. They could also show the consistency 



of the Fs form factors in the K^l and K'^^ modes. 



Experimental results for K]^ 



are also 



available (Batley et al. 2004a Carroll et al. 1980 MakoflE 



et al. 1993). Althou gh the present statist ics |5464 events 
in the experiment of Batley et al. (2004a)] is not compa- 
rable with the charged mode, the branching ratio and 
the form factors have been determined and a value for 
the chiral LEG L3 was obtained: 



La = (-4.1 ±0.2) X 10 



-3 



(4.94) 



The contribution of the decay chain K''^ — > [K^e^Vf.) — > 
■K'^'r:'^e~v^ is suppressed by at least four orders of magni- 
tude in the branching ratio (Chizhov 2007|. 
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F. K, 



A. A" — > TTTT 



Except for the charge-conjugate modes, there are four 
kinematically allowed Kf,^ decays: 



K+ -^ 7r+TT-7r"e+iye, 



K+ 



7r"7r"^0e+i., 



K" -^ TrVTT^e+i^e, K° -^ -k+tt-t:- e+v^, (4.95) 



related by the isospin relatfon (Blaser 1995) 



2r(i^+^7r°7r"7r°e+z/e) = (4.96) 

V{Kl ^ 7r±7r=F7r±e=Fzy^) - T{Kl -> 7r%°7r±e=Fz/e). 



At lowest order in CHPT, only the vector current con- 
tributes. Because of the small phase space, the branching 
ratios are suppressed: 0(10^^^) for the K^ dec ays and 



0(10 ^^) for the K^ decays at leading order (Blaser 
1995[). Experime ntally, there is only an upper bounc 
( [Nakamura et aL[|2010"| 



BR(X^ 



7r°7r°7r°e+i/J < 3.5 x 10" 



(4.97) 



V. NONLEPTONIC DECAYS 



In this section, we first review the status of iiT — > 27r de- 
cays. NLO calculations including strong isospin breaking 
and electromagnetic corrections allow for the extraction 
of the lowest-orde r coup lings Gs, G27 in the nonleptonic 
weak Lagrangian (2.161 from the decay rates. A careful 



treatment of isospin violation is essential for a reliable 
determination of the s-wave tttt phase shift difference 
5q{Mk) — 52{Mk)- Despite a lot of efforts, the accu- 
racy of theoretical predictions of the ratio e'/e still does 
not match the experimental precision. A NLO calcula- 
tion of iiT — >■ Stt decays significantly improves the com- 
patibility between theory and experiment even though 
some assumptions must be made about NLO LECs. Both 
isospin-violating and CP-violating effects seem to be too 
small to be detected experimentally at this time. Some- 
what unexpected, investigation of a cusp effect in the 
Mt^Ot^o distribution near threshold in iiT — >■ Stt decays has 
led to very precise values for the s-wave tttt scattering 
lengths. 



The amplitudes for the two-pion decay modes can be 
parametrized in the form 



A(/i:" 



A{K° 



A{K+ 



= A+^=Au 



1 



= Aoe 



2X0 



V2 
-.A2e'^^ 



{A 



^"tt") 



^+7r°) 



/2 

1 

72 

^00 = ^1/2 - V2 (^3/2 

Aoe}^" - V2A2e'''\ 
3 



3/2 



A~ 



5/2J 



A 



5/2J 



= ^+0 = 

z 
= ^-A+e'^t. 



A 



3/2 



-A 



5/2 



(5.1) 



This parametrization holds for the infrared-finite am- 
plitudes where the Coulomb and infrared parts are re- 
moved from A^ , j4-|_o 

amplitudes Aa/ (A/ = 



2004a I. The 



(Cirigliano et al. 
1/2,3/2,5/2) are generated by 
the A/ component of the electroweak effective Hamil- 
tonian in the limit of isospin conservation. In the SM, 
the A/ — 5/2 piece is absent in the isospin limit. More 
precisely, A5/2 = in the absence of electromagnetic in- 
teractions and therefore A2 = A2 ■ In the limit of CP 
conservation, the amplitudes Ao,A2 and A2 are real and 
positive by definition. In the isospin limit, the phases 
Xi (I — 0, 2) can then be identified with the s-wave tttt 
scattering phase shifts ^J{^/s = M^)- 

In the isospin limit (phase space is calculated with 
physical meson masses), the amplitudes Aq,A2 and the 
phase difference xo ~ X2 can be obtained directly from 
the three K — >■ tttt branching ratios [we use the recent 
compilation of Antonelli et al.\ ( 2010a )] : 



Xo 



A2 

'X2 



(2.704 ±0.001) X 10"^ GeV, 
(1.210 ±0.002) X 10"* GeV, 
(47.5 ±0.9)°. (5.2) 



In the presence of isospin violation and to lowest order 
in the chiral expansion, the amplitudes A/^i are given by 
(omitting isospin violation in the 27-plet amplitudes) 

72. 



^1/2 = ^G27Fo (M^o - ^^^0) 



V^GsFo 



(M|„ - Ml.) {l 
^e^Fj (5,^k + 2Z) 



.(2) 



3V3 



^3/2 — 



10 



G27^0 {Mlo - Mlo) 



A 



5/2 



+ Gsi^o 



0. 



(M 



_R-o 



Mio 



.(2) 



3V^ 



'i^o (5owk + 2Z) 



(5.3) 
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The various couplings are defined in Eqs. (2.161 and combination with real photon emission. Details on this 
(2.18). Following Cirigliano et al. (2004a I, we have ex- decomposition can be found in Cirigliano et al. (2004a I. 



pressed all amplitudes in terms of the neutral pion and 
kaon masses. The parameter i^o can be identified with 
the pion decay constant Ft,: at this order. The effect 
of strong isospin breaking is entirely due to t he lowest- 
order TT^-rj mixing angle e^^' given in Eq. (4.11 1. Electro- 



magnetic interactions contribute through mass splitting 
(terms proportional to Z) and i nsertions of ffewk- As a 
consequence of CPS symmetry (Bernard et al.\ 19851, 



electromagnetic corrections to the octet weak Hamilto- 
nian do not generate a A/ = 5/2 amplitude at lowest 
chiral order. 

The NLO amplitudes in the isospin limit for both 
TTTT and /C — )■ Stt were first calculated by |Kani-| 



K 



bor et al. (1991). Their main conclusion was that the 



octet amplitude Aq is strongly enhanced by final state 
interactions at NLO whereas the 27-plet amplitude A2 is 
only mildly reduced. Since the generically small isospin- 
violating effects are enhanced in subdominant amplitudes 
with A/ > 1/2 because of the A/ = 1/2 rule, a system- 
atic treatment of isospin violation is called for. 

The first complete analysis of isospin breaking in 
K —>■ TT TT amplitudes was carried out by [Cirigliano et oL 
( |2004a ). References to earlier work can be found there. 
The analysis was repeated and extended to X — > Stt am- 
plitudes by Bijnens and Borg ( 2005bJ. In the follo wing, 
we update the analysis of Cirigliano et al. ( 2004a ) with 
new experimental input (Antonelli et al. 2010a) and with 



new information on the LECs involved, both in the strong 



( Cirigliano et al. 2006 1 and in the electromagnetic sector 



( Ananthanarayan and Moussallam 2004) 



As usual at NLO in the chiral expansion, both one- 
loop diagrams with LO couplings and tree diagrams with 
a single insertion of NLO couplings must be taken into 
account. Including the leading isospin-breaking correc- 
tions (proporti onal to (?§), the amplitu des Aa/ have the 
following form ( Cirigliano et al. 2004a I : 



(27) 



Aai - G27i^.(M^„-Af^„)A^^; 
+ GsfJ {Mlo~Mlo 



1(8) 



r(2)A(. 



A/ 



e^Fl 



4(7) , y AZ) , „ Ag) 



(X) 



(5.4) 



Confronting the amplitudes in Eq. (5.4 1 with experi- 



mental rates (Antonelli et al. 2010a), one can extract the 
LO couphngs Gs,G27 and the phase difference Xo — X2- 
Using i nste ad the dimensionless couplin gs g 8,927 defined 
in Eq. (3.8), an update of the analysis of iCirigliano et al. 
(2004a) produces the results in Table IV 



Comments. 

i. The reduction of the octet coupling gs at 0{GfP'^) 
corresponds to the enhancement of the A J =1/2 
ampl itude by roughly 30 % at NLO ( Kambor et al.\ 
The (rather generous) theory errors of gs, 



1991) 



327 at NLO account for the uncertainties of the vari- 
ous LECs involved. They should not be interpreted 
to include also effects of NNLO or higher. In fact, a 
calculatio n of the leading (double) chiral logarithms 
at NNLO ( Biichler 2006 ) confirmed that tttt rescat- 



tering generates an additional enhancement of Aq 
at higher orders ( [Pallante and Pich[ |2000[ |2001| ). 
However, the values of gs, (727 displayed in Table 



IV are the appropriate values for nonleptonic weak 



amplitudes at LO or NLO, respectively. 

ii. At NLO, both gs and (727 receive small shifts from 
isospin-violating corrections. While this could be 
expected for gg,, it results from a cancellation of dif- 
ferent effects in the case of (727- On the other hand. 



as indicated in Table IV the inclusion of isospin 
breaking reduces 327 by roughly 10% at LO. 

iii. The phase difference xo ~ X2 is taken as a fit pa- 
rameter because the phases cannot be reliably cal- 
culated to NLO in CHPT. At NLO, the fitted phase 
difference is rather sensitive to isospin violation. 
The more interesting difference of s-wave tttt scat- 
tering phase shifts in the isospin limit Sq — 82 at 
-/s = Mk will be discussed below. 

iv. At NLO, electromagnetic corrections induce a 
A/ = 5/2 amplitude: 



h/2 = ReA2/RcA+-l 

= (8.44 ± 0.02cxp ± 2.5th) x 10"^ 



(5.5) 



The meaning of the amplitudes A^j can be inferred from 

the superscript X. A^/' ^A/ represent the octet and 27- 1- tttt phase shifts from K ^ tv-k decays 

plet amplitudes in the isospin limit. A^| represents the 
effect of strong isospin breaking, while the electromag- 
netic contribution is split into a part induced by pho- 
ton loops A^j and the parts induced by insertions of 
Z and 5owk vertices [A^ai ^'^cl A^^\, respectively). The 
photon loops are actually infrared divergent: A^^j are 
the infrared-finite (structure-dependent) terms whereas 
the divergent "infrared components" must be treated in 



In the isospin limit, the phase difference xo ~ X2 acces- 
sible in _ftr — )■ TTTT decays equals the s-wave tttt phase shift 
difference 5q{Mk) — 52{Mk) (Watson's theorem). It has 
been a long-standing problem to reconcile the phase shift 
difference extracted from K — >■ tttt decays with other de- 
terminations of pion-pion phase shifts. This problem be- 
came especially acute after the precise determination of 
TTTT phase shifts from combining dispersion theory with 
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TABLE IV Weak couplings gs,g27 and phase difference xo — X2 at LO and NLO, with (IV) and without (IC) isospin violation. 
The theoretical uncertainties are twice the ones assigned by Cirigliano et al. (2004a I. 



chiral order 


isospin 


38 


.927 


Xo — X2 


LO 


IC 


4.96 


0.285 


47.5° 


LO 


IV 


4.99 


0.253 


47.8° 


NLO 


IC 


3.62 ± 0.002cxp ± 0.28th 


0.286 ± O.OOO60XP ± 0.028th 


(47.5±0.9cxp)° 


NLO 


IV 


3.61 ± 0.002cxp ± 0.28th 


0.297 ± O.OOO60XP ± 0.028th 


(51.3±0.8exp)° 



CHPT (Colangelo et al. 2001a). We review here the 



present status of this problem ( Cirighano et al. 2009 1 



• The experimental situation has substantially im- 
proved in recent years for both the K^ and Ks 
lifetimes and for the branching ratios oi K ^ tttt 
decays (Antonelli et al. 2010a). Compared to the 



analysis of Cirigliano et al. (|2004a), the new exper- 



imental information reduces the phase shift differ- 
ence by more than three degrees (with higher sta- 
tistical significance), bringing it closer to the dis- 



persion theoretical value (Cirigliano et al. 2008a). 



In the original analysis of Cirigliano et al. ( 2004a ) 
the differences 7/ = xi Si{Mk) were calculated 
in CHPT although the xi and 5i separately come 
out much too small at NLO in the chiral expansion. 
The analysis can be improved by relying to a lesser 
extent on the NLO calculation of if — >■ tttt ampli- 
tudes altogether. Using less information potentially 
increases the uncertainty but this is compensated 
by a less biased comparison with the data. The re- 
sulting estimate of isospin violation is more robust 
and it leads to a further decrease of the phase shift 
difference by nearly two degrees. 



The main idea of the alternative procedure of 



Cirigliano et al. ( 2009 ) is to use only the isospin- violating 



parts of the NLO amplitudes as theory input and to de- 
termine 5q{Mk) — 52{Mk) directly from the data. In 
contrast to the chiral corrections for the full amplitudes, 
the isospin-violating (IV) corrections are much smaller 
and therefore more suitable for a perturbative estimate. 
The amplitudes are now parametrized as 



A, 



^^g»5o(MK) 



1 

71 



]4 .52(M^) 



AA 



IV 



3 



A+o = -Aae 



1(52 (Mk) 



AA%. 



(5.6) 



All isospin violation is contained in AA^, AAq^J, AAJj^q 
that can be ex tracted from the NLO amplitudes of 
Cirigliano et al. (2004a). Since isospin violation was ne- 



glected in the 27-plet amplitudes in view of the A/ = 1/2 



rule, AAjj^ {n = -\ — , 00, -f 0) scale linearly with the 
lowest-order octet coupling g^,. 

The moduli of the amplitudes in the isospin limit are 
denoted as Aq, A2. These amplitudes together with the 
phase shift difference 5[){Mk) — 52{Mk) are then deter- 
mined directly from the rates. For this purpose, the mod- 
uli are written as 



l^+-l = 



Aq + —A2 e»('52(Mic)-<5o(Mx)) 



\A 



00 1 



An 



^-A^f,^(S:,{MK)-S„{MK)) 



\A, 



+ A^o 
3 



IVg-i5o(MK) 



(5.7) 



■A^^^(S2{MK)-So{MK)) + Ayl^^ g-«o(M^) 



In order to determine Aq, A2 and 5q{M k) ~ S2{M k) from 
the three rates, also the / = phase 5q{Mk) is needed 



as input (Colangelo et al.\ 2001a) 



'5o(Mk) = (39.2 ±1.5)= 



(5.8) 



From the structure of the moduli in Eq. (5.7) it is obvious 



that the precise value of 5q{Mk) has little impact on the 
phase shift difference. 



With the (updated) values of AA^^ (Cirigliano et al. 



2004a ) a nd with gg = 3.6 (see Ta ble IV), the experimen- 
tal rates ( Antonelli et aL[[2010a l give rise to 



Ao = (2.7030 ± 0.0008) x 10"^ GeV, 
A2 = (0.1249 ± 0.0003) X 10"^ GeV, 
5q{Mk) " 52{Mk) = (52.54 ± 0.83)°, (5.9) 

where the errors are purely experimental. The octet en- 
hancement in K — > TTTT decays is characterized by the 
amplitude ratio 



£0 
A. 



= 21.63 ±0.04, 



(5.10) 



again with experimental error only. 

There are various sources of theoretical uncertainties 
associated with the phase shift difference. Whereas the 
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error of Sq in Eq. (5.8) is completely negligible, a gener- 



ous error of 20% for the overall scale gs of the isospin- 
violating amplitudes gives rise to an uncertainty ±1.1° 
for the phase difference. The major part of the error is 
due to unknown effects of 0{e^p'^) that were estimated 
by iCirigliano et al. (20091 in two different ways. Alto- 



gether, the final value for the phase shift difference in the 
isospin limit from K — >■ tttt decays is 



= (52.5±0.8exp±2.8th)°. 



(5.11) 



Before comparing this result with other determinations, 
one should recall that the isospin limit is defined in terms 
of the neutral meson masses. However, a LO estimate 



( Cirigliano et al. 2009 ) suggests that expanding instead 



around the charged pion mass modifies the phase shift 
difference only by much less than a degree. 

The most recent determinations of the phase shift dif- 
ference from TTTT scattering data are 



(47.7 ±1.5)= 



( Colangelo et al.\ 2001a I 



(50.9±1.2)° ( [Kamihski et a/.[|2008[ ), 



(47.7 ±0.4)° 



(Batley et al. 2010c), 



(5.12) 



(47.3 ±0.9)° (Garcia-Martm et al. 2011). 



These values would agree perfectly well with the phase 
difference xo ~ X2 in Eq- (5.2 1, obtained in the isospin 
limit. However, due to the large ratio Aq/A2 \d. K —^ 
27r decays, isospin-breaking corrections to the dominant 
A/ =1/2 amplitude generate sizable contributions to 
A2, modifying also the amplitude phases. The updated 
determination of the phase shift difference from K — > 
27r decays in Eq. (5.11) (Cirigliano et al. 20091 turns 



out to be in reasonable agreement with the tttt results in 
Eq. (5.12), although with a larger uncertainty. 



^7^ 



The CP-violating ratio e'/e constitutes a fundamental 
test for our understanding of flavor-changing phenomena, 
e and e' parametrize different sources of CP violation in 
Kl — >■ tttt: 



%0 



A(Kl -^ TT+TT-) 



A(Ks 
A{Kl 



TT+TT ) 



A{Ks -^ ttOttO) 



= e + e 



2e'. 



(5.13) 



The dominant effect from CP violation in K^-K^ mixing 
is contained in e, while e' accounts for direct CP violation 
in the de cay amplitudes. The present experimental world 
average (Abouzaid et al.\ 2011a Alavi-Harati et a/.]|1999 



2003a[ |Barr et al.\ |1993| [Batley et al.\ |2002| [Burkhardt 



et al. 


1988 


et al. 


2001 



Fanti et aL[ [19991 [Gibbons et a/. [ [19931 [Lai 



Re{e'/e) = \[l- 



Vo. 



V+ 



= (16.8 ±1.4) X 10"\ 



(5.14) 
demonstrates the existence of direct CP violation in K 
decays. 

When CP violation is turned on, the amplitudes Aq, 
A2, A2 acquire imaginary parts. To first order in CP 
violation, e' is given by 



'T2 



J(X2~Xo) 



ReA2 
Revlo 



ImAo 



Im^o 

ReAo ReA2 



(5.15) 



Since IniA/ is CP-odd the quantities ReA/ and xi a-re 
only needed in the CP limit (J = 0,2). e' is suppressed 
by the small ratio ReA2/ReAo ~ 1/22. The phase (j)'^ = 
Xi — Xo +7i'/2 = ( 42.5 ±0.9)° is very close to the so-called 
superweak phase ([Nakamura et al\ 2010 ) 



2{Mk,-Mk,) 



Ks 



Kl 



tan"^ ( ^'J^^^ ::''^'" I = (43.51 ± 0.05)°, 

(5.16) 
implying that cos (0^ — 4'<l) ~ 1- 

To obtain the theoretical SM prediction for e', the 
CP-conserving amplitudes ReA/ are set to their experi- 
mentally determined values. This procedure avoids the 
large uncertainties associated with the hadronic matrix 
elements of the four-quark operators in >C^ , in par- 
ticular ReAg that involves several octet operators with 
a complicated mixing under the renormalization group. 
Thus, one only needs a first-principle calculation of the 
CP-odd amplitudes Im^o and Im^2 ; the first one is com- 
pletely dominated by the strong penguin operator Qq, 
while the leading contribution to the second one comes 
from the electromagnetic penguin Qg. Fortunately, those 
are precisely the operators that are expected to be better 
approximated through our large- A'c estimate of LECs. 

Equation (5.15) involves a delicate balance between 
the two isospin contributions. A naive estimate of IvaAi 
at lowest order in the chiral expansion, i.e., using the 
tree-level formulae in Eq. (5.3 1, results in a large numer- 



ical cancellation leadi ng to unrealistically low values of 
t' /e around 7 x 10"^ ([Bosch et al.\ |2000| JBuchalla et al 



TgMljBuras et aZ.[[2001||Buras and Jaminj |2004[ |Ciuchini| 



et al. 1999 



1993 1994 



1998). The true SM prediction 



is then very sensitive to the precise values of the two 
contributing amplitudes ( [BertoHni et al. 2000 Hambye 



et al. J2000). The one-loop CHPT corrections generate 
an important enhancement (~ 35%) of the isoscalar am- 
plitude (tttt rescattering) and a reduction of A2, destroy- 
ing the accidental lowest-order cancellation and bringing 
the SM prediction of e'/e in good agreement with the ex- 



Pallante et al. 2001). 



perimental measurement (Pallante and Pich 2000, 2001 



Owing to the large ratio ReAo/ReA2, isospin violation 
plays also an important role in e'/e. Small IV corrections 
proportional to the large octet coupling G% feed into the 
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small amplitude A2 generating relatively large contribu- 
tions, which can modify the predicted value of e' in a siz- 
able way. A systematic analysis of isospin-breaking cor- 
rections in e' was undertaken in Cirigliano et al. (2003c) 



where references to earlier work can be found. To first 
order in isospin violation, one finds 



o*(X2-Xo) 



w+ 



Im^ 



(0) 



ReA 



(0) 



(I + A0 + /5/2)- 



lmA-2 



ReA 



(0) 

2 J 



where 



u}+ = RCA2 /ReAo 
Ao 



ImAo ReA[,°^ 



ImA^"^ RcAo 



(5.17) 



(5.18) 



and /5/2 is defined in Eq. (5.5). The superscript (0) on 
the amplitudes denotes the isospin limit. 

IVD.A2 is itself first order in isospin breaking. One usu- 
ally separates the electromagnetic penguin contribution 
to lxa.A2 from the isospin-breaking effects generated by 
other four-quark operators: 



ImAs = ImA^'^P 



ImA 



non— cmp 



(5.19) 



A discussion of this scheme-dependent separation in the 
framework of CHPT can be found in |Cirigliano et al\ 
(2003c). Splitting off the electromagnetic penguin con- 
tribution to ImA2 in this way, one can write e' in a more 
familiar form as 



v2 



ImA 



(0) 



ReA 



(0) 



(1 - rieff) - 



ImA^P 



Re^ 



where 



51 IV 



f^IV - Ao - /5/2, 



RcA, 



(0) 



ImA^°"-^'"P 



Re^ 



(0) 



IthA, 



(0) 



(0) 

(5.20) 

(5.21) 
(5.22) 



The quantity i^cff includes all effects to leading order in 
isospin breaking and it generalizes the more traditional 
parameter fliv- Although Oiv is in principle enhanced 
by the large ratio Re^o /'^^A\ , the actual numerical 



/5/2, which interfere destructively with fiiv to produce 
the final result (Cirigliano et aL |2003c ): 

17off = (6.0 ±7.7) X 10-2. (5.23) 

The small value obtained for iloff reinforces the dom- 
inance of the gluonic penguin operator Qg in £'■ Taking 
this into account and updating all other inputs, the SM 



prediction for e'/e turns out to be (Pallante et al. 2001 



Pichl[2004t 



Re (e'/e) 



19 ±2 11! ±6) X 10" 



(5.24) 



in excellent agreement with the experimental measure- 
ment shown in Eq. (5.14). The first error was estimated 



by varying the renormalization scale /z between Mp and 
rric. The uncertainty induced by rris, which was taken 
in the range TOs(2 GeV) = 110 ± 20MeV, is indicated by 
the second error. 

The most critical step is the matching between the 
short- and long-distance descriptions, which was done at 
leading order in 1/Nc- Since all next-to- leading ultravi- 
olet (OPE) and infrared (CHPT) logarithms have been 
taken into account, our educated guess for the theoreti- 
cal uncertainty associated with subleading contributions 
is ^ 30% (third error). 

The control of non-logarithmic corrections at NLO 
in 1/Nc remains a challenge for future investigations. 
Several dispersive analyses (Bijnens et al. 



1999b 2001 



2001 



Cirigliano et al. 2001[ 2003aj .Cirigliano and Go 



2004 



owich 



2000 2002 Donoghue and Golowich 2000| Knecht et al. 



Narison 2001 ) and lattice calculations 



dBlum et al. |2003| [Boucaud et al.\ |2005[ [Noaki et al. 
2003 ) of {Qs) = {2'it\Qs\K} already exist (most of them 
in the chiral limit). Taking the chiral corrections into ac- 
count, those results are compatible with the value used in 
Eq. ( 5.24[ ). Unfortunately, the hadronic matrix element 
(Qg) = {"^T^lQelK) is more difficult to compute. Two an- 
alytical estimates in the chiral limit, using the so-called 



minimal hadronic approximation (Hambye et al. 



1999 



and the X-boson approach ^Bijnens and Prades 
2000a), find large 1/Nc corrections to {Qe). It would be 



2003) 



interesting to understand the physics behind those con- 
tributions and to study whether corrections of similar 
size are present for physical values of the quark masses. 
Lattice calculations of (Qe) are still not very reliable 
and give contradictory results, often with the wrong sign 



analysis shows all three terms in Eq. (5.21 ) to be relevant ( Bhattacharya et al. 20051 Blum et al. 2003 Noaki et al. 



when both strong and electromagnetic isospin violation 
are included (Cirighano et aL[|2003c). 



2003 Pekurovsky and Kilcup 2001). The usual proce- 



The numerical analysis of [Cirigliano et al. ( 2003c ) 



found large cancellations among the different contribu- 
tions to r^cff. A well-known example are the contribu- 
tions of strong isospin violation via n'^-r] mixing where 
the sum of t] and 77' exchange generates an fijv of the 
order of 25%. However, already at the level of n^-r] 



mixing alone, a complete NLO calculation (Ecker et al. 



2000b) produces a destructive interference in flw, with 
fiiv = (15.9±4.5)x 10"^. Inclusion of electromagnetic ef- 
fects slightly increases fiiv and generates sizable Aq and 



dure in lattice simulations has been to calculate K ~> tt 
and K — > vacuum matrix elements and use lowest-order 
CHPT to recover the physical iiT — >■ 27r amplitudes. 
In addition to all usual lattice artifacts (finite volume 
and lattice spacing, quenched approximation, unphysi- 
cal masses, etc.), this misses completely the crucial role 
of final-state interactions in the isoscalar octet ampli- 
tude. Thus, a truly major effort is needed to compute 
the i^ — >■ 27r matrix elements directly. Some progress 
in this direction has been achieved, relating the phys- 
ical K -^ 2Tr amplitudes to the corresponding matrix 
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elements in a finite Euclidean volume, which are bet- B. Tf — > Stt 



ter suited for lattice simulations ( Lellouch and Liischer 



2001 Lin et al. 2001 1. So far, this technique has been im 
plemented in simulations of the A 7 = 3/2 K'^ — >■ tt+tt' 



amplitude with promising results (Boucaud et al 
Goode and Lightman[ |2010| |Kim and Sachrajda| 



2005 



2010 



Liu 2010[). A first estimate of the Aq amplitude at un- 



physical kinematics has been recently reported (Blum 



et al. 2011) 



More work is needed to reduce the present uncertainty 
quoted in Eq. (5.24). This is a difficult task, but progress 



in this direction may be expected in the next few years. 



3. CP violation in K^-K'^ mixing 

Since Re(e'/e) <C 1, the ratios 7?^_ and 7?„„ pro vide a 
direct measurement of |e| (Nakamura et al\ 2010): 



There are five CP-conserving decays to three pions (the 
K~ decays are not listed separately): 

A^oo = A{KL{k) ^ ^°(pi)7r°(p2)^°(P3)), 
Al_^ = A{KL{k)^TT+{p,)7r-{p2)7r°{p,)), 

Al_o = A{Ks{k)^7r+ipi)n-{p2)n'{p3)), 
Aoo+ = AiK+ik) ^ n°ip,)7:"{p2)TT+ip,)), 
A++_ = AiK+ik)^7T+ip,)7:+ip2)n-ips)). (5.29) 

For the kinematics one uses the variables 

si = (fc - Pif , S2 = (fc - P2f , S3 = {k- psf (5.30) 



with 



So 



;{si +S2 + S3) 



1 



{^K-\ + hJ) = (2.228±0.011)xl0- 



1 



(5.25) 
in perfect agreement with the semileptonic asymmetry 

r{KL -^ n-£+i^) - r{KL -^ n+e-D) _ 2Re(e) 

r{KL -> n-e+i^) + r{KL -^ n+e-D) ^ TTW 

^ (3.32 ±0.06) X 10"^ (5.26) 
The theoretical prediction can be written in the form 



(Ml + M^i + M^2 + M^s) , (5.31) 



where the masses are those of the particles appearing in 
the decay under consideration. 

In terms of the Dalitz plot variables 



52 -Si 



y 



S3 - Sq 

7r+ 



(5.32) 



( Buchalla et al. 1996[ |Buras and Guadagnoli 2008 ) the a mplitudes are often expanded as ( Devlin and Dickey 



|e| = C,k,BKA^X^x{A^X^{l^p)rjuSo{xt) (5.27) 
+ VctSoixc, xt) - r]ccSo{xc)} , 



1979) 



where Xj = mf/ M^, So{xi, Xj) and So{xi) are the Inami 
and Lim ( 1981 ) box functions. 



a 



GIfIMkoM^ 



3V2n^iMK^-MKs 



= 3.7 X lO'^, 



(5.28) 



Vu 



0.5765 ± 0.0065, rj^ 



= 1.43 ±0.23 and 7?™^° 



0.496 ± 0.047 are short-distance QCD corrections (Brod 



and Gorbahn|[20T0{[Bm-as et aL[[T990||Herrlich and Nier 



ste 1996) and k^ = 0.94 ± 0.02 accounts for small long- 



A^oo = 3(ai + as) + 3(Ci - 2(3) y 



^+-0 = {ai + as) - Wi + I33)y 

+(Ci - 2C3) (y' + Ix' 
+(Ci - 26) (y' - Ix' 

o 2 ^ 4 , 

^+--0 = T^^^lsx- -£,3xy, 



distance contributions (Buras et al. 2010) 



i004 



The renormalization-group-invariant parameter B 
measures the hadronic matrix element (-^*'|QAS=2|-f^*') 
in units of its vacuum satur ation approximatioii. In 



(5.33) 
f -"1 + 2^3 ) + ( ^1 ~ 2^^ ~ ^^^ ) ^ 



the large-A^c limit, B = 3/4 ( Buras and Gerard 



1986 



Gaiser et al. 1981). The most precise lattice determi- 
nations, obtai ned with 2 4 - 1 act ive flavors, quote B = 
0.724 ±0.030 dAubin et ai [ioiol and B = 0.749 ± 0.027 



<Ci + C3){y' + ^^' 



-{^i+C3 + Q[y'-^x^ 



(Aoki et al. 



2oTT| : 



A, 



(-2ai +03 



) + (-I3i + 1^3 - V373 



Eq. (5.27) provides a parabohc constraint in the plane 
p-fj that is included in the SM unitarity triangle fits 



( Bona et al. 2006 Charles et al. 2005 ) . The recent pre- 



cise lattice value for B introduces some tension in the fit. 
Determining 77 and p from other observables, it i mplie s 
|e| = (1.90 ± 0.26) X lO'^ ( |Brod and Gorbahn[ |2010| , 
slightly smaller than (5.25). 



-(2Ci+2C3)h/ + -a;2 

+iii+C3-Q(y'-l 



To NLO in the chiral expansion, the K ^ 3tt ampli- 



tudes were first calculated by Kambor et al. ( 1991 ) in the 
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isospin limit. Their analysis was repeated and updated 
by Bijnens et al. (2003), confirming the main observa- 



tion of |Kambor et a/.| ( [19921 11 99 1[ ) that inclusion of NLO 



corrections significantly improves the LO current alge- 
bra amplitudes. One main reason for the much better 
agreement with experimental data is that the quadratic 
slope parameters Cii---i?3 vanish at LO. On the other 
hand, NLO CHPT only provides the leading contribu- 
tions for these quadratic slope parameters, which more- 
over depend on a number of LECs. Some assumptions 
about the combinations of NLO LECs occurring in the 
amplitudes must be made for a comparison with experi- 
ment. 

Experiments on X — > Stt decays provide information 
on the rates and on the Dalitz plot distributions. The 
latter are conventionally expanded up to second order in 
x^y (assuming again CP conservation): 



A{si, 32-, Sz) 



A(so,so,so) 



= 1 + gy + hy + kx 



(5.34) 



[see Eq. (5.33l]. The experi- 



except for Ks — >■ 7r"^7r 
mental data available at the end of 2004 were confronted 
with NLO CHPT by [Bijnens and Borg| ( [20041 [2005a|b[ ), 
including also isospin- violating and radiative corrections. 
Assuming that the additional electromagnetic LECs at 
NLO in the Lagrangian ( 2.18 ) all vanish at a certain scale, 
their conclusion was that isospin breaking in X — >■ Stt is 
in general small. That analysis seems worth repeating, 
not only because of better knowledge of LECs reviewed in 
Sec. [III[ but also because of the more precise experimen- 
tal information on the Dalitz plot slopes in Eq. ( 5.34 |, 
in particular for the K^ modes from NA48/2 [see Batley 
et al. (2010a) and references therein]. 



1. CP violation in iC — > Stt decays 

The decay Ks —> 37r° violates CP. In analogy to K^ 
2iT in Eq. (5.13), one defines the amplitude ratio 



rjooo 



^000 



4^ 

^000 



= e + e' 



000- 



(5.35) 



-y=o 



The parameter Eqqq is a measure of direct CP violation. 
To lowest order in CHPT, there is a simple relation be- 



Wolfenstein, 1980), which implies 



tween Eqqq and e' (D'Ambrosio and Isidori 1998 Li and 



•=000 



« k 



Therefore, it will be very difficult to detect direct CP 
violation in this decay. In fact, ^5 —> 37r° has not been 
observed at all so far. The accurate theoretical prediction 



BR{Ks -^ 37r°) = 1.9 x 10" 



(5.37) 



should be compared with the best upper bound ( |Am- 
brosino et al. 2005[ ) 



BR{Ks -^ 37r°) < 1.2 x 10"^ (90% C.L.), (5.38) 



TABLE V Dalitz slope asymmetries for K — >■ Stt decays: 
comparison between theory ( [Gamiz et al.\ [2003[ ) and experi- 
ment ([Batley et aL\\2007\. 



A? 



(-2.4 ±1.2) X 10"^ 
(-1.5 ±2.2) X 10"'' 



4^ 



(1.1 ±0.7) X lO"'^ 
(1.8 ±1.8) X 10"* 



Ref. 



Gamiz et al. 



Batley et al. 



(20031 



(20071 



corresponding to 



|r?ooo| < 0.018. 



(5.39) 



In the search for direct CP violation in K decays, the 
three-pion decays of charged kaons have played a promi- 
nent role. In addition to rate asymmetries, both the- 
ory and experiment have paid special attention to asym- 
metries in the linear Dalitz plot parameter g defined in 
Eq. KM. 



After a number of conflicting estimates [reviewed by 
D'Ambrosio and Isidori (1998)], the theoretical state-of- 
the-art Dalitz plot asymmetries c ome fro m a NLO CHPT 
calculation of|Gamiz et al. ( 2003 ) [see also Prades ( 2008 1 ; 



Prades et al.\ (2007)]^ In fact, this calculation also in- 



volves an estimate of the dominant contributions to the 
absorptive parts of the relevant amplitudes of 0{GsP^). 
The linear Dalitz plot slopes g^ ^ refer to the decays of 
K^ into three charged pious (C), and one charged and 
two neutral pious (iV), respectively. The CP-violating 
quantities of interest are the slope asymmetries 



aC,n _ 



9 C.N 



9c,N 



9cM 



' 9c,N 



(5.40) 



These asymmetries depend at LO on the imaginary parts 
of the LECs Gh a nd gowk [Eqs. ^JE^ and (pTSl)] only 



Gamiz et al. 



( 2003 1 found that A'^ is relatively insensitive 
to NLO LECs whereas A^ is less stable. The theoretical 
results are compared with the most precise experimental 
data in Table [Vj 

To a good approximation, e'/e also depends on the 
imaginary parts of the LO LECs Gg and gcwk only (see 
Sec. V.A.2). Therefore, one can establish bounds for 
A^ with in the SM, using the experimental value for e'/e 
([Prades[[2008|. 



(5.36) 2. TTTT scattering iengtiis from K ^ 3tv decays near threshold 



Because of the 7r+-7r° mass difference there is a cusp 
in the MttO^o distribution at MttO^o = 2Mt^+ in K ^ Sit 
decays with two tt" in the final state. It w as first seen in 



7r°7r° (Batley e^gZ, 



2006 



2009a ), more recently 



also in Kl — > Stt^ ( Abouzaid et aLl "2008a|. It is due to 
the charge exchange scattering of p i ous in the final state 
([Budini and Fonda[ [19611 [Cabibbo[ [2004| 



K^ 



~{Tr~^TT 



TT TT TT . 



(5.41) 
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FIG. 7 Interference between tree and one-loop amplitudes 
generates a cusp at M^o^o = 4M^4. in Tf — >■ Stt decays. 
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The interference between tree and one-loop amplitudes in 
Fig. [T] generates the cusp via t he sq uare-root singularity 
at A'i^o^o = 4M^_|_. From Eq. (5.41) the effect is seen to 
be mainly sensitive to the combination of tttt scattering 
lengths 



ao - a2 ~ A[tx'^ 



TT TT ) 



thresh- 



(5.42) 



Various approaches have been pursued to extract the tttt 
scattering lengths from K -^ in near threshold. 



Following the original approach of Cabibbo ( 2004 ) , 



based on unitarity and analyticity, a systematic ex- 
pansion of the singular terms of the M^OttO distri- 
bution in powers of the scattering lengths was per- 
formed by Cabibbo and Isidori ( 2005 1 . 



ii. In a related method, unitarity and analyticity were 
combined with CHPT (IGamiz ei aLl 12007). 



iii. A two-loop dispersive representation of i^ — >■ Stt 
amplitudes in the presence of isospin breaking is 
under construction ( Kampf et al. . 2009 ) . 



In the most advanced approach, based on a nonrel- 
ativistic effective field theory (NRQFT), the K ^ 
Stt amplitudes are expanded in powers of the scat- 
tering lengths and of the pion momenta in the K 
rest frame ( Bissegger et al. . 2008 Colangelo et al. 



2006 Gasser et al.. 2011). 



In the NRQFT approach of the Bern-Bonn group, the 
power counting involves a small parameter e that charac- 
terizes the size of pion three-momenta: \p\/AI„ = 0{e). 
In contrast to standard CHPT, the scattering lengths are 
not predicted but extracted from the data. The decay 
amplitudes are then given by a two-fold expansion in e 
and in the scattering lengths (denoted below generically 
as a). At any given order in a and £, only a finite number 
of graphs contribute because each loop is suppressed by 
one order in s. 

NRQFT is manifestly Lorentz invariant and therefore 
frame independent. Analyticity and unitarity are guar- 
anteed as in standard CHPT. However, unlike in stan- 
dard CHPT, the amplitudes are valid to all orders in 
the quark mass es. They have been calculated up to 
0{e* ' " ~ 



2006) 



e ) (Bissegger et al. 



2008 



Colangelo et al. 



A major advantage of the NRQFT approach is that 
photons can be incorporated in a straightforward man- 
ner, allowing for the systematic inclusion of electromag- 
netic corrections (Bissegger et al. 2009). In the Coulomb 
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FIG. 8 NA48/2 Ke4 and cusp results from two-parameter fits 
in the (ao, 02) plane. The smallest contour corresponds to the 
combination of NA48/2 results. The cross-hatched ellipse is 
the CHPT prediction ([XMl of Colangelo et aZ. (2001a b). The 



dash-dotted lines correspond to the recent result from DIRAC 
(Adeva et al. 20111. We thank Brigitte Bloch-Deva ux for 
updating the original figure from Batley et al. (2010c I. 



gauge, only transverse photons appear as internal lines 
in diagrams whereas the Coulomb photons give rise to a 
non-local vertex via the equations of motion. Of course, 
also real photon emission must be included, in particular 
to cancel infrared divergences. However, the nonrelativis- 
tic power counting then shows that finite bremsstrahlung 
effects are small near threshold. 

The production of pionium, bound states of charged 
pions, upsets the nonrelativistic counting. A certain re- 
gion around the cusp is therefore excluded from the data 
analysis such that one- photon exchange of 0{ae^) is suf- 
ficient. 

The decay spectra were calculated to 0(e^£^) for all 
K —>■ Stt channels and, in addition, to ©(e^ae^) for the 
channels of main interes t, K^ -^ 7r+7r'^7r° and Kl —> Stt" 
(Bissegger et al. 2009). The radiatively corrected am- 



plitudes were used i n the analysis of NA48/2 data by 
Batley et al. (2009a) to extract the scattering lengths. 
Fig. |8 sho ws a co mparison of the cusp an alysis of |Bat 
ley et al. (2009a) and the Ke4 studies of Batley et al 



(2010c). The smallest contour corresponds to the com- 
bination of the two methods. The result of the DIRAC 



experiment (Adeva et al. 2011) and the CHPT predic 



tion of Colangelo et al. ( 2001a|b ) are shown as well. The 



corresponding numerical results are displayed in Table 
|VI[ The agreement between theory and experiment is im- 
pressive. 
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TABLE VI Experimental and theoretical results for s-wave tttt scattering lengths. First line: cusp analysis of NA48/2 ( Batley 

The external error is due to 



et al. 2009a I using the Bern-Bonn framework (Bissegger et al. 



the uncertainty in the ratio of the amplitudes for K^ 



2008| [Colangelo et al. 

- and V"^ ->• 7r"^7r"7r". 



20061 



(Batley et al. 2010c I. Third line: analysis of Colangelo et al. (2001a I on the basis of Roy equations and CHPT 



Second line: Kei analysis of NA48/2 



Batley et al. 



Batley et al. 



( 2009a 



(2010c I 



Colangelo et al. (2001a I 



ao - a2 



0.2571(48)etat(25)syst(14)e 



0.264 ± 0.004 



ao 



0.2220(128)stat(50)syst(37)th 
0.220 ±0.005 



132 



-0.0241(129)stat(94)syst(18)oxt 

-0.0432(86)stat(34),y,t(28)th 

-0.0444 ± 0.0010 
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FIG. 9 Z-penguin and box contributions to K 



VI. RARE AND RADIATIVE DECAYS 

Kaon decays mediated by FCNC fall in the set of rare 
and radiative decays. These modes are suppressed in 
the SM and their main interest, other than their own 
understanding, relies on the possible observation of New 
Physics effects. Most of these processes are dominated by 
long-distance contributions, such as K ^ 77, K — > 77*, 
K — >■ TT"/* and others. However, there are also pro- 
cesses governed by short-distance amplitudes, such as 

K -^ 'KUV. 

Long-distance dominated decays have been studied in 
the chiral framework. The leading contributions, mostly 
O(p^), have been evaluated within CHPT. In many of 
the processes estimates of the dominant NLO corrections 
have also been carried out. 



A. K 



■KVV, TTTTUV 



The rare decays Kl 



and K — > tt 1/1/ can be 



predicted with a precision surpassing any other FCNC 
process involving quarks, thus making them a clean 
place to look for non-standard signals. The analysis of 
these processes makes use of the full arsenal of effective 
field theory, from short-distance effective Hamiltonians 
to CHPT. We outline below the main steps of the anal- 



ysis and then quote the final results for the branching 
ratios, with a brief discussion of the uncertainties. 

Let us begin with a qualitative discussion that high- 
lights the main features of these modes. In the SM, 
the K — >■ TTvV decays proceed through FCNC amplitudes 
generated at the quark level by Z-penguin and box dia- 
grams (see Fig. [9]). Separating the contributions accord- 
ing to the intermediate up-type quark running inside the 
loops, the quark-level amplitude has the structure 



A„, (6.1) 



The equa- 
tion above reflects the quadratic ( "hard" ) nature of the 
GIM mechanism. This in turn implies that the top- 
quark contribution, which carries a large CP-violating 
phase, accounts for ~ 68% of A{s -^ dvv) while the 
charm- and up-quark contributions amount to ~ 29% 
and ~ 3%, respectively. So we see that the power-like 
GIM mechanism implies a large suppression of the u- 
quark (long-distance) contribution and complete domi- 
nance of the CP-violating contribution by the i-quark 
term. These properties are often summarized by stating 
that the K -^ 'kvV decays are "short-distance" domi- 
nated. The short-distance domination also implies that 
to a good approximation the interaction is described at 
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low energy by one single local operator, Q13 in Eq. (2.9 1, 



whose hadronic matrix element can be related to the form 
factors appearing in K13 decays. All these features make 
it possible to predict K — ^ ttpT' rates very accurately. 

The above qualitative discussion can be put on 
solid footing by systematically employing effective-theory 
techniques. As usual in the application to weak decays, 
three steps are required, (i) Determine the effective La- 
grangian CcS at the weak scale /i ~ Mw by integrat- 
ing out the heavy gauge bosons and the top quark. In 
this case, the diagrams of Fig. [9] generate the operator 
Q13. The Wilson coefficient C13 is determined entirely 
in terms of At and Xt, and it now includes NLO QCD ef- 



fects (Buchalla and Buras 1999; Misiak and Urban 



and two- loop electroweak corrections (Brod et al 



[1999) 



poll 
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Buchalla and Buras 1998 ). (ii) Evolve £off down to a low 



hadronic scale using renormalization group techniques. 
In this step charm-loop contributions (the c quark is still 
an active degree of freedom) generate corrections to C13 
proportional to XcXc In Xc- These are known up to and in- 



cluding NNLO QCD effects (Buchalla and Buras 1994 



Buras et al. 2005 



tions (Brod and Gorbahn 2008) 
old jj. 



2006 1 and NLO electroweak correc- 
At the charm thresh- 
At the 



rric one integrates out the charm quark 
current level of precision, one needs to keep not only the 
leading dimension-six operator, but also dimension-eight 
operators (Falk et al.. 2001), whose relative importance 



is suppressed by M'^/m^ ~ 15%. (iii) Calculate hadronic 
matrix elements of the operators appearing in £cff at a 
low scale /i ~ 1 GeV. Ccs contains the purely semilep- 
tonic operator Q13, the dimension-eight operators gen- 
erated at the charm threshold, and the \AS\ — 1 four- 
quark operators, which generate long-distance contribu- 
tions to K — )> TTiyi? (they correspond to the u-quark con- 
tributions to Z-penguin and box diagrams). The matrix 
element of Q13 can be related to if ^3 form factors (Mar- 
1996|). This analysis was recently up 



ciano and Parsa 



dated by Mescia and Smith (2007) to include first-order 



isospin-breaking effects to NLO in CHPT as well as long- 
distance radiative corrections for the K^ decay. Con- 
cerning the long-distance contributions and dimension- 
eight operators, the relevant matrix elements were calcu- 



lated in CHPT ( [Isidori et al] [2005| [Luand Wise| p94 ) . 



Owing to the CP properties of the K2 state, the op- 
erator (5i3 only contributes to Kj^ — >■ Tr^iyv through a 
violation of the CP symmetry. The transition ampli- 
tude is completely dominated by direct CP violation, the 
contribution from K'^-K^ mixing being only of the or- 
der of 1%. There exist tiny CP-conserving contributions 
through higher-order short-distance operators and long- 
distance corrections, but their effect on the branching 
ratio is negligible ("Buchalla and Isidori 1998). 



Putting all the ingredients together, the predicted SM 
rates for K — > tti/P decays can be written as 



BR{Kl -^ TT^VV) = KL 



ImXt 
1^ 



X] (l-S,), (6.2) 



BR(A" 



(1 + Aem) 



[ A5 



X 



ReXt 



X 



RcXc 



[Pc + SPc. 



(6.3) 



where Se — V^\e\ [1 -I- Pc/jA'^ X) — p\/rj contains the smal l 
K^-K'^ mixing contribution (Buchalla and Buras 1996), 
with A = T4s, A, p and ij the Wolfenstein CKM parame- 
ters. The overall factors kl,+ encode the hadronic matrix 



element related to Kis data (Mescia and Smith 2007): 



KL 



(2.231 ±0.013) X 10" 
(5.173 ±0.025) X 10"" 



1" (A/0.225)*, 
(A/0.225)^ 



(6.4) 



They are now dominated by experimental uncertainties 
in Kg^. Defining the infrared-safe photon- inclusive rate 



by the cut E^™^ < 20 MeV, the electromagnetic cor- 
rection takes the value Aem = —0.003 (Mescia and 



X 



Smith 2007 ) . Th e top-quark contribution is given by 
1.469 ± 0.017 ( |Brod et al.\ [20TT] [Buras et al. 



2005 



2006), with an uncertainty due to the parametric error 



in rrit, matching scale /it, and higher-order electroweak 
effects. The dimension-six charm contribution has the 
value Pc = 0.38 ± 0.04, with error dominated by the 
parametric uncertainty in rric. Finally, the long-distance 
{u quark) and dimension-eight charm contributions can 
be lumped in SP^^u = 0.04 ± 0.02 ( [Isidori e^ a/!] [2005 ). 
Taking the CKM matrix elements from global fits, one 
arrives at (j Brod et al] [2011[ [Buras et al] [2005[ [2006[ ) : 



BR(A:l -^ TT^vv) = (2.4±0.4) X 10"", (6.5) 
BRiK+^-K+vv) = (0.78 ± 0.08) X 10"^". (6.6) 

In both cases the uncertainty is largely parametrical 
(~ 80% for the Kl mode and ~ 70% for the K+ mode), 
due to CKM input, TOc, rrit and as{Mz). As the deter- 
mination of CKM parameters improves in the next few 
years, we can expect to reach accuracies better than 10%. 
In summary, from a theoretical perspective K — >■ irvi? 
decays offer the cleanest window on non-standard con- 
tributions to the s ^f d transitions (complementary to 
B meson studies). Moreover, within the general context 
of flavor physics, K -^ iivD decays are perhaps the most 
promising place to look for non-standard signals, due to 
their SM-specific suppression induced by CKM factors 
[A{s — > duv) (X V^g]. Even if one considers SM exten- 
sions with the same CKM suppression factor, K — >■ tzvu 
decays still provide an excellent probe, because the SM 
contributions can be predicted very accurately. 

On the experimental side, the charged kaon mode was 
observed (Artamonov et al. :2008), while only an upper 



bound on th e neutral mode has been achieved ( Ahn et al. 
20081 [20T0| : 

(6.7) 



1.15N 



X 10" 



BR{K+ ^TT+vD) ^ {l.nt\ 
BR{Kl^tt^vv) < 2.6 xlO"* (90%C.L. 



New 



( Spadaro 



experiments are under development at CERN 
20lT| and J-PARC (Watanabel [20To| for 
charged and neutral modes, respectively. These exper- 
iments aim to reach 0(100) events (assuming SM rates), 
thus beginning to seriously probe the new-physics po- 
tential of these rare K decays. Increased sensitivities 
could be obtained through the recent P996 proposal for 
a Ar+ — > TT^vD experiment at Fermilab and the higher 



kaon fluxes available at Project-X (Tschirhart 2011). 



Finally, let us mention that the decays AT^.g — )■ mriyp 
and K^ —> tt^tt'^vv share the same feature of short- 
distance domination as the corresponding single-pion 
modes, and could therefore provide another probe of the 
underlying s — > dvv transition within and beyond the 
SM. The calculation of the decay amplitudes requires 
taking the matrix element of the current s7/i(l — 75)^ 
between the kaon and two-pion states, which can be ex- 
tracted from the measured A'^4 decays using isospin sym- 
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cia 1996 



metry ( Chiang and Gilman 2000 Littenberg and Valen- 



19941. T 



or directly calculated in CHPT (Geng et al. 



leoretical predictions can be summarized as fol- 

(6.9) 



lows (Littenberg and Valencia, 19961: 



-13 



1.8 X [(1.37-^)^ + 0.1777^] X 10 
BR{Kl -^ n^ir^i^D) ~ (1.37 - pf x 10"", 
BR(if=^ -^ TT^TT^vv) ~ 7 X [(1.37 - pf + 77^] X lO^^^ 



Experimental searches are still far above the expected SM 
rates and have reached the following 90% C.L. limits: 



BR{K+ ^TT+TT^vv) < 4.3x10"^ ( |Adler et a/.[[200T| ), 
BR{Kl -^ n°n°w) < 8.1 x 10"^ ( [Ogata et al.\^m\ 



(6.10) 



B. K 



. jMjM 



The amplitude for a transition of the type 

K(p)^7*('Zi)7*(92) 



(6.11) 



is determined by a tensor amplitude M^'^{qi,q2) that has 
the most general form compatible with gauge invariance: 



Mf"" ^ 



+ 



Q^Qi 



4 



qi ■q2 



qU2 , 
q\ 



91-92 p ^ 

, 92^92 
92 



M|:«(9?,9l) 



91-92 i, y 



X 6(9?, 9l) 
+ ie'^''P''q^pq2.c{qlql). 



(6.12) 



Thus, for instance, the amplitude into two real photons is 
given by ^ = M''''(gi, 52) ep(9i)ei'(92)- Bose symmetry 
implies that the invariant amplitudes a((7j,q|), 6(gJ,g|) 
and c{q\, g^) are symmetric functions of their arguments. 
If CP is conserved the amplitudes a and h contribute 
to Ki{Ks) — >■ 7*7* while the amplitude c determines 
K^{Kl) -^ 1*1* ■ When one of the photons is on-shell 
{q\ = for instance), AI^^'^ is described by two invariant 
amplitudes. 



M^"^ 



(92V -91 -929'^") 6(0, 92^) 
^^^''^9ip92.c(0,g2), 



-le"^ 



(6.13) 



which also remains valid for both photons on-shell. 

These processes are dominated by long-distance dy- 
namics, in particular those with both photons on-shell 
or when off-shell photons produce lepton pairs. Short- 
distance amplitudes give only tiny contributions. 

Since the photon does not couple directly to neutral 
particles, the only possible local contributions should 
come from field strength tensors, which are absent at 
0{p^). Moreover, the C(p'*) Lagrangians Cf^^Ji^ and 



-Gsp" 



^Q ~i in Eq. (|2.16p only contain two operators with F^" 





FIG. 10 Lowest-order [©(p*)] contributions to Ks 



77. 



terms, but they are coupled to at least two charged par- 
ticles. Therefore, the leading K — >■ 7*7* contribution is 
generated by a non-local 0{p'^) loop amplitude that is 
necessarily finite. In the following and unless explicitly 
stated we assume CP invariance. 



1. Ks 



77 



The photons produced in this decay have parallel po- 
larizations {Ff^yF^^") and, up to one loop, there is no 
short-distance contribution due to Furry's theorem (Gail- 
lard and Lee |1974b ). This decay arises, at leading order, 
from a finite two-pion one -loop amplitude with one ver- 



tex from C^^=2^ or C^^^^l (iD'Ambrosio and Espriu 1986 
Goity 1987), as shown in Fig. TH Its contribution is 



6^(0, 0) = ^aFo(^G8 + ^627) (1-^2)^(^1), 

"(6.14) 
where rp = Mp/Mx and F{x) is given in Eq. ( Al ). The 
charged-kaon loop contribution vanishes because of the 
factor (1 - r\) = 0. 

The decay width is given by 



F {Ks ^ 77) = 



Ml 

64 TT 



6(0,0)1 



(6.15) 



Taki ng th e LO values of the couplings Gs and G27, 
Eq. ( |6.14[ ) results in BR{Ks -J> 77) = 2.0 x 10"^, which 
compares rather well with the present experimental world 



average (Nakamura et al. 2010) 



BR{Ks -^ 77) = (2.63 ± 0.17) x 10" 



(6.16) 



Although the full 0{p^) amplitude has not been cal- 
culated, the dominant effects were implemented through 
(i) the inclu sion of unitarity corrections f rom Ks — >■ tttt -^ 
TT~^Tr~ -^11 (Kambor and Holstein 1994), and (ii) a local 
contribution introducing an unknown coupling constant 
( [Buchalla et aL[|2003[ ), 



&(4+6)(o,o) 



2aFQ 
-Ml 



B{M, 



4aGg 

ttF. 



Ml {I 



;jai, 
(6.17) 
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FIG. 11 Dominant contribution to Ki 



77- 



where B{s) can be found inlKambor and Holstein ( 1994). 
Including the 27-p let contribution in B{s), the exper- 
imental rate (6.16) implies ai = (—1.2 ± 1.3) x 10^'^, 
showing that unitarity corrections are enough to repro- 
duce the measured branching fraction. 



2. Kr 



■11 



This decay produces photons with perpendicular po- 
larizations [Sfiv^oF '-'■'' PP'^) and then it is the amplitude 
c(0, 0) in Eq. ( |6.13| that contributes. Owing to its GIM 
suppression, the short-distance amplitude only gives a 
few-percent contribution to the full width dGaillard and 



Leel |1974b[ [Herrhch and Kalinowski[ |1992| |Ma and Pra 



mudita 1981 Pramudita 1988). Accordingly this decay 



is also dominated by long-distance dynamics. 

The dominant contribution, shown in Fig. |11[ is given 
by the weak transition of the kaon into a non-flavored 
pseudoscalar meson and its correspond ing decay into 
two photons (Ma and Pramudita 1981), the latter be- 



ing determ ined by the anomalous Lagrangian £wzw in 
Eq. (2.15). At lowest order in the chiral SU(3) expansion 
[C'(p'')J, only 7r° and rjg, propagate. However, the ampli- 
tude vanishes exactly due to the Gell-Mann-Okubo mass 
relation (iGell-Mannl 119611 lOkubol 119621) . Therefore, the 



decay starts at 0(jf) where the singlet 771 state is also 
included: 



Although there is still some discussion on the value of 
0, we take 9 = —20° arising in the large-JVc analy- 



ses ( Herrera-Siklody et al. 1998 Kaiser and Leutwyler 
1998[). ^ parametrizes the amount of SU(3) breaking 



(Donoghue et al. 1986) 



c 






(6.21) 



while p carries the information of the breaking of nonet 
symme try through the weak in teractions of the singlet at 
0(p2) ([Donoghue eTarj|l984|): 



3 (r;i| £l^^l-i|X°) 




(6.22) 



If SU(3) and nonet symmetries are exact we have ^ = 
and p — I. Finally, i^^, -Fi,g and F^i-^ are the decay con- 
stants of the pion, ryg and 771, respectively. The values 
of the symmetry-breaking parameters are still not pre- 
cisel y known. We take F^ JF^ = 1.34 an d F^, /F^ = 1.0 



from Kaiser and Leutwyler| ( 1998 ) [see also Feldmann and 



Kroll (2002|J. The SU(3) breaking parameter was esti- 
mated to be ^ ~ 0.17 (IDonoghue et ail 11986), but this 



value was challenged with the claim that it cancels with 
an additional s — )■ dgg contribution (He et al. 2003). 



Hence we will consider, conservatively, ^ ~ 0.0 — 0.2 
Finally, dominance of the pion pole seems to require a 
small breaking of nonet symmetry, p ~ 0.8 (Cheng 1990 
D'Ambrosio and Portoles[|1998a ) 



The decay width is given by Eq. (6.15) with the func- 
tion c(0, 0) instead of 6(0, 0) . If we consider the experi- 



mental determination (Nakamura et al. 2010) 



BR{Kl -^ 77) = (5.47 ± 0.04) x 10" 



(6.23) 



a value of 6* = —20° accommodates ^ = with a value of 
p ~ 0.7 while ^ = 0.2 requires p ~ 0.8. 



3. Ks 



-yi+r 



„(6) 



where 



(0,0) 



Fo 



1 - r2 'SI- r2) 



aFo{Gs-G27)F2{p,^,e), (6.18) 



(l + Ocos6' + 2\/2psin6' 



F„. 



2\/2-^sinf 

Fn 



3(1 -r^,) 



»)8 ^ no 

2V2pcos9- (l+,^)sin^ 



-sin0 + 2\/2^^cos( 



F, 



F,. 



(6.19) 



Vs ^ Vo 

Here 9 is the mixing angle between rjg and rji states: 



The amplitude for processes with one off-shell photon 
decaying into a lepton pair is given by 



^ ^ _ M^''(<Zi, 92) £* (gi) u(fc)7,«(fc'), 



92 



(6.24) 



with (72 



; + k'. At 0{p^) the amphtude 6(0, g|) in 
Eq. (6.13) is uniquely determined b y a one-loop calcu- 
lation of the Ki — )• 77* transition (Ecker et al. 1988 



Sehgal 1973). The dominant octet contribution is 



6(4) (0, 



li) 



-GsaFo{l~rl)H{z), 



(6.25) 



where z = qJ/Mk ^^'^ ^i^) i^ given in Eq. ( |A3| ). The 
result for the spectrum in g| is usually normalized to the 
two-photon width as 




cos a — sm 
sin 9 cos 




(6.20) 



1 dV 



J(l-.)3 



H{z) 



HiO) 



1 



Imn(z), (6.26) 
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with Imn(z) the electromagnetic spectral function asso- 
ciated to the lepton pair: 



spectrum in the variables Zi — qf/M'^, normalized to 
r^^ = T{Kl ->■ 77)- is given by 



-lmn(z) = —(1 + 2'^' 

IT StT 



1 -4- 



' (^ " 4r2) 



e^ 



(6.27) 
where r^ = mi/MK- T^^ = T{Ks -^ 77) is given by 



6^(0,0) in Eq. (6.141 [notice that H{0) = -i F{l/rl)]. 



Integrating the spectrum, one predicts the ratios 



r{Ks -^ j£+e-) _ I 1.2 X 10-2 (i = e), 

T{Ks ^77) ^1 2.8 X 10-4 (^ = m)- 



These processes have not been measured yet. 



4. Kl -^ 7^+£" 



(6.28) 



2 / a \2 



1 d^T 



T^-y dzidz2 Z1Z2 VStt 



|/(Z1,Z2)PA3/2(1,Z1,Z2) 



n l + 2T^\/l-4'f> (6-32) 



where /(zi,22) — c{qf , q2) / c{0 , 0) . A par ametrization 
suitable to analyze the data is provided by (D'Ambrosio 



et al. . 1998b I 



f{zi,Z2)\^^^ = l + ttDlP Yl 



z, - NPJMI 



'^°- n z,,-A'Pp/AP 



(6.33) 



i=l,2 



K 



As in the Kg decay, the normalized spectrum is given 
by Eq. (|6.26P with the ratio H{z)/H{0) substituted by 



fiz) 



c(0,0) 



l + hz + 0{z^ 



(6.29) 



Being of higher order in CHPT, our knowledge of the 
form factor f{z) is unfortunately rather limited. The de- 
cay mechanism in Fig. 11 generates an unambiguous con- 
tribution from the electromagnetic form factor of the n^ 
(r?, 77'), which amounts to a slope by = M'^/M'^ ~ 0.41. 
However, there are other possible contributions that have 



been analyzed only within explicit models (Bergstrom 



erar|[T983l[T990l [D'Ambrosio and PortolesHT997| |Ecker 



et al. 1990 Sehgal 1973), giving slopes in the range 



6c:i0.8-1.0. 

The experimental data can be fitted with a simple 



parametrization such as (D'Ambrosio et al. 1998b) 



^(^)|dip = 1 + "dip; 



which provides the slope }f 



Mj/Ml' 



(6.30) 



'K 



/MJ_ 



A re- 



cent determination by the KTeV collaboration ( Abouzaid 
et al. 2007b) for the electron case finds OfoiP = — 1.73 ± 
0.05 (6°^^ = 0.71 ± 0.02), while a^ip = -1.54 ± 0.10 
(6dip = 0.63 lb 0.04) was previously ex tracted from the 
muon mode ( Alavi-Harati et al. 2001c). 



5. Kl^ It i^ It ^2 

The dominant long-distance amplitude for this process 
is driven by the creation of two lepton pairs: 

AiTd 
A = -.^M'"'{qi,q2)u{ki)-f^v{k[)u{k2hMk2), 
I1I2 

(6.31) 
with qi = ki + k[. If CP is conserved only the ampli- 



The channels ^1 = ^2 = e and li = e, £2 = /i have 
already been observed experimentally. The first one is 
only sensitive to the paramet er amp due to the low invari - 
ant mass of the e+e" pairs (Alavi-Harati et al. 2001d[ ). 
Although the second channel should be more suitable 
to extract the parameter f3r,iri in practice that is still 



|20 



not feasible due to poor statistics (Alavi-Harati et al. 
2003b). Consequently, /3dip = is assumed and aoip is 
extracted. The most precise determination is obtained 
for the channel with muons giving q^dip = —1.59 ± 0.37, 
in good agreement with the results from K^ — >■ 7//+/i~. 
The world average is dominated by the Kl — > ^e~^e~' 
mode and is given by Qdip = — 1.69 ± 0.08. A complete 
lowest-order cal culation of QE D radiative corrections was 
also performed (Barker et al. 120031). Its impact on the 
slope ttDip could amount to a 15% correction, well within 
the present error. 

The high-g^ behaviour of the function ,f{zi,Z2) in the 
model of D'Ambrosio et al. ( 1998b) enforces a constraint 
on its parameters that is violated only in a very mild way, 
namely the sum rule 1 -I- 2q;dip + /3dip = 0. Hence, from 
the world average value of a^ip , we can estimate the phe- 
nomenologically elusive second parameter: /S^ip ~ 2.4. 



c. K^i+r 



by 



The most general amplitude for this process is given 



A{K -^ i+r) = u{k) {iB + C75) v{k'). (6.34) 



If CP is conserved the amplitude B (p-wave ^Pq) deter- 
mines K^ — )■ t^i~ while K2 — )■ t^t~ proceeds via the 
s-wave ^S'o amplitude C . The associated width is 

TiK ^ i+n = ^A (/3||i?|2 + |Cp) , (6.35) 



tude c((7i,(72) in Eq. (6.12) contributes. The invariant with /?£ — yjl — imj/M'^ 



K- 
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FIG. 12 Two-photon contribution to K° -> I'^T . 

1. A's -> I'^r 

The main contribution comes from the amp litude B^~f 
that gives the transition ^5 -> 7*7* -> l'^(.~ ([Ecker and 



Pich 1991) 



This two-loop ampUtude is finite because 
chiral symmetry forbids any CP-invariant local contribu- 
tion at this order. The resuh can be written as 



R's = 



T{Ks^(.+(.-^) 
T{Ks ^ 77) 

27r2|iJ(0)|2A'f 



(6.36) 



/, 



K 



<?,disp 



ilo 



abs 



where H{0) foUows from Eq. (A3) and /^^disp (^£,abs) in- 
dicate the dispersive (absorptive] parts of the two-loop 
diagrams. The i = e case is dominated by the absorp- 
tive contribution: /g^disp — 1-4, ^e.abs — —35. This gives 
-Rg — 7.8 X 10^^ corresponding to 

BR(if5 -^ e+e-) = 2.1 x 10"^", (6.37) 

to be compared with the recent bound by Ambrosino| 
eFaLl ( [2009b| : 

BTl{Ks -^ e+e-) < 9 x 10^^ (90% C.L.). (6.38) 

For £ — fi there is a slight dominance of the dispersive 
contribution, /u.disp — ~2.8, /„_abs — 1-2, which implies 
i?g = 1.9 X 10^6 and 



BR{Ks 



M /^ 



5.1 X 10"^^ 



(6.39) 



to be compared wit h the (almost 30-year-old) bound 
( [Gjesdal efal] |1973[ ) 

BR{Ks -^ ^i'^^i~) < 3.2 X 10"^ (90% C.L.). (6.40) 

There exists a small short-distance contribution to this 
decay through the CP-violating component of the s-wave 
amplitude C in Eq. (|6.35|), which in the SM was esti- 



mated for the muon case as (Buchalla et al. 1996 Isidori 



and Unterdorfer 2004 1 
Gf a{Mz) 



G 



sd 



TT sin 9 



V2m^FKlT^{V:,V^,)Y{xt 



w 



(6.41) 

where Y{xt) is given in Buchalla et al. (1996). Using 
the Wolfenstein parametrization, the short-distance con- 
tribution to the branching ratio is given by 



BR{Ks 



MA* Jsd 



1.4 X 10 



1.7 X 10" 



-12 



-13 



Kb 

0.041 


4 


A 
0.225 



(6.42) 



an order of magnitude short of the CP-invariant contri- 
bution. 



2. Kl 



£+£- 



If CP invariance is preserved this decay is given by the 
amplitude C in Eq. (6.34 1. It is convenient to normalize 
the rate to the Kl — > 77 mode: 



Ri = 



T{Kl^£+£-) 



T{Kl 

Vtt 



-re 



11) 
2 



(6.43) 



(1^^, 



displ 



|-P<?,absl ) 



The absorptive amplitud e gets contributions from several 
available on-shell states (iMartin et al. 1970) , but the 77 



intermediate state dominates completely (Fig. 12 1: 



'abs 



lit 



^1 - 

2p, ""yi + fi, 



(6.44) 



The dispersive part arises from the one-loop diagram 
Kl -^ 1*1* -^ f''~^IJ'~ and local CHPT terms to ab- 



sorb the loop divergence (Gomez Dumm and Pich 1998 



Isidori and Unterdorfer 2004 Knecht ei aL 1999a I: 



^,disp 



4/3, 



In^ 



1-/3^ 



-Lio 



I3i-l 
ft + 1 



3 In 



mi 



x(m), 



'We 
(6.45) 



where the local contribution can be split into its long- 
distance and short-distance components: 



x(Ai) = Xiiil^) + Xsd- 



(6.46) 



X^..y (/i) compensates the scale dependence of the one- loop 
amplitude, while Xsd accounts for the short-distance con- 
tribution sd — >■ 11^ fi~ . 

The Kl — > 1^-^ l-J'^ decay is well established experimen- 
tally. The measured rate (Nakamura et al. 
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R1 



(1.25 ±0.02) X 10" 



(6.47) 



appears to be nearly saturated by the absorptive contri- 



bution (|6.44|, i?^|abs = 1-195 X 10-5. The difference be- 



tween these two numbers provides an experimental mea- 
surement of the modulus of the dispersive amplitude and, 
therefore, up to a two-fold ambiguity, of the local term: 



x(Mp 



3.75 ±0.20 
1.52 ±0.20 



(6.48) 



The Kl —>■ e+e rate is then predicted to be i?£ = 
(1.552 ±0.014) X 10-8 or (1.406 ± 0.013) x 10-^, respec- 
tively. Both values are in agreement with the present 
experimental resul t Rl\exp = (1-65 ± 0.91) x 10-^ (|Am- 
brose et al. 1998), obtained with only four events 



In- 
cidentally, this is at present the tiniest branching ra- 
tio ever measured: BR{Kl — >■ e+e-) = 9^4 x IQ-^^. 
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( Ambrose et al. 



The corresponding branching ratio for t he muon case is 

BRJK l -^ ^l+^J.-) = (6.84±0.11) x 10 

2000| . 

At 0{p^) in the chiral expansion, the Kl — >■ 77 amph- 
tude is given by the tt", 77, 77' exchange mechanism shown 
The details of the weak transition vertex can- 



in Fig. 11 



Therefore, removing the short- 
d, Eqs. (le^iSl) to 0A^ can be 



eel out in the ratio R)^. 

distance contribution %, 

directly applied to the electromagnetic 27 decays of the 

neutral unflavorcd pseudoscalar mesons. Moreover, the 

local CHPT contributions are also the same. Thus, from 

the measured 7r° — )■ e~^e~ and rj — > ij.^ij.~ rates one can 



determine X-yyifJ-) (Gomez Dumm and Pich 1998), again 
up to a two-fold ambiguity. One gets x-yjiMp) = 0.8±1.0 
or-20.7±1.0from7r0 ^ e+e", a.nd x-yyiMp) = 5.5±0.9 
or —0.8 ± 0.9 from 77 — > i^^ n^ ■ The second tt'^ — >■ e^e'' 
solution is clearly excluded. The first one agrees very well 
with the negative 77 — >■ /i~''/i~ solution, but it is only 3.2ct 
away from the positive one. Keeping the two alternative 
possibilities and following the PDG average prescription, 
we get 



X77(^p) = 



-0.1 ±0.7 
3.4 ±2.3 



(6.49) 



to be compared with the theoretical estimate x-yy {M^ 



0.3±0.9 (Knecht et al. 1999a I, evaluated within lowest- 
meson dominance in the large-A''c framework. An al- 
ternat i ve theo r etical estimate was obtained ( D'Amb rosio| 
et al. 1998b Isidori and Unterdorfer 2004p perform- 
ing the 27 loop integration with the form factor (6.331, 
with anip = -(1 + /3dip)/2 = -1.69 ± 0.08 as deter- 
mined from the radiative decays. One gets in this way 
XiiiMp) 



3.3 ± 1.3, also compatible with (6.49). 



The short-dis tance contribution to K^ -^ ii^ IJ^ is well 
known at NLO ( [Buchalla et al.\ |1996[ ) and a NNLO eval- 
uation of the charm-quark contribution Ynl was carried 
out more recently (Gorbahn and Haisch 2006). The 



relative sign with respect to the long-distance ampli- 
tude can be fixed unambiguously in the large- TVc limit 



2004): 



(Gomez Dumm and Pich 1998 Isidori and Unterdorfer 



Xlf = 4.965x10^ [Re{\t)Y{xt)+Re{K)Y^^] 



= -1.82 ±0.04. 



(6.50) 



Adding the long-distance component (6.49), one gets 
x(Mp)_= -1.9 ± 0.7 or 1.6 ± 2.3, to be compared with 
Eq. (6.48). Although the present uncertainties are still 



too large to make a meaningful test of Xsd, a better un- 
derstanding of the long-distance amplitude could well un- 
cover new-physics contributions. 

The longitudinal polarization Pl of either muon i n the 
decay Kl — >■ A'^A'" is a measure of CP violation (|Pais 



and Treiman[ |1968b| |Sehgal[ |1969[ ). Within the SM, the 
main source for Pj^ is indirect CP violation due to K^ — 
K^ mixing. The longitudinal polarization arises from the 



interference of both amplitudes in Eq. (6.34 1: 

Nr - Nl 



Pl^ 



Nr + Nl 



^^Im(SC*), (6.51) 



where Nr, Nl are the numbers of outgoing fi~ with pos- 
itive or negative helicity, respectively and P is the full 
decay width. 

Within the SM, one expects ( [Ecker and Pich[ |1991[ ) 



|Pl| = (2.6±0.4) X 10^^ (6.52) 

This polarization has not been measured yet. 

D. K ^ TTf-y^*^ 

The general amplitude for K{k) — > T^{p)"fiqi)^{q2), 
A{K ^ 7r77) = e;(gi) 4(92) M'^''(fc, 51,92), (6.53) 
contains four Lorentz structures: 

A{z,y) 



Af^" 



Ml 



'iB{z,y) 
M^ 



(92V -91- 92 .gn 
(-fc.(7ifc.92 5^''-gi-g2fc''r 

+ k■qlq^k'' + k■q2k^'q•() 



C{z,y) 
Ml 

D{z,y) 

Ml 



e^'""'qipq2. 



(6.54) 



[e^'''"'{k-q2qip + k-qiq2p)k^ 
+ (fcP^-^7 + A;''£''"^'') fc„(7i^g27] 



where z = {qi + q2Y/Ml and y = k ■ {qi - q2)/Ml. 
Bose symmetry requires the invariant amplitudes A{z, y), 
B{z, y) and C(z, y) to be even in y, while D{z, y) is odd. 
In the limit where CP is conserved, A and B contribute 
only to Kl -^ t^^H, while C and D contribute to Ks — > 
7r'^77. All of them are involved in K~^ — > 7r"'"77. 

The double differential rate for unpolarized photons is 



dydz 



Mk 
297r3 



z" {\A + B 
' . 1 



\C? 



(6.55) 



y 



A(l,<,^) 



m 



\D\ 



The physical region is given by < \y\ < A^/^(l, rJ, z)/2 
and 0< z < (1 - r^)^. 

The processes K — >■ 7r77 have no tree-level contribu- 
tions of 0{p'^). At 0{p'^) the amplitudes B and D are 
still zero, since there are not enough powers of momenta 
to generate the gauge structure. Therefore, B and D 
arise only at O(p^). Notice that both B and D lead to 
contributions also for small z. 

The amplitudes for K -^ nj^* have a related but much 
more involved structure than with both photons on-shell. 
We comment briefly on these processes. 
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FIG. 

K+ 



13 Unitarity contribution from K'^ 



TV n TT to 



TT 77. 



1. K^ 



TT 77 



The leading contribution of 0{p'^) to the dominant 
octet amplitude for K'^ — > 7r+77 was determined in 



Ecker et al.\ ( |1988) 



(4) ^ GsMj^a 
^ ' 27rz 



[{z + l-rl)F{zlrl) 
+{z-\ + rl)F{z)-cz , (6.56) 



where F{z) is defined in Eq. (All and 
c = 1287r2 [3(L9 + Lio) + iVi4 - N15 - 2Ari8] /3 (6.57) 



in terms of the LECs in Eqs. (2.151, (2.161. Note that the 



loop contribution is finite and, consequently, the coun- 
terterm combination is scale independent. The subdom- 
inant 27-plet contribution to the amplitude A^^\z) was 
also determined ( Gerard et al. 2005 1 . In addition, start- 
ing at 0{p'^) but going beyond it. 



C(z) = 



GaW^a 



3z 



' Mk 



i{z-rl) 



(6.58) 



where r^ = MyjMx and T^o = r(7r° -^ 77). No- 
tice that the imaginary part of the -q pole has not been 
included. This amplitude is generated by the Wess- 
Zumino-Witten functional in Eq. (2.15) through the se- 
quence K~^ — >■ 7r"'"(7r'^, 77) — > 7r"'"77. The rj' contribution, 
though fairly suppressed by its mass, was also considered 
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within U(3) CHPT ( [Gerard et al 

Although there is no complete evaluation to O(p^), 
the most important contributions have been estimated. 
The unitarity corrections from if"*" — ^ 7r"'"7r+7r~, shown 
i n Fig . [T3I were determined by |D'Ambrosio and Portoles| 
( [1996 ). They contribute to both amplitudes A{z,y) and 
B{z), where A becomes y-dependent. Local non-resonant 
contributions were also studied and a naive chiral dimen- 
sional estimate indicates that they are small. Vector- 
resonance contributions to the 0{p^) LECs are intro- 



duced throu gh their effect on the amplitude B (jCohen 
et al 



1993): 



2Gi^Mia 



gjw^t 



lim Bv{z). 



(6.59) 




FIG. 14 Predicted spectra in the diphoton invariant mass at 
different values of c for K^ — >■ 7r"''77. 



The vector-exchange contribution to A can also be writ- 
ten in terms of the parameter ay , if we assume that these 
local amplitudes are generated through strong resonance 
exchange supplemented with a weak transition in the ex- 



ternal legs (Ecker et al 1990): 



Av = GsM]^aa+{2, + rl- z)/t:. 



(6.60) 



These vector-exchange contributions were estimated by 
D'Ambrosio and Portoles (1997) and found to be very 



small: a 



-0.2 ±0.3. 



The branching ratio of this decay and its spectrum 
were measured by Kitching et al. (1997): 



BR(ii:+ -^ 7r+77) = (1.1 ± 0.3) x 10" 



(6.61) 



Hence it is possible to determine the value of c in 
Eq. ( 6.56 ) . From the analysis of the normalized spectrum 



at C^^the value c = 1.6 ± 0.6 (xVdof = 0.9) is ob- 
tained. Including the unitarity corrections and neglecting 
local 0{p^) contributions, the fit improves (x^/dof = 0.7) 
yielding c = 1.8 ± 0.6. The diphoton invariant-mass dis- 
tribution is shown in Fig. [M] for different values of c. The 



spectrum observed by NA48/2 is shown in Fig. 15 



The electromagnetic penguin operators induce a CP- 
violating charge asymmetry in K^ —>■ 71^77. This asym- 



metry is estimate d to be tiny within the SM ( Ecker et al. 
19881 [Gaol 120031). 



2. Ks — >■ TT 77 

The leading contribution to this decay is given by 
K'-' — > 7r''(7r'',77) — > 7r"77, where the first part is dom- 
inated by the octet weak transition in the Lagrangian 
(2.16) and the second part arises from the Wess-Zumino- 
Witten anomalous term in Eq. (2.15). The am plitude 
starts at 0{p'^) but includes also higher orders (Ecker 
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FIG. 15 Spectrum for K^ — >■ 7r"'"77 in the diphoton invariant 
mass by NA48/2. From [Moralesl (|2008[) . 
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FIG. 16 Comparison of the normalized spectra in the dipho- 
ton invariant mass for Kl — >■ 7r°77 at Oip^) and 0{p^) in 
CHPT. The data are from |Lai efaj^^m^ . 



et al. 1987b I 






'■^ Mk 



3F„ (z - 



(6.62) 



keeping F.^ ^ F^ to account for chiral corrections. An 
amplitude D{z^ y) would arise at higher chiral orders. As 
in the K^ decay, the rj' contribution was also included 
by Gerard et al. (20051. It is seen, however, that the 



amplitude is dominated by the pion pole. 

In order to eliminate the overwhelming background 
from Ks "> tt'^tt", one restricts the kinematical region 
to z > 0.2. The chiral prediction gives 



BR{Ks -^ ^°77).>o.2 = 3.8 x 10" 



(6.63) 



that compares well with the experimental measurement 
(iLai et aL 1 120041) 



BR{Ks ^ 7T°ll)z>o.2 = (4.9 ± 1.8) x 10" 



(6.64) 



3. Kl 



Tv'^'yy 



The relevance of this channel extends beyond its own 
interest as it provides a CP-conserving contribution 
via the two-photon cut to the decay Kl — >■ 7r°i?+f~ 



( Donoghue et al. 1987 



Sehgal 



1993a 



Seh^ 



Ecker et al. 1988 Heiliger and 
1988) that competes with direct 



and indirect CP-violating contributions. 

The absence of 0(p'*) counterterms for this process in- 
dicates that the only contribution is a finite one-loop re- 
suit at this order ([Cappiello and D'AmbrosIo 1988 Ecker 



tude A{z) is 
AW(z) = 



etal. 1987b Sehgal 1990). The octet part of the ampli- 



GsMla 



K<- 



[{z - rl) F {z/rl) 
-{z-\-rl) F{z)\ , (6.65) 



where F{z) is given in Eq. ( Al ). This result gives rise to 
BiK{Kl -^ 7r°77) — 6.8 x 1 0~', significantly smalle r than 
the present PDG average (Nakamura et al. 2010) 



BR(i^L ^ 7r°77) = (1.27 ±0.03) x 10" 



(6.66) 



This indicates that higher-order chiral corrections should 
be considered. 

The unitarity co rrections from Kl — ^ 7r^7r~7r° con- 
tribu ting at 0{p^) (Cappiello et al. 1993[ Cohen et al. 



1993) give rise to both A and B amplitudes. The inclu- 
sion of the 2 7-plet contribution was carried out by |Cap- 
piello et al. (1993). As in the K^ case, those amplitudes 



also get local contributions from both vector exchange 
and non-resonant parts. From dimensional analysis the 
latter seem to be negligible, but the situation of vec- 
tor resonance contributions is still not settled. This is 
parametrized by Oy, defined as in Eq. (6.591, that was 
estimated theoretically to take the value ay = —0.7 ±0.3 
(D'Ambrosio and Portoles 19971. The measured spec- 



trum clearly 

ter: NA48 (JLai 



favors a non-zero value for this parame- 

= -0.46 ± 0.05 



et al. 2002) gets 



et al. 2008b) obtains 



while KTeV ^^ouzaid 

—0.31 ± 0.09. The predicted spectrum is compared with 

the measured distribution in Fig. 



16 



On the basis of the experimental situation at the time, 
Gabbiani and Valencia ( 2002 ) argued that the complete 
set of three 0{p^) counterterms is needed to describe 
both branching ratio and spectrum. However, the latest 



exper imental analyses ( Abouzaid et al. 2008b Lai et al. 



2002) find a satisfactory representation in terms of a^ 



only, supporting the dominance of vector exchange. 



A. K -> TT-yl+r 

These decays are obviously mediated hy K ^ 7777* 
and are closely related to the case with both photons on- 
shell. The study o f J^+ ->■ 77+7^+^" and Kl-^ 7r°7^+^- 
was carried out by Donoghue and Gabbiani ( 1997[[T998 ) 
and|Gabbiani|(|1999|). 
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The decay amplit ude h as a more complex Lorentz 
structure than in Eq. (6.54|. Nevertheless, at 0{p'^) these 
decays have the same properties as K'^ -> 7r+77 and 
Kl — >■ 7r°77, i.e., both have finite one- loop contributions, 
the charged channel has the same scale-independent 
counterterm c and the neutral one has no local contri- 
butions. Taking c = 1.8, one obtains at this order: 



BR{K+ - 
BR(A'l 



TT 76 e' 
.0 



^) = 1.4 X 10" 
7r"7e+e") = 1.0 x 10" 



(6.67) 



These predictions have to be compared with the recent 
experimental measurements (Batley et al.\ 2008a) 



BR{K+ -^ 7r+7e+e") = (1.29 ± 0.13) x 10"^ (6.68) 



and (Abouzaid et al. 2007a) 



BY{.{Kl -^ 7r°7e+e-) = (1.62 ± 0.17) x 10" 



(6.69) 



As can be seen, the agreement is rather good in the 
charged-kaon channel. 

Unitarity corrections from K — > tttttt, extending be- 
yond 0{p'^), were performed for both decays. We com- 
ment on them in turn. 

• K'^ — ^ 7r+7e+e^: Local 0{p^) counterterms were 
found to be negligible in K~^ — > 7r+77. However, 
neglecting them also in this case and using c = 1.8, 
one finds BR{K+ -^ n+'ye+e-) = 1.7x 10"^, which 
disagrees with the experimental result. Indeed, the 
analysis of the spectrum (for rrieej > 260 MeV) car- 
ried out in.Batley eiaL|(|2008a|) gets c = 0.90±0.45. 



Hence local contributions of 0{p^) and beyond 
might be non- negligible after all. 



K]^ — )■ 7r°7e"'"e~: Together with unitarity correc- 
tions, Donoghue and Gabbiani (1997) included ex- 
plicit contributions from vector resonances. Their 
prediction BR{Kl -^ 7r°7e+e") = 2.3 x 10"® is in 
slight disagreement with the experimental measure- 
ment quoted above. Abouzaid et al. (2007a) con- 

-0.46 (in 



eluded that employing a value of Uy 
accordance with Kl — >■ 7r'^77), the branching ratio 
goes down to BR{Kl -^ 7r°7e+e") = 1.51 x 10"®, 
in better agreement with the measured width. 



E. K 



■ni+r 



The FCNC transitions K -> ttI+I- {£ = e, /x) are dom- 
inated by single virtual-photon exchange {K — ^ t^J*) if 
allowed by CP invariance as in Ks and K^ decays. This 
contribution is CP violating for K^ —>■ 7r°£+£" and con- 
sequently this process has become a point of reference for 
studying the CP-violating sector of the SM. 



1. Ks, K^ 



Tii+r 



The amplitude for K{k) -^ 7r(p)^+(p+)£ (p_) is de- 
termined by an electromagnetic transition form factor 




FIG. 17 Unitarity K — >■ tttttt contribution to K^ , Ks — >■ TT7*. 



in the presence of the nonleptonic weak interactions 
p'Ambrosio et a/.l|1998al|Ecker ei aL||1987a|): 



A 



Gp a 

47r 



V{z){k+pYu,{p_)^^v,{p+), (6.70) 



where z = q^ /M'j^ and q = k — p. The spectrum in the 
dilepton invariant mass is given by 

dV Gla^Ml 



^K 



dz UniAn) 



p/^/l 



\V{z)\\ 
(6.71) 



where A = A(l, z, rj) and Arj < z < (1 — jv)^. 

Because of gauge invariance V{z) = at 0{p'^). The 
first evaluation at C'(p^) of the dominant octet amplitude 
was carried out bylEcker et al.\ il987a.\j: 



V+{z) 



Gh 



Gp 

for K"^ -> 7r='=7*, and 

Vs{z) -- 
for Ks -^ 7r°j*. Here 



<^{z) + <^{z/rl)+w+' 



Gg 
Gf 



[2^{z)+ws] 



$(z) = ^-- + -(i--)g(z), 
^ ' 18 3z 3z V a) ^ '' 



(6.72) 



(6.73) 



(6.74) 



wher e G{z ) is gi ven in Eq. ( A2 ) . In terms of the LECs in 
Eqs. (2.15) and (2.16), the local 0{p'^) contributions are 



647r2 

3 

32^2 



.-A^r5 + 3LS) + -ln^, 



WS 



(2iV[4 



N- 



157 



- In — ^ 
3 M^ 



(6.75) 



K 



Although a complete study to 0{p^) has not been per- 
formed yet, the unitarity corrections from K — ^ tttttt, 
shown in Fig. |17| were determined by |D'Ambrosio et al.\ 
(1998a). The full result can be decomposed in a polyno- 



mial contribution (linear in z up to this order) plus the 
unitarity loop corrections. 



V,{z)^a, + hz + Vr{z) 



ii^+,S), (6.76) 



38 



the latter being given by 

T.,r7r/ ^ _ Oj + Pijz - Zp) /rl 



1 



<^{zlrl) 



1 
(6.77) 



with zo = r^ + 1/3, ry = Mp/MK and 

i73, /3+-2(ei+e3-^3)> 



= /3i~^/33 + V3' 
4 



as = 



^73, 



Ps = o ^3' 



(6.78) 



in terms of th e pa rameters in Eq. (5.33). The 
parametrization (6.76) includes the 0{p'^) contributions. 



The pol ynom ial part incorporates the local counterterms 
in Eq. \Q.lb\ through a^^^ = G^IGpil/i - w+) and 
ag = — G8/Gi?(l/3 — ws)^ and tiny contributions to 
the slopes b] ' from the O(p^) kaon loop. 



Vi{z) in Eq. (6.761 is expected to be an excellent ap- 
proximation to the complete form factor of 0{p^). It 
only assumes that all contributions except the two-pion 
intermediate state can be well approximated by a linear 
polynomial for small values of z. The predicted rates can 
then be expressed in terms of a^ and hi in Eq. (6.76): 



BR 



K± 



BR'^i 



BR' 



Ks 



BR^, 



[0.15 - 3.31 a+ - 0.90 b+ + 60.51 a\ 
+ 16.36 a+h+ + 1.77 6^ ] x 10"^ , 
[1.19- 19.97a+- 6.565+ + 120.16a^ 
+ 69.42 a+b+ + 10.59 h\ ] x 10"^ , 
[0.01 - 0.55as - 0.176s + 43.76a| 
+ 11.83 asbs + 1.28 fe| ] x 10"^° , 
[0.07 - 3.96 as - 1.34 6s + 86.90 a| 



50.21 as5= 



7.66 fe| " 



X 10" 



(6.79) 



If flijfei ^ 1 as expected, the polynomial contribution 
dominates over the unitarity-cut loop corrections coming 
from K — ^ TTTTTT. In particular, due to the strong sup- 
pression of Ks — > Tr+vr^Tr*^, the latter are tiny for the Ks 
decay. Thus, these decays are very sensitive to the chiral 
LECs. 

The values of «+ and 6+ have been fitted from the 
K+ -^ T:+i+(r spectra. One finds a+ = -0.578 ± 0.01 6 
and 6+ = -0.779 ±0.066 for £ = e ( [Batley et am2009b[ ), 
while the muonic m ode gives a+ = —0.5 75 ± 0.039 and 
b+ = -0.813±0.145 (IBatley ei aLl|2011a|). The resulting 



experimental branching ratios are 



BR(/s:± -^ Ti^e+e-) 

BR(X± -^ TT^fl+fJ,-) 



(3.14 ±0.10) X 10"^ 

(9.62 ±0.25) X 10"^ (6.80) 



In Figs. 18 and 19 we show the spectra of K^ — > 7r+i?+i'^ 
{£ — e, ^) and the comparison with the theoretical pre- 
dictions for the chiral form factors (6.76 1 and for a linear 



form factor [without V^^'iz) in Eq. (6.76) 



For the Ks decay only branching ratios are available. 
By using the ratio bi/ai = 1/ry given by vector meson 
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FIG. 18 Spectra for K'^ — ^ n^e'^e in the dilepton invariant 



mass. The data are from Batley et al. (2009b I 
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FIG. 19 Spectra for K'^ — >■ tv^ fi'^ jj. in the dilepton invariant 
mass. The data are from Batley et al. (2011a I. 



1998a). it is found that 



dominance (D'Ambrosio et al 

1.06lg;2i fo'" the electron (IBatley et al 



\as\ 



\as\ = 1.54 



+0.40 
-0.32 



for the muon case ( Batley et al. 2004b I . 



2003 ) and 



Notice though that the vector meson dominance ratio 
bi/ai quoted above fails in the charged kaon decay. The 
measured rates are 



BR(ii:s -^ 7r"e+e") = (5.8 1^^) x 10"^ 
BR(i^s -^ TT°n+fi-) = (2.9 1};^) x 10-^ 



(6.81) 



where the large errors reflect the low statistics. 

Other parametrizations going beyond 0{p*) have also 
been put forward, focusing on the polynomial part of 
V{z). Friot et al. (2004b) assume a minimal narrow- 



resonance structure, while a related vector meson domi- 



nance approach is introduced in Dubnickova et al. ( 2008 1 . 



Unfortunately, the size of present data samples is insuf- 
ficient to distinguish between models. 

In addition to the dominant K — ^ nj* amplitude, the 
form factor V{z) receives a short-distance contribution 
from the operator Q^y in Eq. (2.9 1. Although negli- 



gible in the branching ratios and spectra, its interfer- 
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FIG. 20 CP-conserving contribution to Kl — >■ tt"^ 



ence with the long-distance ampHtude leads to a CP 



-> TT 



.± 



l+i~ (Ecker 



violating charge asymmetry in K^ . .. _ , ^ 

1988). This asymmetry is tiny within the SM 



et al. 



( [D'Ambrosio et a/.[[T998a||D'Ambrosio and Gaol|2002[ ) 



2. Kr 



^H+r 



Kl — > TT*' (£+f )j=i with the lepton pair in a vector or 
axial-vector state is CP violating and, accordingly, also 

Kl ^ TT 



■n'^^* . There are three main contributions to the 
^i'^i^ decay that, in principle, could be of the 



Ecker et al. 1988): 



same order of magnitude (Donoghue and Gabbiani 1995 



i. A direct CP-violating transition induced by the 



short-distance Lagrangian (2.5). The relevant con- 
tributions come from the operators Qyy and Qta 



in Eq. (2.9) and are proportional to ImAj. 



ii. An indirect CP-violating amplitude due to the 
K^ — K'^ oscillation, proportional to the CP- 
violating parameter e: 



A{Kl ^ 7r°^+r)|c 



= eA(/^5^7rVr). (6.82) 



7r°77 



iii. A CP-conserving contribution from Kl 

through 77 — >■ £+£~ rescattering, as indicated in 
Fig.[^ 

Let us first consider the CP-violating contributions 
(iBuchalla et al, 2003; iD'Ambrosio et al.\ |1998al iFlynn 



and Randalirri989b, .Heiliger and Sehgal[ |1993a[ [Isidori 



et al. 2004). Together with a vector contribution V{z), 
as in Eq. (6.70), from indirect CP violation and the Q^v 



amplitude, the Qja operator generates axial-vector and 
pseudoscalar terms: 



y'-(z) 
y-(z) 

A{z) 
P{z) 



±e [as + bsz], 

I ImAt f^ (z), 

^/ Za 

. 47r yjA 
I — 1= — 
^/2a 

. Sit in a 



V2c 






(6.83) 



where f^'^{z) are the if £3 form factors. In V'""^(z) 
we have neglected the tiny unitarity corrections from 
Ks — > tt+tt^tt'' and have made explicit a possible am- 
biguity in the relative sign with respect to the short- 
distance contributions. 

The pseudoscalar amplitude is helicity suppressed and 
can be neglected in the electron case. Using bs/ag = 
1/ry, the total CP-violating contribution to the rate is 



(Buchalla ei aL 2003 Buras ei aL 1994a I 



BR{Kl ^7r°e+e-)c 



(6.84) 



10"^2 X 



15.7 |asp± 6.2 las 



/ImAt 



Vio-4 



2.4 



/ImAt 
I 10-4 



For the mu on case (including t he dependence on m^), 
one obtains ( Isidori et al.\ [2004 ) 



BR(ifi^7rV+M")cPv = 



(6.85) 



IQ-^^ X 



3.7|asP±1.6|a5|(^ 



1.0 



/ ImAt 
UO-4 



In both cases the error in the numerical coefficients is 
estimated to be 10% — 20%. For an expected value of 
05 ^ 1, the largest contribution comes from the indirect 
CP-violating component, particularly in the case £ = e. 
Based on several assumptions on the matching between 
short- and long-distance dynamics, Buchalla et al. ( |2003 ) 
concluded that a positive interference between the in- 
direct and direct CP-violating components is the most 
natural setting. Using the values of \as\ extracted from 
the analyses of Ks — > 7r°^+^- and considering a positive 
interference, one finds 



BR{Kl -^ 7r°e+e-)cpv 
BR{Kl -^ 7rV+M")cPv 



(3.1 ±0.9) X 10"", 

(1.4 ±0.5) X 10-". (6.86) 



Now we consider the CP-conserving contribution from 
Kl — >■ 7r''77. Neglecting the y dependence of the 27 am- 
plitudes A, B in Eq. (6.54), the combination A + B gives 
|77)j=o and the amplitude B provides |77),/=2- However, 
the transition to \£^£~) is helicity suppressed (oc mi) for 



the initial J = state ( [Donoghue et al. . 1987[ ). Accord- 
ingly, for ^ = e only the amplitude B [arising beyond 
0{p'^) in the chiral counting] could give a non- negligible 
contribution, while for i = fi both amplitudes can be sig- 



nificant (Ecker et al. 1988 Flynn and Randall 



Isidori ei aL 2004 Morozumi and Iwasaki 1989). 



1989a 



Through a naive dimensional analysis of the ampli- 
tude B{z), [Ecker et aJ] ( [1988[ ) concluded that the CP- 
conserving contribution to Kl — >■ 7r''e'''e- was much 
smaller than the CP-violating one, making this decay a 
good testing ground of the SM. This was later confirmed 
by a two-loop calculation, including unitarity corrections 
from Kl — >■ 7r+7r-7r° and ignoring 0{p^) local contri- 
butions (Isidori et al. 2004), which found values below 
10-42 for the CP-conserving Kl — >■ 7r°e+e- rate. A local 
contribution to B{z), parametrized through the coupling 
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ay [see Eq. (6.59)] determined from K^ —> tt^j^ data, 



results in smaller values around 10"^'^. The contribution 
from B{z) turns out to be very small because of phase 
space and angular momentum suppression. Therefore, 
the Kl — > TT^^e+e^ decay is dominantly CP violating. 
The CP-conserving contribution to the muon channel, 

BR{Kl -^ 7rV+Ai")cPc = (5.2 ± 1.6) x 10"^^ (6.87) 

is of the same order but slightly smaller than the CP- 
violating component (6.86). The interference between 



CP-conserving and CP-violating amplitudes can then 
generate a sizable transverse polarization of the muons 
that could be within reach of the next generation of ex- 
periments (Ecker et al. 1988). 



The experimental upper bounds on these decays are 
gettin g close to the SM predic tions. At 90% C.L., we 
have (lAlavi-Harati et aZ. 1 12004) 



BR{Kl -^ 7r°e+e-) < 2.8 x 10" 



and ( Alavi-Harati et al. 2000b ) 



BR{Kl -^ 7rV+Ai~) < 3.8 x 10 



10 



-10 



(6.88) 



(6.89) 



The amplitude for K{pk) — !• ''^i{pi)''^2{P2)l{<l) 

A{K -> W7) = e;(<z) M''(g,pi,p2) (6.90) 

is decomposed into dimensionless electric {E) and mag- 
netic (M) components defined by 



M„, 



E{z^) Mjz,) , p ^ 

[ Zl P2m - Z2 Pl^ \ + ,,3 e^ypap^p'^q , 



M 



K 



Ml, 



(6.91) 

with Zi = q-pi/M'^ {i = 1, 2) and Z3 = Z1 + Z2 = E*/Mk, 
where E* is the photon energy in the kaon rest frame. 
Summing over photon helicities, there is no interference 
between both amplitudes and the double differential rate 
for an unpolarized photon is given by 



'^'^ - ''^ (ii?(zoi^+iM(z.)r 



(6.92) 



dzi dz2 (At:)^ 

X [ziZ2 (1 - 2(zi + Z2) -r\- r^) - r\zl - r^zl] , 

with fj = MttJMk- The electric amplitude can be de- 
composed as the sum of inner-bremsstrahlung (IB) and 
direct-emission (DE) components, E = E^^ + E^e, while 
the magnetic amplitude contains DE only. Ei-b is deter- 
mined by the Low theorem ( Low , [1958 ) : 



Emizi) 



MkZiZ2 



A{K -)-TriTr2). 



In general, E'ib dominates the photon spectrum for low 
values of the photon energy (in the kaon rest frame), but 
it is obviously absent for Ks, Kl —>■ 7r°7r°7. 



The DE amplitudes explore the electromagnetic struc- 
ture of the hadron interaction. They are customarily de- 



composed in a multipole expansion (D'Ambrosio et al. 
19921): 



E^E = Eiiz^) + E2{z3) ■ zi2 + O {ZI2) , 
M = Miizs) + Ahizs) ■ zi2 + O {ZI2) , (6.94) 

with Z12 — Zl — Z2. This decomposition is particularly 
relevant for initial states with definite CP (neutral chan- 
nels), as even electric and odd magnetic multipoles are 
CP odd while the others are CP even. Hence CP in- 
variance would require the amplitudes Ei , M2 , E^, . . . to 
vanish for Kl — ^ 7r+7r~7, whereas Mi, £^2, M3, . . . vanish 
for Ks — >■ 7r+7r^7. For K^ — )■ 7r^7r°7 all of them are al- 
lowed. Moreover, the limited phase space (|zi2| < 0.17) 
indicates that higher-order multipoles are suppressed. 

The decays K -^ Tr{pi)Tr{p2)i'^{k+)i~{k-) where the 
lepton pair is produced by an off-shell photon have am- 
plitudes 

AiK ^ Tn:i+i-) = ^V^ Mk-h''vi{k+), (6.95) 



where q ~ k^ + k^ and 






SliZi) 



M, 



K 



-P2m 



M{zi) ^ p ^ 

^fj.upijPiP2q ■ 



Ml, 



(6.96) 



1. K^ 



TT TT 7 



The IB contribution to this decay is suppressed by the 
A/ = 1/2 rule because it is proportional to the K~^ — >■ 
tt'^tt'^ amplitude, a A/ = 3/2 transition. However, this is 
still the largest contribution to the rate. The dominant 
DE contribution is 0{p'^) in CHPT ([Bijnens et oL 



1992 



Ecker et al. 1992). At this order, the electric amplitude 
has a finite one-loop contribution of which the octet part 
turns out to be tiny (D'Ambrosio and Isidori 1995; Ecker] 



et al. 1994): 



E„„„ — 






{Ml - Ml) [Kk{z+) + hKr^{z+)] 



where the functions hij{z) are given in Eq. (A4). A size- 



(6.97) 



able 27-plet contribution to the DE electric amplitude, 
from an enhanced two-pion loop, has been recently iden- 



tified (Mertens and Smith 2011). There is also a local 



scale-independent contribution 






{Nu - Afi5 - A^i6 - iVir) (6.98) 



in terms of the LECs in Eq. (2.16). Model-dependent 



(6.93) estimates ([D'Ambrosio and Portoles 1998b Ecker et al. 



1993) indicate that the vector meson contributions to 

i?^t are suppressed whereas an axial- vector meson ex- 
change piece remains. This amplitude could be mea- 
sured from the interference with the IB contribution. The 
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0{p'^) magnetic amplitude has a reducible component 
from the Wess-Zumino-Witten anomalous Lagrangian in 
Eq. ( |2.15 1 and a direct non-anomalous contribution from 
odd-intrinsic-parity operators in Eq. (2.16): 



m(4) 



' ^^ ^ [1 - 167r2 (37V29 - 7V3o)] . (6.99) 



27r2 F, 







A factorization estimate of the couplings fD'Ambrosio 



and Portolcs^ ^1998bJ indicates that the different contri- 
butions are of the same order and interfere co nstructively. 
Some 0{p^) analyses were also performed (D'Ambrosio 
and Isidorij |1995[ |Ecker et 11] |1994| but the lack of 



knowledge on the LECs at this order makes it difficult to 
assess their relevance. The role of form factors induced 



by vector meson exchange was also explored ( Cappiello 



and D'Ambrosio 2007 D'Ambrosio and Gao 2000). 



The differential cross section is usually written in terms 
of T^, the kinetic energy of the charged pion in the kaon 
rest frame, and of W"^ = {q ■ Pj<-)(g •p-i-)/(M^+M^), fa c- 



torizing the IB contribution ([D'Ambrosio et of. 1992), 



'2-r± 



d^r 



'2-r± 
IB 



d'T 



dT* dW^ 



dT: dW^ 



[1 + 2 cos {Si -Sl±: 



(Yi 



Yh 



)W^ 



YeW 
(6.100) 



where r+ (r~) corresponds to the decay of A'+ {K^). 
Here is a possible CP-violating phase, 6\ are the strong 
rescattering phases for a final tttt state of isospin / and 
orbital angular momentum L Ye, Ym are the DE elec- 
tric and magnetic amplitudes, respectively, normalized 
to K^ — >■ 7r='=7r°. It is usually assumed that Ye and Y]\i 
are nearly independent of T* and are therefore constants 
over the Dalitz plot. The linear term in W^ corresponds 
to the interference between the amplitudes -Ede and Eib ■ 



A recent measurement by NA48/2 (Batley et al. 2010b) 
establishes that Ym /Ye = —11 ± 3, showing that among 
the DE amplitudes the magnetic contribution dominates. 
Charge asymmetries in this channel are interesting 
observables for CP-violating effects ([D'Ambrosio and 



Isidori 



1998|). A bound on the fully integrated asym- 
(r+ - r-)/(r+ -l- r-) has recently been 



metry 5V 

published ([Batley et al.\ 



90% C.L. translates into 



2010b[ ): \5T\ 
< 0.56. 



smt 



< 1.5 X lO"'' at 
The asymmetry 



constructed with the partially integrated decay widths 
(Colangelo et al. 1999) is also of interest: 



dA 



w 



dT+/dW^-dr-/dW^ 
dW^ ~ dT+fdW^ + dT-fdW^' 



(6.101) 



Batley et al.\ ( |2010b[ ) obtain Aw = (-0.6 ± 1.0) x lO^^. 



2. Kr 



TV TT 7 



The IB component is proportional to the amplitude 
for Kl — > TT+TT^. Consequently, it violates CP and it 
is rather small. Assuming CP invariance, the dominant 
contribution is given by the Mi amplitude because E2 is 





FIG. 21 Diagram (a) has anomalous and resonance- 
dominated contributions to Kl — >■ 7r^7r~7. Diagram (b) gives 
a direct weak transition in terms of the couplings in Eq. (2.161. 



suppressed by phase space. These features make this de- 
cay very suitable to study th e DE magnetic amplitude 
([D'Ambrosio and Gao 



2000 D'Ambrosio and Isidori 



1995 1 [D'Ambrosio and Portoles[|1998a; Ecker et aL[|1992 



1994[ [Lin and Valencia! |1988[ |Picciotto[ [1992| ) 

The leading E2 amplitude arises at ©(p*) in CHPT. 
There is no local contribution and therefore E2 is given 



by a finite one-loop result (D'Ambrosio and Isidori 1995 
Ecker et aLl|1994[): 



E, 



(4) 



-Kk{z+) 



Ml) [Kk{z-) 
hK7^{z+)], 



hKriiz^) 

(6.102) 



with hij{z) given in Eq. (A4). However, as commented 
above, this contribution is even smaller than the CP- 
violating IB component: l^it^p/E'iBl -- 10^2 The27-plet 
contrib ution to the E2 amplitud e has also been deter- 
mined ( [Mertens and Smith[|2011| ). 

The Ml multipole at lead ing 0{p'^) is gener ated by a 
constant local contribution (jEcker et al. 1994) in terms 



of the LECs in Eq. (|2.16|) 



16 e ^«^^- 



Fn 



{N29 + Nzi). 



(6.103) 



The magnetic amp litude was also studied at 0{p^), 
including both local (D'Ambrosio and Portolls} 1998a[ 



Ecker et al. 1994 Lin and Valencia 1988 ) and one-loop 



contributions (D'Ambrosio and Portoles 1998a). A piece 



of the local term is completely fixed by the anomalous 
WZW Lagrangian in Eq. (2.15), as shown in Fig. 21 v. 



27r^Fo 



(6.104) 



in terms of 9, ^ and p defined in Eqs. (6.20), (6.21) and 
(|6.22|), where 



i^i 



1 



1 



3(l-i) 



fl 



F^ 



■ m 



' C) cos 6* -I- 2\/2p sin 61 
d-V2 



F, 



no 



sinfl 
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F^ 
F 



Vs 



^pcosO- (1 + sin 6* 



sin 6* + \/2 



F, 



Vo 



cost 



(6.105) 
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This expression includes higher orders beyond 0{p^). 

The first Z3 dependence comes with resonance dom- 
inated local counterterms that depend both on strong 
vertices (with a weak transition in the external kaon 



measurement of the equivalent of the parameter qmi in 



Eq. ( |6.107[ ) for the amplitude Ei 
find]^i7r< 0.21 (90% C.L.). 



Abouzaid eral](|2006b[) 



leg) and pure weak vertices (D'Ambrosio and Portoles 



1998a). A similar situation arises in the one- loop ampli- 
tude MiLp = M^zw+Mn- where the first term originates 
from the WZW Lagrangian with a weak transition in the 
external kaon leg, and the second is determined by ver- 
tices of the weak Lagrangian (2.16). Unfortunately, due 



3. Ks 



TV TT 7 



to the presence of weak LECs, it is not possible to provide 
a model-independent theoretical prediction. 

The phenomenological analysis was originally carried 
out with only the magnetic amplitude and with a p-pole 
dominated form factor (Lin and Valencia 1988): 



Contrary to the K]^ decay, this channel is dominated 
by the IB amplitude because there is no suppression in 
the Ks — > 7r+7r~ decay. From the remaining contribu- 
tions, Ml and E2 are CP violating. Therefore, the dom- 
inant DE amplitude is E i with a leading contribution of 
C'(p^). The octet piece (D'Ambrosio and Isidori |1995 



D'Ambrosio et al. 1993 Mertens and Smith 201 ID con 



tains a scale-independent local amplitude 



G&Ml ~ 
M = e ^A# ffMi 



1 



ai/a2 



E. 



(4) _ 



4ze 



GgA'fi 



27r2Fo ^^"^ V^ ' (Af2-Af|) + 2Mi^3. 

(6.106) 

Both p arameters were determined by [Abouzaid et al.\ 
( |2006bD : I^Afil = 1.20 ±0.09 and 01/02 = (-0.738 ± 
0.019) GeV^. A model- independent parametrization is 
provided by a power series expansion in the variable z^: 



Fo 



[Nm - A^i5 - iVi6 - N, 



17J 



(6.110) 



with the same combination of LECs that appears in 
K^ -^ 7r*7r°7 in Eq. (6.98). In addition, there is a finite 
one-loop contribution 



E, 



(4) 



. GsMk ,2 



Ml)[^AK,{-z:i) 



M = e 



G^Ml 
2tt^Fo 



^9 Ah (l + c^s). 



(6.107) 



+ Kk{z+) + hK^{z+) + {z+ ^ z_)] . (6.111) 



The values obtained from the expression (|6.106|) corre- 
spond to I^Mil 



1.30 ± 0.12 and c 



-1.74 ± 0.08. 



The 27-plet amplitude has been determined by [Mertens] 
and Smith (2011). Depending on the unknown combina- 



tion of LECs in Eq. (6.110), the ratio of IB and DE con- 
It was noticed ( [Ecker et al.\ |1994[ ) that the slope is tributions amounts to LiB/rDE ^ 10^ and IB therefore 



rather large compared with naive expectations based on 
resonance saturation. Indeed, a quadratic approxima- 
tion also shows rather large values for both linear and 
quadratic slopes ( Alavi-Harati et al. |2001e). Although 



dominates the experimental branching ratio (Nakamura 



et al. 2010) 



W£l{Ks -^ 7r+7r-7, E* > 50MeV) = (1.79±0.05) x 10" 



the value of c can be easily accommodated for a reason- 
able range of the parameters involved in the theoretical 



(6.112) 



analyses (D'Ambrosio and Portoles 1998a Ecker et al. 



1994 ) , the predicted photon spectrum is slightly at vari- 



4. K -^ Tnr£+e~ 



ance with the experimental determination. The ratio be- 
tween the width due to the DE magnetic amplitude and 
the total width was also measured^ToB/XTDE + Tib) — 



0.689 ±0.021 for E^ > 20MeV ( [Abouzaid et aZ.[|2006b ). 
The full branching ratio is 



BR{Kl -^ n+TT--f) = (4.15 ± 0.15) x 10"^ (6.108) 

Both IB and Ei amplitudes are CP violating. It is 
usual to parametrize this violation through 



\V+- 



J<l>+- 



A{Kl -^ TT+TT 7)|lB+_Ei 



= ^4 



A{Ks 

-7 -"n^ 



TT+TT 7) 



The lepton pair in this decay is produced by an off- 
shell photon and in consequence the dynamics is that 
of the K — > TrTr7* process. The widths are smaller 
than those for a real photon and only the decays 
Kl, Ks — » TT+TT "e+e" have been observed e xperimen- 
tally ( [Abouzaid et aL[[2006a[ [Lai et al.\\2m^ . This de- 
cay provides a useful tool to look for the CP-violating 
interference between electric and magnetic amplitudes 
because the lepton plane furnishes a measurement of 
the photon polarization vector. Bounds on the decays 
Kl — » TT°TT°£+.^~ have also been determined for £ = e 



(Alavi-Harati et al. 2002) and i = /j, (Abouzaid et al. 



where r]^_ is defined in Eq. (5.13). The experimen- 
tal determination of these quantities is already rather 
old ([Matthews et al.\ [19951 [Ramberg et al.\ [1993|: 



(6.109) 2011b). As the /z^/x case is strongly suppressed by the 
available phase space, we only consider e+e~ production 



(2.36 ± 0.06) X 10 



(44 ± 0.4f~aEd 



K- 

dictions ([Tandean and Valencia 2000[). More recent is a 



/e < 0.3 at 90% C.L., in accordance with SM pre- 



The decay amplitude is given by Eqs. ( 6.95 ) and ( 6.96 1 . 
A first study of the decay Kl -^ Tr+TT~e+e~ with con- 
st antekctricaiidmagneticamplitudes was carried out 
by Sehgal and Wanninger ( 1992 ), while the calculation at 
leading 0{p'^) of the electric amplitudes (keepin g a con- 
stant magnetic contribution) was performed in Elwood 
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et a^.| ( |l995[ ). The complete results up to C(p^) for both 
e~ and K~^ 



K 



TT^TT 



TT'^Tr'^e'^e 



can be found 



Pichl (2001) 



i. Kl — > TT+TT e+e 

For this decay £2 = —£i{pi ^ P2) and we will only 
quote the amplitude £1 . The leading tree-level con- 
tribution to the electric amplitudes is provided by 
j{0 _ j^o mixing with an ofF-shell photon radiating 
off the pion legs. Therefore, it is proportional to 
the indirect CP-violation parameter e: 



^1 



(2) 



-iiee 



GsMkFq 
2M'^zi + q- 



{Ml -Ml). 



(6.113) 



At 0{p^) a local scale-dependent contribution 



arises, 



^1 



(4) 



2ie 



G^Mk 
3Fn 



q^ [Nl^ - Nl^ - 3{N._ 



16 



N17)] , 

(6.114) 
together with a divergent one-loop amplitude 
(Pichl 2001). At this order, the magnetic ampli- 



tude is constant and coincides with the case of a 
real photon in Eq. ( |6.103[ ). In |Pichl| ( |200T| ) the mag- 
netic amplitude was modified by including a p-pole 
parametrization analogous to (6.106). 



TT+Tr'^e^e 



For this channel there is no relation between the 
two electric amplitudes as in the Kj^ decay. It is in- 
teresting to observe that at leading 0{p^) the octet 
contribution vanishes in the isospin limit and only 
the 27- plet represen tation gives a non-zero contri- 
bution (|Pichl[ 2001[). At 0{p^) there are local scale- 



dependent amplitudes 



(4) _ 



■--let 



^(4) 



C M 
-I e -^-^ [-GM^z^iNu - A^i5 - N,e - iVir) 

-4q'iNl,-N{,)], 

-» e ^^1^ [6M|.Z2(iVi4 - A^i5 - A^i6 - N^) 

-2q'{Nl, + 2NI,) + 6q'{Nl, - N„)] , (6.115) 



reducing to the scale-independent result ( 6.98[ ) for 
an on-shell photon. There is also a divergent one- 
loop contribution at 0{p'^). As in the previous case, 
the magnetic amplitude is constant and coincides 
with the case of a real photon in Eq. (6.99). 



The branching ratios of K]^ and Kg decays were mea- 
sured by N A48 (iBatley et al.\ 1201 Ibl iLai et al.\ 120031) : 



BR{Ks -> 7^+7^"e^ 



-) 



(3.08 ±0.20) X IQ-^, 
(4.93 ±0.14) X 10-^ 

(6.116) 



For the K]^ decay Pichl (2001) studied the impact of 



the weak LECs appearing in the amplitudes. However, as 



pointed out before, the main interest of this process is the 
study of the interference between electric and magnetic 
amplitudes provided by the measurement of the photon 



polarization ( Ecker and Pichl 



Heiliger and Sehgal[ 1993c Sehgal and van Leusen 



2001 



Elwood et al. 



2000 Sehgal and Wanninger 1992). As we have seen 



1996 



1999 



before, the electric amplitude of K]^ — >■ 7r"'"7r~e"'"e~ is 
generated by the indirect CP- violating mixing of K^ and 
K'^ . Its interference with the CP-conserving magnetic 
amplitude produces an asymmetry in the distribution of 
the angle (j) between the e+e" and tt+tt" planes in the 
kaon center-of-mass system: 



Ac 



Jo del) ' 

^/2 dr 

del) ' 



rr 
JO 



' Jn/2 d4> ' 
' + Jtt/2 d$ ' 



(6.117) 



This observable was measured^ by Abouzaid et al. 



(2006a) as A 



(14 ±2)%. Ecker and Pichl (2001) 



observed that the theoretical prediction of Acp with the 



leading electric and magnetic amplitudes in Eqs. (6.103) 



and (6.113) is short by almost a factor two compared 
with experiment. With the inclusion of an energy de- 
pendence in the magnetic amplitude along the model in 
Eq. (6.106 1, agreement with the experimental determina- 



tion is restored. 

G. Other decays 

1. K'^ — >• 777 

The CP-preserving amplitude for this process has 
parity-violating and parity-conserving contributions: 



A{K'^^^^^)=el{q^)e*p{q2)e*^{q^)[ 



M^^-^ + M?^^ 



PC 

(6.118) 

The first contributes to K^ — > 777 while the parity- 
conserving amplitude governs Ks — >■ 777. The Lorentz 
structure of both amplitudes is rather complex and can 
be fou nd in Heiliger et al. (1994) and Ho and Tandean 
(2010). Both amplitudes are strongly suppressed by an- 



gular momentum. This is due to the fact that any two 
of the three photons in the final state must have J > 2 
as a result of gauge invariance and Bose symmetry. 

Estimates of the branching ratios have recently been 
carried out (Ho and Tandean |2010 ) by using naive di- 
mensional analysis to predict the order of magnitude of 
the LECs in the appropriate chiral Lagrangian of 0{p^^) 
contributing at leading order: 



10"^^ < BR{Kl 
IQ-^^ < BR{Ks 



777) < 10-", 

777) < 10"^^. 



(6.119) 



^ The analogous asymmetry for the Ks decay is compatible with 
zero (IBatley et aL||2011b|l. 
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There is a recent experimental upper bound on the first 
process ( Tung et al.\ 2011): 

BR(ifi ^ 777) < 7.4 X 10"^ (90% C.L.). (6.120) 



2. Kl 



'IT 



£+£- 



This decay is relevant in order to det ermine the back - 
ground subtraction for Kl -^ 7r'^£+£~ (Greenlee 19901. 



Its theoretical description as a bremsstrahlung radiation 
off the leptonic legs in the decay Kl -^ ^t^t~ appears to 
be in agreement with the experimental branching ratios 
(lAlavi-Harati ei aLl|2000al|2001a|): 



^'K{Kl -^ 77e+e- E* > h MeV) = (5.8 ± 0.4) x 10" 



BR(A'i -> 77M+/^' 



m^^ > 1 MeV) 



(lOJ 



X 10" 



3. Kl 



■ "fUV 



This decay exhibits an interesting interplay between 
long- and short-distance contributions. Both of them 



have been considered (Richardson and Picciotto 1995): 



1) A long-distance amplitude consisting of two pieces, 
a weak transition in the Kl leg followed by constituent 
quark triangle diagrams, and a resonance-dominated am- 
plitude; 2) A short-distance contribution given by ds — > 
^vV. The results obtained within this approach give 



BR(i^i^7i/17)|id - 10- 
BR{Kl ^ -fiyV)\,d - 10- 



(6.122) 



showing that the process is short-distance dominated 
(iMaand Okadal|1978|. 



The CP-conserving and CP-violating short-distance 



contributions were estimated to be ( Geng et al. 
BR{Kl ^7'^'^)lcpc 



2000) 



1.0 X 10-^3, 

1.5 X 10"l^ (6.123) 

somewhat at odds with the results mentioned above. 
Given the small size of these branching ratios, this pro- 
cess is a good benchm ark to look for effects beyond the 
SM ([Jiang et a/. [[20031). 



4. Ks -^ etK4^2 

If CP is conserved, these decays proceed via 

Ki-,^*j*^£+e^e+q. (6.124) 

To leading order, ©(Ggp"*), only one-loop diagrams with 
pions contribute to the amplitudes. The branching ratios 
were estimated by Birkfellner (1996): 

BR(7^s ^ e+e~e+e") = 7x10"", 

BR{Ks ^ e+e-ii+^i-) = 8 x 10"^^ 

BR{Ks -^ ^i+pL-pL+^-) = 1 X 10"". (6.125) 

No experimental limits are available. 



5. Kl,Ks ^7r°7r°7 



The absence of charged particles in these decays im- 
plies that they do not have an IB contribution. More- 
over, in the CP limit the amplitude for the Kl decay 
is purely electric and the one of Ks is purely magnetic. 
In addition, Bose symmetry requires both multipoles to 
be even. A consequence of all these features is that the 
branching ratios are very small. They have still not been 
observed experimentally (Abouzaid et al.\ 2008c I: 



BR{Kl -^ 7r%°7) < 2.43 x 10""^ (90% C.L.). (6.126) 

For the Kl decay the leading multipole amplitude E2 
is at least of 0{p^), because the 0{p'^) oct et contribu- 
tions, both lo cal and loop pieces, vanish (Funck and 
Kambor[[T993l). The 27-plet contributions to the 0{p'') 



(6.121) one-loop am plitude give a tiny addition ( [Mertens and 



Smith 
10- 



-10 



2011 1. A naive chiral dimensional estimate of the 



ocal contributions predicts BR{Kl 



r-O^O. 



(D'Ambrosio and Isidori 1995 Ecker et al. 1994 



fjL 



Mertens and Smith 2011). However, other esti mates pre- 
dict a larg er branching ratio ~ 10~* (Heiliger and Seh- 
gal 1993b ). The Ks decay has not been studied in detail 



but it is expected to give a much smaller branching ratio 
([Heiliger and Sehgal[ [T993b| . 



6. Kl^Ks 



■K^ix'^^'y 



The decay Kl — » 7r°7r°77 was s t udied at leading or- 
der in CHPT ( [Dykstra et~ai] [1991[ [Funck and Kambor 
1993). It was found to be dominated by the pion-pole 



contribution, showing a strong dependence on the cut on 
the photon energies. 

The Ks decay has no local contribution at O(p^). The 
amplitude at this order is therefore given by a finite one- 



loop piece (Funck and Kambor 1993) that provides an 



unambiguous prediction: 

BR{Ks -> 7r°7r°77, m^^ > 20 MeV) = 4.7 x 10"^ 

(6.127) 
There are no experimental bounds for these decays. 



7. K ^ 37r7 

The experimental status of if — >■ 37r7 decays is still 
rather meager. Only the two channels with a charged 
kaon have been detected experimentally. None of the 
decay modes of a neutral kaon have been seen. 

Like in many other radiative K decays, there is 
both an electric and a magnetic amplitude. The mag- 
netic contributions are severely suppressed for all chan- 



nels (D'Ambrosio et al. 1997), especially for the ones 



where the non-radiative decay amplitudes are dominantly 
octet: K~^ — ^ 7r°7r°7r+7, K^ — ^ 7r+7r"'"7r~7 and Kl — >■ 
7r+7r"7r°7. For these channels, the rates can be predicted 
with great accuracy by using the concept of "general- 
ized bremsstrahlung" (D'Ambrosio et al. 1996). In this 
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way, optimal use can be made of available experimen- 
tal information on the non-radiative amplitudes. How- 
ever, as long as the minimal photon energy E"™'" is small 
(i?""" < 40 MeV), the difference between "generalized 
bremsstrahlung" and the leading Low contribution (Low[ 



1958 ) will not be accessible experimentally in the foresee- 



able future. The genuinely radiative contributions (direct 
emission) would only matter for large 
the expense of the number of events. 



emission) would only matter for large i?""", of course at 



The situation is different for the transition ^5 — >■ 
'Tr~TT^j. At lowest chiral order, the amplitude can 
only proceed through a A/ — 3/2 transition (via 
bremsstrahlung) and it is therefore suppressed by the 



TT+TfO^ 



the NLO 



A/ = 1/2 rule. Similarly as for K~^ 
contributions generated by octet operators become rela- 
tively more important. The electric amplitude is domi- 
nated by NLO weak LECs, with a large theoretical un- 
certainty (D'Ambrosio et al. 19971. 



The comparison between theory and experiment is dis- 
played in Table [VlTl For the two charged K decays and 
for Kl — ^ Tr+TT^TT^T, the theoretical accuracy is only lim- 
ited by the precision with which the parameters of the 
non-radiative amplitudes are known. The branching ra- 
tio for Kg — > 7r"'"7r~7r°7 has a much larger theoretical un- 
certainty. Experimental information on this mode would 
therefore be especially interesting, but the expected rate 
is unfortunately very small. 



VII. CONCLUSIONS AND OUTLOOK 

Kaons have been at the center of many ground- 
breaking developments in particle physics, which are 
worth repeating here. From the introduction of internal 
"flavor" quantum numbers (strangeness), to parity vio- 
lation, meson-antimeson mixing, quark mixing, CP vi- 
olation, and the suppression of flavor-changing neutral 
currents, kaon physics has played a key role in the con- 
struction of what we now call the Standard Model. 

In this review we have summarized the status of the- 
oretical predictions for all SM allowed kaon decays with 
branching ratios greater than 10~^^. Theoretical pre- 
dictions rely on the weak operator product expansion to 
organize the short-distance dynamics, the renormaliza- 
tion group to evolve the relevant operators down to the 
hadronic scale, and non-perturbative methods such as 
CHPT and lattice QCD to deal with low-energy strong 
dynamics. The accuracy of theoretical predictions ranges 
from excellent to just fair, depending on whether the de- 
cay amplitude is dominated by short- or long-distance 
effects. In all cases we have compared theoretical predic- 
tions to existing data and we have discussed the impact 
of future measurements. 

It is remarkable that to date kaons still provide unique 
information on physics at vastly separated energy scales: 



• On one hand, kaons are an excellent probe of low- 
energy strong interaction dynamics {E ~ ^qcd)- 
While many kaon decays are by themselves a valu- 
able testing ground for chiral effective theory tech- 
niques, specific modes such as K — > tt'kH.v and 
K — > TTTTTT allow one to extract reliable informa- 
tion on low-energy tttt scattering, thus testing our 
understanding of purely strong dynamics. 

• On the other hand, kaon decays encode informa- 
tion about flavor dynamics at the electroweak scale 
E ~ O(10^~^) GeV through their sensitivity to 
virtual exchange of heavy SM particles {W , Z , t) 
and possibly non-standard particles and their fla^ 
vor structure. Theoretical predictions for K^ — K^ 
mixing, short-distance dominated rare decays, and 
semileptonic branching ratios have reached a level 
of accuracy that already puts very stringent con- 
straints on extensions of the SM, in ways often 
complementary to other flavor observables and elec- 
troweak precision tests. Here we have not discussed 
in a systematic way the impact of kaon decays on 
extensions of the SM, which can be found in re- 



views such as Artuso et al. ( 2008 ) and Antonelli 



et al. (2010b) in the broader context of the low- 



energy /high^ntensity frontier. 

While by now kaon physics is a mature field, there 
are a number of challenges and exciting prospects for 
the near future. On the theoretical side the main chal- 
lenge will involve improving our control of strong inter- 
action effects. We anticipate that in the coming years 
lattice QCD (in synergy with CHPT) will play an in- 
creasing role in determining the long-distance dynamics 
relevant for kaon decays. This will not only increase the 
accuracy of predictions for semileptonic modes, but will 
also put the SM predictions for nonleptonic decays on 
firmer ground, in both cases enhancing the constraining 
power of these decays with respect to models of New 
Physics. On the experimental side, the world-wide pro- 
gram aimed at rare decays will face the challenge of mea- 
suring branching ratios at the level of 10~^^. This pro- 
gram will produce high-statistics results for other chan- 
nels as a by-product, and has a great potential to discover 
effects of New Physics in the "golden modes" K — >■ ttvv. 
Even in the event of agreement with the SM predictions, 
these measurements will provide essential information on 
the flavor structure of SM extensions at the TeV scale. 
In both scenarios, kaons will play a central role at the 
low-energy frontier of our search for physics beyond the 
SM. 
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TABLE VII Comparison between theory (D'Ambrosio et al. 19971 and experiment (Nakamura et al. 20101 for K — > 8717 
decays. The theoretical prediction for BK{Ks — >■ 7r^7r~7r"7) depends on NLO weak LECs and is a rough estimate only. 





B^^CMeV) 


BR (theory) 


BR (PDGTO) 


K+ -> 7V°7v\ + -f 


10 


(3.78 ±0.05) X 10-** 


(7.6 1^;'^) X 10"'' 


K+ ^7r+7r+7r-7 


5 


(1.26 ±0.01) X 10~* 


(1.04 ±0.31) X 10""' 


Kl — >■ TT TT^TT 7 


10 


(1.65 ±0.03) X 10-^ 




Ks — > 7r''"7r~7r"7 


10 


2 X 10-i« 
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Appendix A: One-loop functions 

The following functions occur in the one-loop ampli- 
tudes for rare and radiative K decays: 



{Airf 



1 — - arcsin 



z < 4 






\J '^1 z — 1 arcsin ( ^ J 



z < 4 



m \ z > 4 



mz) 



1 



2(1 -z)2 
-2z 



rlj Krl 



^'^)-<k 



(A2) 



, (A3) 



~C2o{p\p\M^,M^)], (A4) 

with p2 = M^ (for h^^K and hKrj) or p^ — M|- (for /itttt)- 
C20 is defined by the three-point integrals for q^ ~ 0; 



(2^)4 [P - Mf] [{£ + (7)2 - Mf] [{£ - p)2 - M|] 



^5^''C2ob^b±g)^Mf,M2) 



(A5) 
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